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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 
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BIVARIATE  LOCATION  PROBLEMS 

BY 

DAWN  PETERS 
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Chairman:  Ronald  H.  Randles 

Major  Department:  Statistics 

Affine  invariant  signed-rank  tests  are  proposed  for  the  one-  and 
two-sample  bivariate  location  problems.  The  test  statistics  proposed 
are  modifications  of  the  one-  and  two-sample  tests  of  Blumen  and 
Mardia,  respectively.  The  asymptotic  null  distributions  of  the 
proposed  statistics  are  obtained  for  the  class  of  elliptically 
symmetric  distributions.  In  addition,  the  asymptotic  distributions 
under  certain  contiguous  alternatives  are  obtained  for  the  proposed 
statistics,  as  well  as  for  Blumen’ s sign  statistic  and  a slightly 
modified  version  of  Mardia's  sign  statistic.  Comparisons  are  made 
between  the  proposed  statistics  and  several  competitors  via  Pitman 
asymptotic  relative  efficiencies  and  Monte  Carlo  results.  The  signed- 
rank  tests  proposed  appear  to  be  robust.  They  perform  better  than  the 
leading  competitors  when  the  underlying  distribution  is  light- 
tailed. For  heavy-tailed  distributions,  Blumen's  test  and  Mardia's 
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test  perform  better  than  the  proposed  one-  and  two-sample  statistics, 
respectively.  However,  the  proposed  statistics  clearly  dominate  the 
one-  and  two-sample  Hotelling’s  T in  both  light-tailed  and  heavy- 

p 

tailed  settings  and  perform  virtually  as  well  as  Hotelling's  T when 
the  underlying  population(s)  is(are)  bivariate  normal. 
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CHAPTER  1 
INTRODUCTION 

In  this  study  we  investigate  test  statistics  for  the  one-  and 

two-sample  bivariate  location  problems.  The  rank  procedures  we 

propose  are  modifications  of  existing  procedures  used  with  one  and  two 

samples.  For  the  one-sample  problem,  assume  Y,,...,Y  are 

_1  _n 

independently  and  identically  distributed  as  Y,  where  Y is  from  a 
continuous  bivariate  population.  Let  Y = (Y^.Y,^)',  1 = and 

0 = (QyQ^)'  = (median  (Y.^),  median  (Y^))'.  We  are  concerned  with 
the  test  of 

H : 0=0  against  H : 0*0.  (1.1.1) 

o _ „o  a _ _o 

Letting  = (Z.^,  )'  = (Y^  - 0q),  i = 1,...,n,  we  see  that  we  can, 

without  loss  in  generality,  consider  the  test  of 

H : 0=0  against  H : 0*0.  (1.1.2) 

o a 

Probably  the  most  well-known  procedure  to  test  the  hypotheses  of 
1.1.2  is  that  proposed  by  Hotelling.  Hotelling’s  T2  is  defined  as 

T2  = T2(Y1  , ...,Yn)  = n Y'lf1  Y, 

where  (1.1.3) 

- 1 n i n 

Y = - l Y. , and  I = — I (Y  -Y)(Y  -Y)'. 

~ n i=1  _i  n-1  i=1  -1  ~ -i  - 


1 


2 


If  l,  the  variance-covariance  matrix  of  Y,  is  positive  definite,  then 
2 

under  HQ  T is  asymptotically  Chi-square  with  two  degrees  of 

freedom.  If  Y is  a bivariate  normal  random  vector,  then  the  null 
2 

distribution  of  T is  a multiple  of  the  F distribution  with  two 
numerator  degrees  of  freedom  and  n-2  denominator  degrees  of  freedom. 
This  procedure  is  the  likelihood  ratio  test  of  1.1.2  when  the 
underlying  distribution  is  bivariate  normal.  Desirably,  T2  is 
invariant  with  respect  to  nonsingular  linear  transformations.  That 
is,  if  L is  a 2x2  nonsingular  matrix,  then 

T2(Y. Y ) = T2(LY. LY  ). 

- 1 _n  _ 1 _n 

For  the  one-sample  problem  we  shall  call  this  invariance  property 

affine-invariance.  (In  the  two-sample  setting  affine-invariance  will 

refer  to  statistics  which  are  invariant  with  respect  to  translations 

of  the  observations  in  addition  to  the  nonsingular  linear 

transformations  discussed  here.)  This  appealing  invariance  property 

ensures  that  rotations  and  reflections  of  the  observations  about  the 

origin  as  well  as  scale  changes  will  not  affect  the  test  statistic. 

Many  of  these  transformations  of  the  observations  are  equivalent  to 

changes  in  the  covariance  structure  of  the  population  sampled. 

2 

An  early  competitor  of  Hotelling's  T is  a bivariate  sign-test 
due  to  Hodges  (1955).  Hodges  counts  the  number  of  observed  vectors 
falling  on  each  side  of  a given  line  through  the  origin.  His 
statistic,  H,  is  the  maximum  of  the  counts  obtained,  maximizing 
over  all  lines  through  the  origin.  H is  distribution  free  under 


3 


HQ  for  the  class  of  continuous  distributions  such  that  0 satisfies 

P(aY1+bY2  > a01  + b02)  = 1/2,  (1.1.4) 

for  all  (a,b)  * (0,0).  Additionally,  H is  invariant  with  respect  to 

nonsingular  linear  transformations.  The  null  distribution  of  H was 

tabulated  by  Hodges  for  small  samples.  Later,  Joffe  and  Klotz  (1962) 

simplified  Hodges’  representation  of  the  null  distribution  yielding  an 

expression  for  the  asymptotic  null  distribution  as  well. 

Another  affine-invariant  statistic  was  proposed  by  Blumen 

(1958).  His  sign  test  statistic  is  obtained  by  conditioning  on  the 

vectors  Y^ , . . . ,Y  and  their  respective  projections  through  the  origin, 

namely  -Y^,...,-Y  . The  angles  between  adjacent  vectors  are  adjusted 

so  that  each  of  these  standardized  angles  is  ir/n  radians.  The 

distances  of  the  observed  points  from  the  origin  are  also  standardized 

so  that  each  vector  has  unit  length.  This  standardization  process  is 

illustrated  in  Figures  1.1  through  1.3.  The  three  vectors  Y.)  , Y , 

and  Y^  shown  in  Figure  1.1  are  displayed  again  in  Figure  1.2  along 

with  their  projections  through  the  origin,  -Y1 , and  -Y  . In 

Figure  1.3  the  angles  between  adjacent  vectors  are  standardized  as  are 

the  lengths  of  each  vector.  The  standardized  vectors  are  denoted 
s s s 

by  Y^  , Y2  , and  Y^  . This  standardization  method  is  used  to  preserve 

affine-invariance.  (This  fact  will  be  established  clearly  in  Section 

2 of  Chapter  2.)  To  obtain  Blumen' s statistic,  the  center  of  gravity 

of  the  standardized  versions  of  the  vectors  Y, Y is  found.  His 

~ i ~n 

test  statistic  is  proportional  to  the  squared  distance  from  the  origin 
to  the  center  of  gravity.  The  further  the  center  of  gravity  is  from 


Figure  1. 


4 


VS 


Figure  1.1  Observed  Vectors 


2 Observed  Vectors  and  Their  Projections  through  the  Origin 
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Figure  1.3  Standardized  Vectors 
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the  origin,  the  greater  the  evidence  that  H is  true.  Blumen’ s 

a 

statistic  is  distribution  free  under  H for  the  class  of  distributions 

o 

given  by  1.1. 4.  Blumen  shows  that  the  asymptotic  null  distribution  of 
his  statistic  is  Chi-square  with  two  degrees  of  freedom.  The  null 
distribution  of  Blumen' s statistic  was  tabulated  by  Klotz  (1964)  for 
small  samples.  Our  proposed  signed-rank,  affine-invariant  statistic 
is  a modification  of  Blumen' s sign  statistic.  We  will  present  a more 
formal  discussion  of  Blumen' s statistic  in  the  following  chapter  as  an 
aid  to  the  introduction  of  our  modification.  Subsequently,  additional 
properties  of  Blumen' s test  will  be  derived. 

Several  one-sample  statistics  have  been  proposed  which  are  not 
affine-invariant.  Bennett  proposes  two  procedures.  The  first  method 
(Bennett,  1962)  considers  the  signs  of  the  components  of  (Z.^,  Z^.), 
i = 1 , . . . ,n.  Letting 


S11  ' I[Z1i  > °'  Z2i  > 0]' 


£10  ‘ j,  I[ZU  > Z2i  < «• 


S01  ■ Z Itz,l  < Z2i  > °^' 

1 = 1 


and 


S 


00 


n 

Z ICZ  < 0,  Z < 0], 
i = 1 1 


Bennett's  first  statistic  is  given  by 


(1.1.5) 


The  asymptotic  null  distribution  of  B1  is  Chi-square  with  two  degrees 
of  freedom  and  hence  B1  can  be  used  for  a large  sample  test. 

Bennett's  procedure  was  introduced  as  a general  p-variate  procedure, 
P-2.  His  suggested  bivariate  procedure,  given  here,  is  identical  to 
the  large  sample  test  procedure  proposed  by  Chatter jee  (1966). 
Chatterjee  also  develops  a small  sample  conditionally  distribution- 
free  procedure  based  on  ^ , S^,  SQ1  , and  S . To  do  this, 
Chatterjee  uses  the  result  that  conditional  on  S„  +S„„  and  S +S  , 

11  00  10  or 

and  are  independent  and  distributed  as  binomial  random 
variables  with  parameters  (S^+Sqq,  1/2)  and  (S1  +SQ1,  1/2), 
respectively.  The  second  procedure  proposed  by  Bennett  ( 1 96 i4 ) is  a 
signed-rank  test  generalizing  Wilcoxon's  univariate  signed-rank 
test.  Bennett  bases  his  test  on  the  univariate  signed-ranks  of 


{Z 


1 1 


V a"d  fZ21 


9 • • • 9 


Z^}.  Let 


U+  = 1 Rank(|Z11|)I[Z11  > 0]  - j n(n+1) 


V+  - Z Rank  ( | Z | )I[Z  > 0]  - -j  n (n+1 ) , 
i-1 


and 


(1.1.6) 
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where  Rank(|z_|)  denotes  the  rank  of  |Z  | among  Iz.J lZinl» 

i = 1.  2.  Bennett's  signed-rank  statistic  is  a quadratic  form  based 
+ + 

on  U and  V . His  statistic  is  asymptotically  Chi-square,  under  HQ, 

with  two  degrees  of  freedom.  In  addition  to  the  quadratic  forms 

proposed  by  Bennett  and  Chatterjee,  Bickel  (1965)  and  Sen  and  Puri 

(1967)  also  consider  quadratic  functions  of  both  the  coordinate-wise 

sign  statistics  and  the  coordinate-wise  signed-rank  statistics.  The 

class  of  statistics  examined  by  Sen  and  Puri  is  the  most  general, 

allowing  functions  of  the  coordinate-wise  ranks  to  be  used  in  place  of 

the  ranks  themselves.  Both  Bickel  (1965)  and  Sen  and  Puri  (1967) 

present  asymptotic  properties  of  the  statistics  considered.  None  of 

the  quadratic  forms  discussed  are  affine-invariant.  Hence,  the 

performances  of  these  procedures  vary  with  changes  in  the  underlying 

covariance  structure  of  the  population  and  with  changes  in  the 

direction  of  the  shift  of  0 from  the  origin. 

Another  generalization  of  the  univariate  Wilcoxon  signed-rank 

statistic  is  given  by  Killeen  and  Hettmansperger  (1972).  They 

consider  the  perpendicular  projection  of  Z.,  i = 1,...,n,  onto  the 

line,  l(t),  which  is  a counterclockwise  rotation  of  t radians  of  the 

horizontal  axis.  Letting  Z1  (t) , . . . , Z (t)  be  the  projected  points, 

Killeen  and  Hettmansperger  then  calculate  the  Wilcoxon  signed-rank 

statistic,  W(t ) , based  on  the  relative  positions  of  Z.(t) Z (t) 

1 n 

on  the  line,  l(t),  whose  center  is  at  the  origin.  The  test  statistic 

suggested  is  max  W(t).  Although  this  statistic  has  high  Bahadur 
0<t<2ir 

efficiency,  neither  its  small  sample  nor  its  asymptotic  null 
distributions  have  been  found,  and  hence  its  utility  is  questionable. 
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Brown  and  Hettmansperger  (1985)  propose  a statistic  based  on  the 
gradient  of  Oja's  measure  of  scatter  (Oja,  1983).  Letting 
A(Y^,  Yj,  0)  denote  the  area  of  the  triangle  formed  by  Y^ , Y ^ , and 
9 as  vertices,  define 


n n 

T(0 ) - Z I A(Y  , Y , 0).  (1.1.7) 

i-1  J-1  ~ ~ 

i*J 


The  value  0 that  minimizes  T(0)  is  Oja's  generalized  median  for 
bivariate  samples.  Brown  and  Hettmansperger  propose  the  use  of 


Qn(0)  - VT(9), 


(1.1.8) 


the  vector  of  partial  derivatives  of  T(0).  The  generalized  median  of 

A A 

Oja  is  the  value  0 such  that  |Qn(0)|  is  minimized.  The  proposed  large 
sample  test  statistic  for  the  hypotheses  of  1.1.2  is  a quadratic  form 
of 


S 


n 

E 

i-1 


), 


(1.1.9) 


where  Q (-Y  ) is  calculated  using  the  observations  Y Y , and 

2n  _i  _1  n 

their  projections  through  the  origin,  -Y, ,...,-Y  . 

_1  _n 

Also  invariant  with  respect  to  linear  transformations  is  a 
statistic  proposed  by  Dietz  (1982).  She  suggests  using  the  sign  and 
signed-rank  procedures  proposed  by  Bickel  after  a transf ormation  is 
applied  to  obtain  affine- invariance . Each  of  the  tests  proposed  by 
Brown  and  Hettmansperger  (1985)  and  Dietz  (1982)  are  quite  difficult 


10 


to  apply.  The  asymptotic  efficiencies  of  their  procedures  have  not 
been  determined. 

We  now  consider  bivariate  tests  for  location  in  the  two-sample 
setting.  Many  of  the  procedures  proposed  are  natural  extensions  of 
the  one-sample  procedures  discussed  earlier  in  this  section.  Assume 
^1  ’ * * ’ * jn  an<^  ^ , . . . , Y are  independent  random  samples  from  continuous 
bivariate  populations  with  distribution  functions  G1 (x)  and 
G^(y)  = G (y-0)  , respectively.  We  shall  consider  tests  of 

H : 0=0  against  H : 0*0.  (1.1.10) 

O M 3 ^ w 

The  normal  theory  affine-invariant  test,  extending  the  one-sample 
2 

Hotelling’s  T , takes  the  form 

= 2H2(X-Y)'S_1(X-Y), 

where  (1.1.11) 

- 1 m - 1 n 
N = m+n,  X = - I X.,  Y=-  E Y , 

- m . , _i  n _j’ 


and 


. m n 

S = — { E (X  -X)(X  -X)'  + E (Y  -Y) (Y  -Y) . 
i=1  ~ ~ -1  ~ j=i  ~ ~ 

2 

The  asymptotic  null  distribution  of  T^  is  Chi-square  with  two  degrees 
of  freedom.  When  both  populations  are  bivariate  normal,  a small 
sample  test  can  be  obtained  using  the  exact  null  distribution  of 


11 


, which  is  a multiple  of  the  F distribution  with  two  numerator 
degrees  and  N~3  denominator  degrees  of  freedom. 

A generalization  of  Blumen's  one-sample  statistic  is  given  by 
Mardia  (1967).  Mardia's  affine-invariant  sign  statistic  uses  the 
centered  sample  observations, 


X1  Z,..«, Xm — Z , an d Y ^ — Z , . . . , Y^ — Z , 


where  Z refers  to  the  arithmetic  mean  of  all  m+n  = N vectors.  The 
standardization  process  of  Blumen's  is  applied  to  the  N centered 
vectors.  The  center  of  gravity  of  the  standardized  versions  of  the 
vectors  from  the  first  sample  is  then  obtained.  Again,  the  squared 
distance  from  the  center  of  gravity  to  the  origin  is  the  basis  of  the 
test  statistic  used.  The  precise  definition  of  Mardia's  statistic  is 
given  in  Chapter  5.  Our  proposed  two-sample  signed-rank  statistic 
will  be  seen  to  be  a modification  of  Mardia's  statistic.  It  is  a 
natural  extension  of  our  one-sample  procedure  to  two  samples. 

A distribution-free  procedure  designed  to  detect  broad 
alternatives  was  proposed  by  Vincze  (1961).  Vincze  suggests 
projecting  each  vector  onto  a randomly  selected  line  through  the 
origin.  The  Kolraogorov-Smirnov  two-sample  test  statistic  is  then 
applied  to  the  projected  vectors  using  their  relative  positions  on  the 
line  selected.  As  a precursor  to  Killeen  and  Hettmansperger ' s one- 
sample  procedure,  Vincze  also  suggests  that  it  may  be  desirable  to 
maximize  the  Kolmogorov-Smirnov  deviations  of  the  projected 
observations  over  all  lines  through  the  origin.  However,  the  null 
distribution  when  maximizing  has  not  been  found. 
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Chatterjee  and  Sen  (19614)  extend  the  Wilcoxon  Rank  Sum  Test  and 
Mood's  Median  Test  to  the  bivariate  setting.  The  rank  sum  test  uses 
the  centered  Wilcoxon  Rank  Sums  of  the  marginal  variates: 

m N+1 

S - Z Rank (X  ) - 

J-1  J 


and 


(1.1.12) 


m 


,N+1 


S = E Rank (X  ) - m(^), 
j-1  J 


where  RanktX^)  denotes  the  rank  of  , among 

{Xil , . . . ,Xim,  Y11 , . . . »Yin) , i - 1 , 2.  The  large  sample  test  statistic 
is  a quadratic  form  based  on  S^  and  S^.  A small  sample  conditional 
test  based  on  S^  and  S^  is  also  described.  This  test  is  also  examined 
by  David  and  Fix  (1961)  in  a discussion  of  the  comparison  of  two 
regression  lines. 

To  describe  the  extension  of  Mood's  median  test,  let 


M1  = median  of  {Xn  , . ..  ,Xlm>  Y^} 


and 


M2  = median  of  {X21 X^,  Y21  ’ * ’ * ’ Y2n}  * 


Define 


m 

"l 1 • I[XU  s V x21  s V 
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m 

"o,  ■ l I[xu  > M1  ’ x2i  s V 


m 

M10  • I[x,l  SV  X2i  > V 


M00  ' j,  I[XU  > V X21  > V 


and 


m 


C‘  ^ I[X,i  SV  X21  S V * jf,  I[I,j  SM1-  X2J  SM2]- 


The  bivariate  median  test  statistic  of  Chatterjee  and  Sen  is  given  by 


CS, 


N 

2mn 


(M  “M-. )‘ 
11  00 


(M10~H01  l 
n||]-c 


where  [|x|]  denotes  the  greatest  integer  less  than  or  equal  to  x.  A 
conditional  small  sample  test  and  a large  sample  Chi-square  test  are 
discussed . 

Suguira  (1965)  and  Bhapkar  (1966)  consider  multivariate 
multisample  extensions  of  the  Wilcoxon  Rank  Sum  Statistic.  In  the 
bivariate  two-sample  setting,  these  large  sample  tests  are  functions 
of  the  component-wise  Wilcoxon  Rank  Sum  Statistics.  Puri  and  Sen 
(1966)  and  Tamura  (1966)  discuss  a more  general  class  of  Chernoff- 
Savage  type  statistics  for  the  multivariate  multisample  location  and 
scale  problem.  As  in  the  one-sample  setting,  score  functions  of  the 


coordinate-wise  ranks  may  be  included.  Both  articles  discuss 
asymptotic  properties  of  the  given  class. 

Finally  we  mention  that  Brown  and  Hettmansperger  (1985)  have  also 
discussed  a two-sample  test  based  on  Oja's  measure  of  scatter.  The 
large  sample  test  statistic  used  is  a quadratic  function  of 

m 

?2  ’ VV’ 

where  QN  is  the  gradient  function  defined  by  1.1.7  and  1.1.8  as 
applied  to  the  combined  sample  of  N observations. 

With  the  exception  of  Hotelling’s  T ^ , Mardia's  test  statistic, 
the  test  statistic  of  Brown  and  Hettmansperger , and  the  two-sample 
procedure  we  shall  propose,  the  other  procedures  mentioned  are  not 
affine- invariant . 

The  one-sample  statistic  we  propose  is  defined  in  Chapter  2.  Its 
invariance  with  respect  to  nonsingular  linear  transf ormations  is  also 
verified.  In  Chapter  3 we  derive  the  asymptotic  null  distribution  of 
the  one-sample  test.  The  asymptotic  distributions  under  Pitman 
alternatives  of  our  test  and  the  test  of  Blumen  are  then  determined. 
Comparisons  of  several  one-sample  procedures  are  made  in  Chapter  4 via 
Pitman  asymptotic  relative  efficiencies  and  Monte  Carlo  results.  The 
introduction  of  our  two-sample  test  and  the  verification  of  its 
invariance  with  respect  to  nonsingular  linear  transformations  is  given 
in  Chapter  5.  Chapter  6 presents  the  asymptotic  distribution  of  the 
proposed  two-sample  test  statistic  under  the  null  hypothesis  and  under 
a sequence  of  Pitman  alternatives.  Also  given  are  the  corresponding 
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asymptotic  distributions  of  a "centered"  version  of  Mardia's  two- 
sample  test.  Finally  several  two-sample  test  statistics  are  compared 
in  Chapter  7 through  asymptotic  relative  efficiencies  and  some  Monte 
Carlo  studies. 


CHAPTER  2 

A BIVARIATE  ONE-SAMPLE  TEST  FOR  LOCATION 


2.1  Definition  of  the  Test  Statistic 

In  order  to  motivate  the  test  statistic  we  are  proposing  we  will 

first  describe,  more  formally  than  in  Chapter  1,  the  test  statistic 

proposed  by  Blumen . We  assume  that  YlP...,Y  are  independent  and 

identically  distributed  as  Y,  where  Y is  from  a continuous  bivariate 

distribution  with  median  9 satisfying  1.1.4.  Let  a' a'  denote  the 

1 rn 

angles  between  the  positive  half  of  the  horizontal  axis  and 
Y1 Yn.  respectively.  Here  0 S < 2it,  1 = 1 n.  Let 


and 


4>',  if  OS  4>'  < ir, 
if  it  < $'  < 2tt, 


1 , 

if 

VI 

O 

*i' 

$ IT, 

-1, 

if 

TT  < 

♦i' 

t= 

OsJ 

V 

(2.1.1  ) 


for  i = 1,...,n.  For  each  i,  let  RankC^)  denote  the  rank  of  <j> 
among  * • • • » <J>n • (That  is,  Rank(<f>  ) equals  the  number  of  angles 
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4>1 , . . . , 4>n  that  are  less  than  or  equal  to  <j>  .)  If  we  denote  the  center 
of  gravity  of  the  standardized  observations,  as  described  in  Chapter 
1.  by  (C  , S ),  then  the  center  of  gravity  can  be  defined  by  the 

D D 

following: 


1 


n 


TrRank(<j>.  ) 


C = - E a . cos  (- 
B n i-1  1 


±-> 


and 


n irRank(<|>.) 

Sn  = — E a.  sin( — ) 

B n . . i n 

l=i 


Blumen’s  sign  statistic  is  given  by 


B = 2n (Cg  + Sg) 


(2.1 .2) 


(2.1.3) 


With  Blumen’s  statistic  as  background,  we  can  now  introduce  our 
proposed  signed-rank  statistic.  We  begin  as  Blumen  did  by  considering 
the  observed  vectors  and  their  projections  through  the  origin.  The 
angles  between  adjacent  vectors  are  standardized  as  in  Blumen' s 
statistic,  but  the  lengths  of  the  vectors  are  not.  Instead,  we  use  a 
measure  of  the  length  of  each  vector  that  is  invariant  with  respect  to 
linear  transformations  of  the  observations.  Let 


Y'  E_1  Y 
ti  L ii 


for  i = 1 , . . . ,n  , 


where 


(2.1.M) 


l - — 


n 

Z 

i-1 


Y Y ' 

:i:i* 


n 
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The  matrix  £ is  a consistent  estimation  of  E[YY'],  the  covariance 

matrix  of  Y under  H0>  provided  Y is  symmetrically  distributed  about 

1 1*1 

its  median.  We  now  weight  the  i terms  in  the  averages  CD  and  SD  of 

B B 

Blumen's  statistic  by  the  rank  of  among  D and  define 


(CMB*  SMB 


) by 


and 


'MB 


n T;Rank(4). ) 

- I a cost — )Rank(D.) 

n i=1  1 


1 


n 


MB 


= - l 
n i = 1 


a^sint 


irRank  ( ^ ) 


)Rank  ) 


(2.1.5) 


Our  modification  of  Blumen's  statistic,  which  we  shall  denote  by  MB, 
is  given  by 

HB  ’ 6n_1  IC„*4>  ■ <2- 1-6) 


2.2  Affine-Invariance  of  MB 

In  this  section  we  will  verify  that  our  proposed  statistic  is 

invariant  with  respect  to  nonsingular  linear  transformations  of  the 

observations.  If  L is  a 2x2  nonsingular  matrix,  we  will  show 

that  MB( Y Y ) = MB( LY LY  ),  where  MB(Z, Z ) denotes  the 

-l  ~n  _l  _n  _1  _n 

proposed  statistic,  MB,  applied  to  the  observations  Z„,...,Z  . This 

_1  _n 

result  is  established  via  a sequence  of  three  lemmas: 

Lemma  2.2.1.  MB  is  invariant  with  respect  to  a rotation  of  the 


observations . 
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Proof.  Let  4>  ^ , . . . , be,  as  defined  in  2.1.1,  the  angles  formed 

by  the  vectors  projected  above  the  horizontal  axis  with  respect  to  the 

positive  half  of  the  horizontal  axis.  For  simplicity,  let  us  call 

♦ . <t>  the  projected  angles  associated  with  Y ,...,Y  , 

1 n _ 1 _n 

respectively.  Now  consider  a clockwise  rotation  of  the  observations 

such  that  a single  observed  vector  crosses  the  horizontal  axis.  (See 

Figures  2.1  and  2.2.)  Let  Y^,...,Y  be  the  vectors  corresponding 

# * 

to  Y i » • • • , Y^ , respectively,  after  the  rotation.  Let  4>1  , . . . , be  the 

£ * 

projected  angles  associated  with  Y1f...,Y  , respectively.  As  can  be 

seen  in  Figures  2.3  and  2.4,  the  relationship  between  the  rank  of 

* 

4>i  and  the  rank  of  $ is  given  by 


* 

Rankt^)  = 


if  Rank(<()i)=1  , 


iRank^  )-1 , 


otherwise. 


Since  a rotation  of  the  vectors  Y^,...,Y  can  be  represented  as 
RY1,,,',RIn’  for  some  Positive  definite  2x2  matrix  R,  and  since 


n 

E (RY  ) (RY . ) ' = R E R' 
i=1  ~ -1 


(RYi)'(RZR')“1 (RYi)  = Yj'  E_1  Y = D , 


and 


(2.2.2) 
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Figure  2.1  Observed  Vectors 


Figure  2.2  Rotated  Vectors 
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Figure  2.3  Projected  Angles  of  Observed  Vectors 


Figure  2.H 


Projected  Angles  of  Rotated  Vectors 
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for  i = we  see  that  the  standardized  distances,  D , are 

1 

unaffected  by  such  a rotation.  (See  also  Theorem  2.2.7.)  Let 
* * 

a. a be  the  sign  functions  associated  with  the  rotated 

I n 

# * 

observations,  Y^,...,Y  , respectively,  defined  as  in  2.1.1.  The  only 
sign  function  which  is  changed  through  a rotation  is  that  associated 
with  the  vector  which  is  rotated  across  the  horizontal  axis.  That  is, 


* 


a 


i 


if  Rank(<|> ) =1  , 
otherwise . 


Given  the  above  notation,  our  signed-rank  statistic  when  applied 

# # 


to  the  rotated  observations,  Y Y becomes 

1 n 


-1  *2  *2 

MB  = 6n  (C  +S  ) , 


where 


and 


* 

* 1 n # TrRank(<f>.) 

C - — I a cos( — )Rank(D  ), 

n l n i 


* 

S 


* 

1 n # irRank(i}>.) 

- Z a sin( — — )Rank (D  ) . 

* i ^ = ^ i n i 


(2.2.3) 


If  the  kth  vector  is  the  one  which  is  rotated  across  the  horizontal 
axis,  that  is  Rank(<f>k)  = 1,  then  using  the  relationships  described 


above  we  can  write 
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# 1 n -it (Rank  ( 4>  )- 1 ) 

C = - I aicos(  - )Rank(D^)  + — (-a^cosC^-)  )Rank  (Dk) . 

i*k 

Noting  that 


cos(tt)  = -cos(O) 


tt  (Ran k ( 4>  )-1) 
-cos( - ) , 


we  have 


* 

C 


n 


n 

I 

i = 1 


a^cos ( 


ir(Rank(4)  )-1 ) 

)Rank(Di  ) 


and,  likewise, 


» 

S 


n 


n 

I 

i = 1 


irCRankU  )-1 ) 

a.sin( )Rank (D,  ). 

ini 


We  will  now  make  use  of  two  elementary  trigonometric  identities: 


cos(b-d)  = cos(b)cos(d)  + sin(b)sin(d) 


and 


(2.2.4) 


sin(b-d)  = sin(b)cos(d)  - cos  (b  )sin  (d ) . 


irRank  (<f>  ) 

Taking  b = and  d 

n 


j-j-,  we  can  write 


* 

C 


. n irRank(<J>  ) 

- E a.  {cos( )cos (— ) + 

n 1=1  i n n 


irRank(<j>  ) 

sin( )sin(-)  }Rank  (D,  ) 

n n i 


and 
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* 1 n irRank  (<J>  ) irRank  ( <f>  ) 

s =n  1 ai{sin( )cos(J)  - cos( — )sin(-)  }Rank  (D.  ). 

“ i_i  ri  n n n 1 


That  is, 


and 


°0S(S,CMB  * Sln(S)SMB 


S ' C03(S)SMB  ' Sln(iT)CMB' 


where  and  SMB  are  defined  by  2.1.5.  Hence 


c 2+S  2 = (c°s(^)CMB  + sin(J)SMn)2  + (cos&S^  - sin(^)CwrJ2 


n MB' 


n'  MB 


n MB' 


• °°32#<cmb*s«b>  * si"2<£>(<4*4> 


2 2 
°MB  SMB* 


Clearly,  further  clockwise  rotations  of  the  observations  will  not 

affect  the  value  of  MB.  Since  any  counterclockwise  rotation  of  the 

observations  can  be  obtained  by  an  appropriately  chosen  clockwise 

rotation,  we  see  that  the  invariance  of  MB  with  respect  to  clockwise 

or  counterclockwise  rotations  of  the  observations  is  verified.  □ 

Lemma  2.2.5.  MB  is  invariant  with  respect  to  reflections  of  the 

observations  about  the  horizontal  axis. 

Proof.  A reflection  about  the  horizontal  axis  results  in  a 

reversal  of  the  ordering  of  the  projected  angles.  If  a* / denote 

1 Yn 
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the  projected  angles  of  the  reflected  vectors,  and  4>  ,...,4)  denote 

1 n 

again  the  projected  angles  of  the  original  vectors,  then,  for  all  i, 

* 

Rank(<j>^ ) = n + 1 - Rank(<J>  ) . 

Although  the  standardized  distances,  D^,  are  not  altered,  it  is  easily 
seen  that  each  of  the  sign  functions  change.  The  modification  of 
Blumen's  statistic  applied  now  to  the  reflected  vectors  is 

-1  *2  *2 

6n  (C  +S  ) , where  here 

* 1 n ir(n+1  -Rank(<f>.  ) ) 

C = - E -a  cos( — )Rank(D.) 

^ ^ n 1 


and 


* 1 n irCn+l-RankC^) ) 


S = - I -a,sin( 
i = 1 


n )Rank(D  ). 


Again  making  use  of  the  trigonometric  identities  given  in  2.2.H,  we 
write 


* 1 n irfn+1')  irRank(<f>  ) TrRank(4>.) 

c ■ "SI^alIoos(_Sr^)cos( S )+3ln(— )sln( i-JIRanktDp 


. -{ oo3 ( 1*2111) Cmd  . sln(i(n^l,s  , 


n MB 


n MB ' 


and 
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S 


* 


ir(n+1  ) 


n 


so  that 


Lemma  2.2.6.  MB  is  invariant  with  respect  to  changes  in  scale. 
Proof.  Consider,  without  loss  in  generality,  multiplying  the 
first  element  of  each  bivariate  observation  by  the  positive  constant, 
b.  To  produce  this  scale  change  we  multiply  each  observation  by  the 
2x2  positive  definite  matrix. 


unaffected.  It  is  clear  that  the  sign  functions  are  also  unaffected 
as  each  observation  cannot  be  transformed  across  the  horizontal 
axis.  The  ordering  of  the  angles  is  unaffected  as  well,  as  can  be 
seen  by  considering  the  effect  of  a change  in  scale  on  the  slopes  of 
the  observed  vectors.  The  slopes  of  the  transformed  vectors  will  all 
be  a constant  positive  multiple  of  the  slopes  of  the  associated 
original  vectors.  Trivially,  MB  is  thus  unaffected.  □ 

Theorem  2.2.7.  MB  is  invariant  with  respect  to  nonsingular 


b 


S = 


0 


As  illustrated  by  2.2.2,  the  standardized  distances,  D , remain 

l 


linear  transformations  of  the  observations. 
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Proof.  Let  us  represent  a nonsingular  linear  transformation  by  a 

2x2  nonsingular  matrix,  L.  This  linear  transf ormation  applied 

to  Y ,...,Y  results  in  the  transformed  observations  LY  ..  LY 
~n  _1  ’ _n* 

Using  the  singular  value  decomposition  theorem  (see,  e.g.,  Rao,  1973, 
p.  42),  we  can  represent  L as  PAQ',  where  P and  Q are  2x2  orthogonal 
matrices  and  A is  a 2x2  diagonal  matrix  with  positive  diagonal 


elements.  Hence,  it  suffices  to  show  that  MB  is  invariant  with 

respect  to  the  transformations  defined  by  P,  A,  and  Q.  That  is,  we 

need  to  show  that  MB(Y  ,...,Y  ) = MB(PY PY  ) = MB ( AY AY  ), 

- 1 ~n  -1  ~n  _ 1 _n 


for  any  2x2  orthogonal  matrix  P and  any  2x2  diagonal  matrix  A with 
positive  diagonal  elements.  Firstly,  note  that  the  transformation,  A, 
produces  a change  in  scale  as  discussed  in  Lemma  2.2.6.  It  remains, 
therefore,  to  consider  only  orthogonal  transformations.  A 2x2 
orthogonal  matrix,  P,  can  be  represented  by  either 


P1  " 


or 


cosg 

sing 

-sing 

cosB  y 

cosg 

sing 

sing 

-cosg 

for  some  6 between  0 and  2ir.  It  is  easily  established  that  P^ Y 
is  a clockwise  rotation  of  Y through  an  angle  of  B radians.  Upon 
noting  that 
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P 


2 


0 


-1 


0 


-sinB  cosB 


cosg 


sinB 


we  see  that  P^Y  involves  a clockwise  rotation  of  Y through  an  angle  of 
6 followed  by  a reflection  about  the  horizontal  axis.  Hence,  by 
Lemmas  2.2.1  and  2.2.5,  the  invariance  of  MB  with  respect  to 
nonsingular  linear  transformations  is  established.  □ 

The  effect  of  a nonsingular  linear  transformation,  on  the 
observations  themselves,  has  already  been  illustrated  through  the 
three  lemmas.  It  may  be  worth  noting  again  that  such  a transformation 
amounts  to  a change  in  the  variance-covariance  structure  of  the 
distribution  of  the  observed  random  vectors.  Because  of  the 
invariance  property  of  MB,  we  will  later  be  able  to  restrict  attention 
to  simplified  variance-covariance  structures.  Note  finally  that  the 
preceding  results  apply  to  Blumen's  statistic  as  well.  That  is,  we 
have  simultaneously  demonstrated  the  invariance  of  Blumen's  statistic 
to  nonsingular  linear  transformations. 


CHAPTER  3 

ASYMPTOTICS  FOR  THE  ONE -SAMPLE  PROBLEM 

3.1  Asymptotic  Null  Distribution  of  MB 
In  this  section  we  will  find  the  asymptotic  null  distribution  of 
MB  under  a class  of  elliptically  symmetric  distributions.  We  first 
express  MB  in  a more  convenient  fashion.  Let 

1 n 

Hn(t)  = - I I(<t>i  £ t) 


and 


Fn(t)  = £ I I(D  S t) 


represent  the  empirical  distribution  functions  of  the  projected  angles 

4>1 » • • • » <t>n  and  the  standardized  distances  D1f...,D  , respectively.  We 

can  rewrite  the  expressions  for  C and  given  in  2.1.5  as 

Mb  MB 

n 

C = E a cos(ttH  U . ))F  (D  ) 
mb  i n i n i 


and 


(3.1.1) 


n 

S = E a sin (irH  (*  ))F  (D  ). 
MB  , . i n i n i 
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In  finding  the  asymptotic  null  distribution  of  MB  = 6n_1 (C2  +S2  ),  we 

MB  MB 

construct  a random  vector  which  has  the  same  asymptotic  null 
distribution  as  n',/2(C1||1  SMB>.  Let,  for  1 - 1 n. 


£o  ' Ethh']> 


and 


* -1 

D = y ' T.  Y 
i _i  o :i’ 


H(t)  = PH  (*i  * t)f 
o 


F (t ) = PH  (Di  £ t). 
o 


Define 


CMB  ” 1 aicos(1THU<  ))F(D, ) 


i = 1 


i"  ' i‘ 


and 


(3.1.2) 


n 

- l aisin(irH(^1))F(D1). 


It  will  be  shown  that  n'172^,  S^B)  and  n"1/2(CMB>  S^)  have  the 
same  asymptotic  null  distribution  under  certain  conditions.  To  verify 
this  we  will  need  to  show  that,  for  all  i, 


H 


lHn(V  * H ( <|>i  ) | = o(1) 


Fn(Di)  - F (D 


H 

0 

■ °p(1>- 


and 
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These  conditions  will  be  shown  to  hold  using  the  following  two 
lemmas . 


Lemma  3.1.3.  Suppose  X X are  independent  random 

nl  nn 

variables  with  distribution  functions  F ...... F , respectively,  and 

nl  nn  K J ’ 

1 n 

empirical  distribution  function  F (t)  = - I I [ X . < t].  Let 


Fn(t)  = n Z Fni(t)*  then 
i = 1 


n 


i = 1 


Slip  | F ( t ) - F ( t ) | ■*  0 a.s.,  as  n ■*  ». 

*-®<t<°° 


Proof.  See  Theorem  1 of  Shorack  and  Wellner  (1986,  p.  106).  (If 
Fn1  “ •••  “ Fnn  = F»  then  this  result  reduces  to  the  more  familiar 
result  of  Glivenko  and  Cantelli.  See,  e.g.,  Serfling,  1980,  p.  61.) 

Lemma  3.1.^.  Let  X^,...,X  be  independent  and  identically 
distributed  as  X,  where  X has  distribution  function  F 

a, 

a = (o^ a )'.  Assume  g(X,a)  is  a real  valued  function  with 

P(g (X , a)  £ t)  = G (t)  and  that  G (t)  is  continuous  in  t.  Let 
a a 

^ A A ^ ^ 

<*n  = (on1 , . . . ,a  ) ' be  an  estimator  of  a.  Suppose  that 
(1)  n^^(an~a)  » 0^(1),  as  n •>  », 

and  (3.1.5) 

(ii)  sup | g (X ,a+n  1/2s)  - g(X,a)|  > o (1),  as  n + ■, 
seg  p 

for  any  bounded  sphere  8 in  RP.  Let 


G .(t)  = r 


n ,a 
_n 


n 

I 

i=1 


I[g(x. 


“n> 


t] 
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denote  the  empirical  distribution  function  of  the  g(X. , a )'s.  Then 

i _n 

for  every  positive  integer  i, 

G ~ (g(XlP  an))  - G (g(X  a))  = o (1),  (3.1.6) 

n ,a  ~ P 

_n 

as  n ■*  ». 

Proof.  First  we  will  show  that  conditions  (i)  and  (ii)  imply 
g(x»«n)  ■ g(x»  <*)  = op(1).  Note  bx  (i)  that  for  any  6 > 0,  there 
exists  a bounded  sphere,  B,  such  that 

1 /2  ~ 

P(n  (an-a)  i 6)  < 6,  for  all  n.  (3.1.7) 

Let  6 > 0,  6'  > 0,  and  e > 0 be  arbitrary.  Write 

P(|g(X,an)  - g(X,a)|  > e) 

= P(|g(x*a  ) - g(X,a)|  > e,  n^2(a  -a)  e 8) 

~n  - _n  _ 

+ P(  | g (X  ,ot  ) - g(X,a)|  > e,  n1/2(a  -a)  t 6) 

~n  ~ _n  _ 

^ P(|g(X,a  ) - g(X,a)|  > e,  n^2(a  -a)  e B) 

~n  - ~n  _ 


+ 6 . 


Letting 


-1  /? 

Wn(s)  = g(X,  a+n  s)  - g(X,  a) 


so  that 
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W (n 
n 


1 /2 


g(X,  otn)  - g(X,  a)  , 


we  have 

P(|g(X,an)  - g(X , a) | > e) 

£ P( I wn (n 1 /2(an-a) ) | > e,  n1/2(an-a)  e (J)  + 6 

£ P(sup |w  (s) I > e)  + 6 
seg 

£ 6 ' + 5 , for  large  n , 

by  condition  (ii).  Using  this  result,  we  will  complete  the  proof. 

Let  6 > 0 and  e > 0 be  arbitrary.  Write 

A 

|G  . (g (Xj , an))  - G (g(X  , a) ) | 
n*«n 

s lG  - (g(X1,  an))  - G (g(X  , a ))| 

"*“n 

a 

+ lGa(g(V  2n})  “ Ga(g(Xi*  «})l • 

A 

Using  the  fact  that  g(X,  a ) - g(X,a)  = o (1)  and  G (•)  is  continuous, 

~n  _ p a 

we  see  that  the  second  term  in  the  sum  converges  in  probability  to 
zero.  We  need  only  show  that 

lG  - (g(Xi,  an))  - Ga(g(X.,  an))| 


also  converges  in  probability  to  zero.  Let 
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1 n * * 

G *(t)  = - Z I[g(X.,  a ) S t],  for  a e RP , 

n, a i=1  1 ~ 

ys)  - sup  |g  (t)  - o (t)|, 

-®<t<°°  n,a+n  s 

and 

«n  ‘ 

P 

It  suffices  then  to  show  that  Q -*■  0.  Now 

n 

P(Q  > e) 
n 

- P(Q_  > e,  n1/2(a  -a)  e B)  + P(Q  > e,  n1/2(a  -a)  t B) 
n _n  _ n _n  _ 

1 /2  * 

S P(Qn  > e,  n (an-a)  e B)  + <5, 
by  3.1.7.  Also, 

P(Q„  > e,  n *(a  -a)  e B) 
n _n  _ 

1 /?  * i ~ 

= P(T  (n  (a  -a))  > e,  n (a  -a)  e B) 
n ~n  _ „n  _ 

£ P(sup  T (s)  > e). 

„ n _ 
seB 

P 

So,  we  can  conclude  our  proof  by  showing  sup  T (s)  -*■  0.  Recalling 

seB  n 


that 
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1 n 

G -1/2  ^ = n 1 I^8(X.,  a+n  1/2s)  £ t], 

n,a+n  s i=l  ~ 


we  have  for  s e B, 


L(t)  * ° -1/2  (t)  S U(t)’ 
n,a+n  s 


where 


1 11 

L(t)  = - l I [sup  g ( X , a+n“1/2s)  < t] 
i=1  seB  1 ~ 


and 


1 n 

U(t)  = — £ I[inf  g (X . , a+n"1/2s)  < t]. 

i=1  seB  1 ~ 


Hence , we  can  bound 


sup|G  (t)  - G (t ) | 

seB  n,a+n"1/2s  “ 


by 


(3.1.8) 


max 


{ | L.  ( t ) - Ga(t)|  , |U(t)  - G (t ) I }. 


Finally,  we  have 
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sup  I L(t)  - G (t) I 

Cl  ' 

-®<t<»  _ 

S sup  | L ( t ) - P(sup  g (X , a+n  1 /2s  < t)| 
-®<tO  seB 


+ sup  | P (sup  g (X , a+n  1/2s)  < t)  - G (t)|. 

-co<t<®  seg  ~ ® 


The  first  term  in  the  above  sum  converges  to  zero  in  probability  by 

lemma  3.1.3,  noting  that  the  distribution  functions  denoted  in  the 

lemma  as  Fn-| » • • • *Fnn  are  identical  and  equal  to  the  distribution 

-1  /2 

function  of  sup  g(X,  a+n  s) , here.  Since,  by  ii, 
seB 


sup  g(X,  a+n  1/2s)  + g(X,  a) 
seB  ~ 


and  G^Ct)  = P(g(X,  a)  £ t)  is  continuous  in  t, 


P (sup  g (X , a+n  1/2s)  < t) 
seB 


G (t), 
a 


for  all  t.  By  Polya’s  theorem,  the  convergence  is  uniform  in  t. 
(See,  e.g.,  Serfling,  1980,  p.  18.)  That  is,  the  second  term  in  the 
sum  above  convergences  to  zero.  Treating  sup  |u(t)  - G (t)|  in  an 

-®<t<®  ® 

analogous  fashion,  we  have,  by  3.1.8,  that 


sup  T (s)  = sup  sup  |G  (t)  - G (t ) | = o (1). 

seB  ~ seB  -°°<t<®  n,a+n  ^s  ® ^ 


37 


In  our  problem  conditions  (i)  and  (ii),  of  lemma  3. 1.4  take  the 

form 


Hq 

(i*)  n 2(E  ^ - eJ)  = 0 (1),  as  n -*•  ® 

o p 

and  (3.1.9) 

Hq 

(ii')  sup|Y'(E  1 + n"1/2s)Y  - Y'  E_1  Y|  = o (1), 
seB  ~ ° ~ ~ o - P 

~ 1 n 

as  n -*•  »,  where  recall  E = - I Y^.',  Eq  = E[ YY '] , and  here  B is  a 

set  of  2x2  symmetric  matrices  whose  elements  are  uniformly  bounded. 

We  establish  condition  (i')  by  showing  that  each  element  of 

n1/2(I  1 - 1 ) is  asymptotically  normally  distributed.  The  following 

lemma,  found  in  Serfling  (1980,  p.  122),  will  be  used. 

Lemma  3.1.10.  Suppose  that  X = (X  ,...,X  )', 

_n  nl  nk 

✓ 1 /2 

p = (p^ , . . . ,p^) ',  and  n (Xn~p)  is  asymptotically  K-variate  normal 

with  mean  0 and  covariance  matrix  E.  Let  g(X)  = (g,(X),...,g  (X))  be 
- 1 _ m _ 

a vector  valued  function  for  which  each  component  g^X),  i = 1,...,m, 

is  real-valued  and  has  a nonzero  differential  at  X = p.  Then, 

1 /2 

n (8(X  ) - g(p))  is  asymptotically  m-variate  normal  with  mean  vector 
0 and  covariance  matrix  DED',  where  the  (i,  j)  element  of  D is 

3gt(X) 

Di,j  = axj  lx  = p- 

Remark  3.1.11.  A sufficient  condition  for  g.  to  have  a nonzero 

1 9gi 

differential  at  p is  that  the  first  partial  derivatives 

~ dX  . 

J 
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1 £ j £ K,  are  continuous  and  not  all  zero  at  y.  (See  Apostol,  1957, 
pp.  110  and  118.) 

We  now  prove  that  conditions  (i')  and  (ii')  hold. 

Lemma  3.1.12.  Let  Y , . . . ,Y  be  i.i.d.  as  Y,  where  Y = (Y  . Y )' 

~ 1 _n  _ i 2 

has  a bivariate  continuous  distribution  such  that  the  covariance 

2 2 

matrix  of  (Y1  , Y^,  Y?)  is  positive  definite.  Then, 

n (E  E "* ) = 0 (1).  as  n ■*  <*> 
o p 

and 


sup|Y'(Zo1+n  1/2s)Y  - Y'  e”1  Y|  = op( 1 ) , 

as  n ->■  »,  where  E , E , and  6 are  as  in  3.1.9. 

Proof.  Write 


\ 

/ ^ 2 

n 

/r ii 

0 

1 

V2’ 

= 

n 

11  2 i 

ti,Y«  / 

and 


°1  2^ 

/e[^] 

e[y,iy2i] 

r 

°22/ 

lE[11 1*2:1 1 

ECY2i]  J 

Then  we  can  write 


22 


r r -r 
11  22 


12 


r r -r 
11  22 


2 

12 


2 

1 2 


-r 


1 2 


r r -r 
1 1 22 


1 1 


r r -r 
1 1 22 


2 

1 2 


2 

1 2 


\ 

/ 
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and 


.-1 


22 


°1 1 °22~° 


1 2 


12 


°1 1 °22_al 2 


12 


\ 


°1 1 °22_°1 2 


'11 


°1 1 a22_°12 


By  the  multivariate  version  of  the  Lindeberg-Levy  Central  Limit 
Theorem  (see,  e.g.,  Serfling,  1980,  p.  28), 


1 /2 

n (r-o)  -*•  N (0,  W), 

H 5 ~ 
o 


(3.1.13) 


where 


_ = (ri1’  r22’  P12)  ’ l (an’  °22’  °12)"’ 


" ■ E[((I?i'  Y2f  4-  yhy2i)-;'». 


and  N^( 0,  W)  refers  to  the  three  dimensional  normal  distribution  with 

mean  vector  0 and  covariance  matrix  W.  We  now  show  that 

1 /2  * * « * 
n (pi  j~°i j ) is  asymptotically  normal  where  r^  and  refer  to  the 

(i,j)  elements  of  E 1 and  E 1 , respectively,  1 £ i < j £ 2.  Taking 

° r2? 

i = j = 1,  define  g(r)  = . Then 

r p -r 
11  22  12 


-r 


22 


9g  | _ 

3rlr=o  , 2 .2  !r  = 

11  _ _ ( r r -r  ) 

v 11  22  12; 


-o 


22 


(°Ua22-°?2>2 


* 0, 


12 


3g  | = 

9r  'r  = a = . 2 ,2  ’ 

22  . . <°,,o22-o12> 


0 
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and 


2°12°22 


11a22"°12) 


2 ‘ 


By  lemma  3.1.10,  remark  3.1.11,  and  3.1.13,  n1/2(r*  -a*  ) is 
asymptotically  normal  and  hence  bounded  in  probability.  In  the  same 
manner,  it  is  easily  shown  that  n1 /2(r*2-o*2)  and  n1 /2(r22~°22)  are 
also  both  bounded  in  probability. 

In  verifying  the  second  result  of  the  lemma,  we  let  M > 0 be  the 
bound  on  each  element  of  the  matrices  in  6.  Then  for  s e B, 


|Y'(I_1+n"1/2s)Y  - Y 'E~ 1 Y 
o _ _ o _ 1 


Irn-’^SYI  . n-,/2|Y^3„^S22.2Y1Y2Sl2| 


-1  /2W  ,,,2  „2 


where 


s = 


2I ) , 

1 1 

S12 

12 

3 22 

That  is 
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sup|Y'(E  ^+n  1/2s)Y  - Y'E_1Y| 
seB  ~ ~ o ~ 


1 n“1/2M(Y2+Y^2|Y1Y2|)  - 0, 


as  n + ®.  _____ 

We  are  now  prepared  to  show  that  n~1/2(C  . S ) and 

MB’  MB' 

-1/2  - 

n ^Cmb ’ SMB^  have  the  same  asymptotic  distribution  under  HQ.  As 
mentioned,  we  are  restricting  attention  to  the  class  of  elliptically 
symmetric  distributions.  The  density  function  of  an  elliptically 
symmetric  bivariate  random  vector,  Y,  is  of  the  form 


fe(y)  = K|zr1/2g((y-0)'E'l(y-0)),  y e R^, 


, -1 


'0  1 


(3.1.14) 


where  g is  a real  valued  function,  K is  a constant,  0 = E[Y],  if  it 
exists,  and  E is  proportional  to  the  covariance  matrix  of  Y,  if  it 
exists.  We  will  assume  for  simplicity  that  Z is  actually  equal  to  the 
covariance  matrix  of  Y.  For  this  class  of  densities  we  have  the 
following  result. 

Lemma  3.1.15.  Let  Y^ Y^  be  independently  and  identically 

distributed  as  Y,  where  Y is  an  elliptically  symmetric  bivariate 
vector  with  mean  0 and  positive  definite  covariance  matrix  Z.  Assume 
also  that  the  covariance  matrix  of  (Y2,  Y^,  Y2)  is  positive 
definite.  Then,  under  H : 0=0, 


-1/2, „ o 

n CMB~CMB: 


P 

) - (0,  0), 


as  n ♦ ®. 


2 


Proof.  Consider  first 


eh  [("'1/2(cMb-eHb»2j 
0 


= E [n-1{  Z a ( cos ( it H ( 4>  ) ) F (D  ) - cos(ttH(  <fr.  ) ) F(D*) ) }2] 
O i = 1 1 n i n l li 


•1  " 2 


= E [n  ' Z a *1  ( cos (irH(<J>))F(D.  ) - cos(ttH(  <J> . ) ) F(D*) ) 2] 

q i = 1'L  nini  i i 


“1 

+ E [n  Z a a (cos(irH  (♦  ))F  (D.  ) - cos  ( ttH  ( 4> . ) ) F(D . ) ) 
o i*j  J nini  i i 


x (cos(TrHn(«})j))Fn(Dj  ) - cos(itH(4>j))F(D*))]. 

(3.1.16) 

Since  Y is  elliptically  symmetric,  Y and  -Y  have  the  same  likelihood 

iinder  HQ.  If  we  let  YP,...,YP  denote,  respectively,  the  projections 

of  the  vectors  Y1,...,Yn  above  the  horizontal  axis,  then  under  HQ 

ai » • • • »an  are  independent  of  the  unsigned  projected  vectors, 

7P,...,YP.  It  follows  that  under  HQ,  a^,...,a  are  independent  of  the 

projected  angles  cf> $ and  the  standardized  distances  D .....D 

in  I n 

and  D , . . . ,D  , since  these  are  all  functions  of  Y?,...,YP  only. 

1 n _ 1 _n 

Therefore,  a^...^  are  under  HQ  i.i.d.  with  Pta^l)  = P(ai=-1) 

= 1/2,  Eh  [a ± ] = 0,  and  the  expectation  given  by  the  second  term  in 
o 

the  sum  on  the  right  hand  side  of  3.1.16  is  zero.  The  first  term  in 
that  sum  is  equal  to 


Eh  [{cos(TrHn(<t)i))Fn(Di)  - cos  (ttH(  (j>1 ) ) F(D*)  }^]. 


(3.1 .17) 


We  will  show  that,  for  all  1 S i $ n. 


cos(TrHn(4,1))Fn(Di)  - cos  (ttH(  4>1 ) ) F(D*)  = o (1),  (3.1.18) 

as  n 4 ®.  We  then  apply  Lebesgue  Dominated  Convergence  Theorem 
(noting  that  the  expression  is  3.1.17  is  bounded  in  magnitude  by  *0  to 
show  that  3.1.17  and  hence  3.1.16  converge  to  zero,  as  n + 

Clearly,  it  suffices  then  to  show  that 


IW  " H(Vl  = °P(1) 


and 


VV  - F<D1>I  - Vn- 


(3.1.19) 


as  n + Lemma  3.1.3  (or  Glivenko-Cantelli ) establishes  immediately 
that 


sup  | H ( t ) - H (t ) | * o (1),  as  n (3.1.20) 

-”<t<”  p 

taking  Fn1,...,F  of  the  lemma  to  each  be  equal  to  H.  Showing 
. * 

'Fn(Di)  “ F(Di)l  = °p(1),  as  n -*■  ”,  requires  the  use  of  lemma  3.1.4. 
Conditions  (i)  and  (ii)  of  Lemma  3.1.4  are  precisely  those  established 
in  Lemma  3.1.12.  This  establishes  that  the  expression  in  3.1.17 
converges  to  zero,  as  n + ”,  and  hence,  by  3.1.16, 

-1/2  - 2 

EH  [ {n  (Cmb“Cmb) ) ] = o (1),  as  n + ».  By  the  Markov  inequality 
o p 

-1/2 

then,  under  Hq  n = 0p^^‘  same  fashion  it  can  be 

shown  that  n1  ^(S^-S^)  = o (1),  under  Hq.  CD 
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- 1 /O 

The  asymptotic  null  distribution  of  n (C  „ , S ) can  now  be 

MB  MB 

found  by  finding  the  asymptotic  null  distribution  of  n~1/2(C  , S ). 

MB  MB 

Theorem  3.1.21.  If  the  conditions  of  Lemma  3.1.15  hold,  then 

under  Hq:  0 = 0,  MB  = 6n  1 * (c^g+s^g)  is  asymptotically  distributed  as 

a Chi-square  random  variable  with  two  degrees  of  freedom. 

Proof.  We  find  the  asymptotic  null  distribution  of 
-11/2- 

(6n  ) (C  , S ) . Let  us  write 

Mb  Mb 


jf,  ei  and  5mb  ■ V 


where 


(3.1.22) 


C.  - 


aicos(TfH(*1))F(D  ) 


and  Si  = a^sinCirHCcJ).  ))F(D  ) , 


i = 1,...,n.  By  the  multivariate  version  of  the  Lindeberg-Levy 

Central  Limit  Theorem,  under  H , 


( 


CMB_nEH  ^Ci-* 
o 


SMB~nEH  '-Si-' 
o 


{nVarR  [C^}172  ’ {nVarH  [S  ]}1/2 

o o 


(3.1.23) 


is  asymptotically  normal  with  mean  vector  0 and  covariance  matrix 

/ " 

1 correlation^  , S ] 


T = 


. (3.1.24) 


correlation!!^ , S ] 
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Since  EH  [a^]  = 0 and  a is  independent  of  both  $ and  D*, 
o 1 * 

EH  ^Ci^  = Eh  ^Si^  = Usin8  the  form  of  the  density  in  3.1.14,  a 
o o 

simple  change  of  variable  argument  will  show  that  A and  D are 

i i 

independent  under  HQ,  for  i = 1 n.  Hence,  for  all  1 S i S n, 

Var  [C  ] = E [a^cos2(irH(<J).  ))F2(D*)] 
o 1 ^o  1 1 1 

= Eh  [cos2(7rH(^.))]EH  [F2(DJ] 
o o 

and  (3.1.25) 


Var  [S  ] = E [sin2(irH(<{>  ) ) ]E  [F2(D*)]. 
o 1 Ho  1 Ho  1 

ft 

Under  HQ,  both  Ht^)  and  F(D..)  are  uniformly  distributed  on  (0,  1)  for 
all  i.  Therefore, 


2 'r  2 

E^  [cos  (irH(<^))]  = J cos  (iru)du  = 1/2 


E^j  [sin  (irH(<^))]  = J sin  (iru)du  = 1/2 
o 0 


and 


(3.1.26) 


Ejj  [F2(D  )]  = j u2du  = 1/3. 
o 0 


Also 


16 


c°vH  [5  , S ] - [C  S ] 
o o 

* Eh  [a^cosCirHC^) ) sin (rrH( <(> . ) )F2(D*)  ] 

= Er  [cos(ttH(  <J>^ ) )sin(irH(  4>^) ) ]En  [F2(D*)] 
o o 1 

f | 2 

= J cos(iru)sin(iru)du  J u du 

0 0 


= 0*  (3.1.27) 

So,  using  3.1.23  - 3.1.27,  (6n  S^g)  is  asymptotically 

bivariate  normal  under  HQ , with  mean  vector  0 and  covariance  matrix 

I2,  the  identiy  matrix  of  order  2.  Lemma  3.1.15  shows  that  under  HQ, 
-1  1 /2 

(6n  ^ (°mb > SMB)  is  also  asymptotically  bivariate  normal  with  mean 

vector  0 and  covariance  matrix  I 2-  It  follows  that 

ns  - ((en-V-Cng)2  ♦ ((bn"1)1^)2 

is  asymptotically  Chi-square  with  two  degrees  of  freedom  under  HQ.  f I 

3. 2 Asymptotic  Distribution  of  MB  Under  Contiguous  Alternatives 
As  the  first  step  in  obtaining  Pitman  asymptotic  relative 
efficiencies,  we  will  find  the  asymptotic  distribution  of  MB  under  a 
sequence  of  alternatives  approaching  the  null  hypothesis.  In  doing 
so,  we  will  consider  a specific  class  of  elliptically  symmetric 
distributions  which  has  a density  function  of  the  form 


2 

} for  y e R , 
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-[(y-e)'if1(y-0)]v 


c(v)exp{- 


h(v) 


(3.2.1) 


where  v,  c = c(v)  and  h = h(v)  are  positive  constants.  Since  MB  is 
invariant  with  respect  to  nonsingular  linear  transformations,  we  can, 
without  loss  of  generality,  take  I = 1^,  the  identity  matrix  of  order 
2.  Letting 


c 


vf(2/v) 

2ir(r(1/v))2 


and 


,2r(1/uK 
^ f(2/v);  ’ 


v 


(3.2.2) 


where 


r(w)  = J uW  "'e  udu  for  w > 0, 
0 


it  can  be  verified  that  the  expression  in  3.2.1  is  a valid  density 
with  covariance  matrix  equal  to  I For  v = 1 the  density  is  that  of 
a bivariate  normal  random  vector.  Taking  v > 1 and  0 < v < 1 produces 
densities  that  are  lighter  tailed  and  heavier  tailed,  respectively, 
than  that  of  the  bivariate  normal. 

With  these  underlying  distributions  we  will  apply  Le  Cam’s  third 
lemma  to  find  the  asymptotic  distribution  of  MB  under  the  alternatives 

Hn:  0 * en  “ (0’  Kn“1/2)'.  (3.2.3) 

(Due  to  the  affine- invariance  of  MB,  we  shall  see  that  the  form  of  the 

alternatives  considered  results  in  no  loss  of  generality.) 

Lemma  3.2.4  (Le  Cam’s  third  lemma;  see,  e.g.,  Prakasa  Rao , 1987, 

p.  104).  Let  p be  a measure  dominating  the  measures  Q and  P . 

n n n 
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Let 


dP 
n 

Pn  ” dy 


and 


dQ 
n 

dy 

n 


be  the  Radon-Nikodym  derivatives  of  and  Q , respectively,  with 
respect  to  y . Also  let 


\ W 

if 

Pn 

> o, 

L = 

1 1 

if 

q 

= p = 0, 

n 

n 

n 

( n 

if 

Qn 

> p = 0. 

n 

Let  S be  a random  variable  and  suppose  that  under  P , (S  , logL  ) is 
n n n ° n 

asymptotically  bivariate  normal  with  mean  vector  y = (y^ , y^) ' and 
covariance  matrix 


where  y = -1/2  Y00.  Then,  under  Q , S converges  in  distribution  to 

a normal  random  variable  with  mean  y^Y^  and  variance  Y^  . 

We  will  take  P^  to  be  the  joint  probability  function  of 

Yl»...,Y  when  0=0,  and  Q , the  joint  probability  function  when 

9=0^.  Taking  y^  to  be  the  Lebesgue  measure  on  Rn,  it  follows  that 

Pn  and  q^  are  the  joint  density  functions  associated  with  Y,j,...,Y 

under  0=0  and  0=0,  respectively.  S will  be  an  arbitrary  linear 
~ _ _n  n 

combination  of  (6n  and  ( 6n 


Hence,  if  we  can  show 
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that  (S^,  logL^)  is  asymptotically  bivariate  normal  under  H , with 

parameters  satisfying  the  relationship  specified  in  Le  Cam's  third 

lemma,  then  will  be  asymptotically  normal  under  the  sequence  of 

(contiguous)  alternatives.  That  is,  we  will  have  shown  that  every 

linear  combination  of  (6n~V/2CMB  and  ( 6n_1 ) 1 /2SMB,  is  asymptotically 

normal  under  H , implying  that  (6n  V/2(C._,,  SUD)  is  itself 

asymptotically  bivariate  normal  under  Hn . 

So,  let  us  denote  the  joint  density  function  of  n independent 

observations  from  3.2.1,  with  Z = I„,  by  p when  0=0,  and  by  q 

2 n - n 

-1/2  ' 

when  0 = 0 = (0,  Kn  ) . Let  P and  Q denote  the  corresponding 

~n  n n 

distribution  functions.  Let 


L 

n 


y 


"n 


n 


I { [ Y ,'Y  ]V 
i=1 


C(Ii‘?n)'(Ii-?n)]Vj] 


(3.2.5) 


Transforming  to  polar  coordinates  via  Y^  - Z^cosfJ).'  and 

# “1  /2 
^2^  = Z^sin<(>',  Z > 0,  and  letting  A = Kn  , we  may  write 

l0gLn  = h 1 {zfV-(Z2-2Zisin<|.i'A+A2)v}.  (3.2.6) 


We  will  first  find  an  expression  which  has  the  same  asymptotic 

distribution  as  logL  under  P . To  accomplish  this  we  use  the 

n n 

following  lemma. 

Lemma  3.2.7.  Suppose  for  each  n M,  X 1 < i -*•<*>  are 

ni  n 

independent  positive  random  variables.  Then 
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T 

n 


K 

n 

Z 

i-1 


P 

-*•  1 


if  for  all  t > 0 

(i) 

and 

(ii) 

Proof.  See  Chow  and  Teicher  (1978,  p.  328). 

An  asymptotic  null  representation  for  logL^  is  given  by  the 
following  lemma. 

Lemma  3.2.8.  If  v > 1/2,  then 

K n 2\>- 1 k2W 

logLn  = — Z 2vsin(<(> ')Z^V  - ^ + op  ( 1 ) , (3.2.9) 

vnh  i=1  n 


K 

n 

Z P(X  > e)  = o(1) 
. . ni 
i = l 


K 

n 

Z E[X  ,I[X  . < 1 ]]  -*>  1 

i = l ni  ni 


where 

W = Ep  [2v(v-1  )Z2v_2sin2(4> ')  + vZ2v'2] 
n 

2 

= K vttc.  (3.2.10) 

Proof.  Letting 

hi(A)  = (Z2  - 2Zisin(<}> ')  A + A2)V, 


we  have 
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dh.(A) 


dA  " v(Zi  ' 2Z1sin(<(i ') A + A2)V_1(2A  - 2Z.sin<})i') , 


and 


d^CA) 


dA 


— = v(v-1)(Z2-2Zisin(4.')A+A2)v"2(2A-2Zisin(<()i'))2 
+ 2v(Z2-2Zisin($')A+A2)V"1 . 


Expanding  h^A)  about  0 via  Taylor's  theorem  yields 


h (A)  = Z2v-2vsin(<J).')Z2v_^  A+A2R 
i i 1 i ni 


where 


* *2  v-2  # 


Rni  = 2v(v-1)(Zi-2Zisin(<|>i')A1+Ai  ) ( A -Z  sin(*') ) 


2 # *2  v_ ^ 

+ v(Z1-2Zisin(<t>i')Ai  + A1  ) , 


for  some  A^  between  0 and  A,  i = 1,...,n.  Hence 


1 n 2 

logL  = ^ I (Z7V-h. (A)) 
n h i i 


a n 2 1 a2  n 

£ l 2vsin(0,;)zf"1  - E R 

h i=1  11  h i=l  ni 


(3.2.11) 


(3.2.12) 


We  need  to  show  that  for  v > 1 /2 
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*2  n 

T £ R„ . 

h i-1  ni 


n k2w 


(3.2.13) 


Write 


ni 


2 * *2  v-1 

v(Z^-2Z^sin(4>.T)  A^  + A^  ) x [1 


2(1-v)(A*-Z.sin(4.'))2 

2~ 

(Z2-2Z.sin(4>')A*+A*  ) 


•] 


2 * #2  V_1 
v(Z^-2Z^sin(4>  ')A^+A^  ) 


A -2Z  sin(4>')A  +Z2sin2(<(> ') 
x [1  - 2(  1 -v)  (— i ■ 1 0 — )]. 


Z2-2Z.sin(<t>')A*+A* 


Since 


(3.2.111) 


2 2 
0 £ (Z2sin2(4.')-2Z1sin4.'A*+A*  ) S (Z^Z^inU'jA^+A*  ) 


(3.2.15) 


it  follows  that  for  1/2  < v < 1, 


2 * *2  V_1 

0 < Rn.  £ v(zJ-2Z  sin(*')A  +A  ) 


< v(Z2cos2  (<(>')  )V_1 


(3.2.16) 


Here  we  are  using  the  fact  that 


Z2-2Zisin(<^')A*+A*  = Yu  + (Y2i~Ai)2 


> Y 


1 i 


= Z2cos2(<t> ') . 
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Using  3.2.14  and  3.2.15,  if  v > 1 and  n is  large  enough  so  that 
Kn  1 /2  < 1 , then 

Rn.  < v(Z^+2Z1+1)V_1  + 2v(v-1)(Z^+2Z.+1)V_1.  (3.2.17) 

Using  3.2.16  and  3.2.17,  we  can  for  large  n bound  Ft  . by 

m 

v(Z^cos2(4>'))V'1  , if  1/2  < v < 1, 

B.(v)  = (3.2.18) 

v(2v-1)(Z2+2Z1+1)V_1 , if  v > 1. 


We  will  show  that  B^(v)  is  integrable  with  respect  to  P^.  Using  again 
the  transformations  = Z^ost^')  and  Y2[  = ZjSint*'),  the  joint 
density  of  Z^  and  <j>'  under  is,  from  3.2.1, 

2v 
z . 

- -r-  (3.2.19) 

f (z  1 , ' ) = cz^  0 < zl  < ®,  0 £ <*>£  < 2ir , 


where  c and  h are  given  in  3.2.2.  The  marginal  densities  of  Z and 
4> ' under  P^  are  given,  respectively,  by 

2v 

_ ZJ_ 

u 

g1  ( z ^ ) = 2-ircz^e  , 0 < z < « 


and 


82(<M  " 2ir ’ 


(3.2.20) 


0 <<*>'<  2ir. 


We  will  be  using  the  fact  that  and  <J>'  are  independent,  as  is 
readily  seen.  We  find 


Ep  [Zf] 
n 


2v 

z 

2irc  j za  + ^e  ^ dz 


a +2  x 

2irc  f 2v  h . 

= 27  J x e dx 


a+2 

trc  .a+2.  2v  , 

- r(— ) h < =0, 


(3.2.21 ) 


if  a > -2.  Also, 


Ep  [ | cos  ( 4> ' ) | 2v  2]  = j |cosu|2v  2du 
n 0 


H /2  9 0 IT 

p—  { J (cosu)  v du  + J (-cosu ) 2v  2du} 
0 tt/2 


? 11  f 2 ? 9 

£ j (cosu)^v  ^du 
17  0 


i \ yCC! « 

0 2w  ( 1 -w) 


sin2u) 


i) 


(1-w) 


(v-  ^)-1  -1 


w dw  < ® , 


(3.2.22) 


if  v > 1/2.  We  see  then,  using  3.2.18,  3.2.21,  and  3.2.22  that  B7v) 
is  integrable  with  respect  to  P if  v > 1/2. 
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We  will  now  use  Lemma  3.2.7  to  show  that  for  v 


n 

Z 

i-1 


ni 


P 

n 


or  equivalently  that 


n 

Z 

i-1 


ni 


n 

Z 

i-1 


R . 
ni 

Wn 


P 

n 

-*■  1 . 


Taking  Y^.  = we  will  first  show 


n 

Z 

i-1 


£ e) 


o(1), 


for  all  e > 0.  We  have 


I Pp  (Y  i e)  - I P (-Si  2 „) 
i-1  n ni  i-1  n W£ 


n B (v) 

* Z PP  ("1 7T  * n) 

i-1  n e 


B.  (v) 

npp  (~w r - n) 

n 


> 1/2 


- 0, 
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B.(v) 

as  n -*■  ®,  since  — is  integrable.  As  for  condition  (ii)  of  Lemma 
3.2.7,  write 


E E [Y  . I[Y  < 1]] 
. ..  P ni  ni 
i = 1 n 


n R R 

,£,  epH^ 

i = 1 n 


"ep  I[%< 1J] 

n 


W EP  CRni  I[  Wn  < 1 
n 


Since  R £ B (v)  and  E_  [B.(v)]  < ®,  we  can  use  Lebesgue  Dominated 
ni  1 P i 

n 

Convergence  Theorem  to  obtain  the  result  that 


lim  I Ed  [Y  I[Y  ( < 1]] 
. , P ni  ni 

n-*-®  i = i n 


lim  W EP  ^ni1^  nW  < 1 -1-1 
n+®  n 


^ Ep  [2v(v-1)Z2v_2sin2U')+vZ2v~2] 


1. 


where  we  are  using  3.2.11  to  find  the  stochastic  limit  of  R , as 
n -*■  ®.  In  conjunction  with  3.2.12,  this  establishes  the 
representation  of  logLn  given  by  3.2.9.  The  expectation  in  3.2.10  is 
easily  established  using  3.2.21  and 


j sin2(<J) ')d<j> ' = 1/2. 
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gi 


P 

n 


The  asymptotic  P distribution  of  IokL  easily  follows  and  is 
n n 

ven  in  the  next  lemma. 

Lemma  3.2.23.  If  v > 1/2,  logL^  is  asymptotically  normal  under 

P - 1 

with  mean  y and  variance  -2y,  where  y = K Wh  . 

Proof.  Let 


* -1/2  n 2-1 
L = n Z L - K Wh 

n i-1  1 


where 


(3.2.24) 


L 


i 


2vKh 


sin($')Z 


2v-1 

i 


i 


* 

By  the  central  limit  theorem  L^  is  asymptotically  normal  under 
with  mean 


P 

n 


Ep  [L±]  - K^"1 


and  variance 


(3.2.25) 


Varp  [L  ]. 
n 


We  find 


Ep  [L1]  = 2vKh_1Ep  [sin(4> ')Z^V_1  ] = 0,  (3.2.26) 

n n 


- T 

2 IT  l 


sin(<|»i')d<()' 


0, 


since 
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and  <{>'  and  are  independent.  Also, 

Varp  [L  3 = Ep  [ L2]  = (2vKh_1)2E  [sinV']E  [Z2(2v_1)] 
n n n n 1 

= (2vKh_1)2(l)(^  r(2)  h2) 

= 2vK2iro,  (3.2.27) 

using  3.2.21  and 


2tt 


J sln2(*')d*' 


1 /2. 


By  Lemma  3.2.8,  logL^  is  under  P^  also  asymptotically  normal  with  mean 

2-1  2 
-K  Wh  and  variance  2vK  ire.  Furthermore, 

-K2Wh_1  = -K2h_1Ep  [2v(v-1)Z2v_2sin2(<(.')+vZ2v_2] 
n 

= -K2h_1Ep  [Z2v_2]Ep  [2v(v-1)sin2(<K)+v] 
n n 1 

- (-K2h~1)(irchv“1)(2v(v-1)(^)+v) 

„2 

= -K  ircv 

= - jj(Varp  [Li ] ) . (3.2.28) 

n 

This  completes  the  proof  and  implies  that  Q is  contiguous  to  p (see. 

n 

e.g.,  Prakasa  Rao,  1987,  p.  104).  I I 

The  following  theorem  gives  the  asymptotic  distribution  of  MB 

under  Q . 
n 
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Theorem  3.2.29.  Suppose  Y, f...,Y  are  independently  and 

_1  _n 

identically  distributed  from  a density  function  given  by  3.2.1,  with 
v > 1/2,  and  E = I^.  Then,  under  the  alternatives  given  by  3.2.3,  MB 
is  asymptotically  distributed  as  a noncentral  Chi-square  random 
variable  with  two  degrees  of  freedom  and  noncentrality  parameter  given 
by 


smbm 


vK/3  r(^  +i ) (r(-) )1  /2b  (^) 
2v  v v 2 


(r(-)) 

V 


5/2 


where 


(3.2.30) 


00 


r(w)  = J 

o 


w-1  -u 
u e 


du , 


w > 0, 


and 


f v ,2v 

J u (1-u)  du. 

0 


Proof.  Applying  Le  Cam's  third  lemma  (Lemma  3.2.4),  we  need  to 

find  the  asymptotic  distribution  of  (S  , logL  ) under  P , where 

n n n 

S„  ' h<6n'',1/2cMBU2(6"'1),/2sMB 

-11/2  -11/2 

is  an  arbitrary  linear  combination  of  (6n  ) C.m  and  ( 6n  ) Swo. 

MB  MB 

By  3.2.24,  3.2.9,  and  Lemma  3.1.15,  the  asymptotic  distribution  of 

(S  , logL  ) under  P is  the  same  as  that  of  (S  , L ),  where 
n n n n n 
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S*  - M6n-V/2CHBU2(6n-V/2iHB. 


n 1 


* * 

We  will  then  work  with  (S  , L ) . Let  us  write 

n n 


:mb  ‘ i-£,  Ci 


a"d  §HB  ' .E,  V 
1 = 1 


where 


and 


C.  = aicos(irH((>i))F(D*) 


Si  = aisin(nH(«fr  ))F(D*). 


(3.2.31) 


We  have  already  established  that 


and 


Ep  CC  ] = Ep  [S  ] = 0 
n n 


Varp  [C±]  = Varp  [S  ] = ^ . 
n n 

(See  the  proof  of  Theorem  3.1.21.)  Therefore,  using  the  multivariate 
version  of  the  Central  Limit  Theorem,  in  conjunction  with  Lemma 

ft  * 

3.2.23,  we  have  that  (S  , L ) and  hence  (S  , logL  ) are  both 

n’  n n 6 n 

asymptotically  bivariate  normal  under  with  mean  (0,  y)  and 
covariance  matrix 


12 


-2y 


r 


(3.2.32) 
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where 


Y12  = C0VP  C(V6  VV6  V«  Li] 

n 


= Ep  [(£  /5  C +1  /&  S )L  ] 
n 


(3.2.33) 


To  find  Y^,  we  need  Ep  [C^L^]  and  Ep  [S^L^].  Refer  to  the 

n * n 1 

definitions  of  , and  D^,  given  in  Chapter  2,  and  recall  that 

E = I_.  Then 
o 2 


* -1 
d4  = y:e  y. 

i _i  o _i 


2 2 2 

Y1i+Y2i  = V 


i 1|  •••  yii| 


and  (3.2.34) 

F(t2)  = Pp  (D*  < t2)  = Pp  (Z2  < t2)  = Pp  (Z  S t), 
n n n 

for  t >_  0.  Noting  that  the  <J>'  defined  in  conjunction  with  MB  is 
identical  to  that  used  in  the  preceding  lemma,  we  use  2.1.1  and  3.2.20 
to  find,  for  0 £ t ^ it, 


H(t ) ■=  Pp  (<t>1  < t)  = Pp  (0  < <}>'  < t)  + Pp  (ir  < < ir+t ) 

n n n 

t TT  + t IT 

2tt  2tt  2tt 

= - . (3.2.35) 


Letting  G(t)  *=  Pp  (Zi  £ t), 
n 
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Ep  [CiLi]  = 2vKh_1Ep  [a1cos(irH(*1))F(Dpsin(*')ZjV"1] 
n n 

= 2vKh"1Ep  [a1oos(*1)sin(*p]Ep  [Z*V-1G(Z  )], 
n n 

* 

using  3.2.3M  and  the  independence  of  {D.  , Z. } and  {a^  4^,  4 

P . Since 
n 


Ep  [a^cos  (<|Osin(<t>')] 
n 


Ep  [cos  (4) ' )si n ( 4>  ')I[0  £ <p ' S it]] 
n 


+ Ep  [-cos  (4 it  )sin(4')I[ir  < 4> ' < 2ir]] 


and 


2n 


2~  j co  s(4>')sin(4>')d4)'  = ^ “ J _cos  ( <}>  '-ir)sin(4 ' )d<j> ' = 0 


we  have 


Ep  [CjL  ] = 0 
n 


Finally , 


Ep  CS.Li]  = 2vKh_1Ep  [G(Z.)Z^v_1]Ep  [a.sin(irH(  4 ) )sin(4 ')  ]. 
n n n 


(3.2.36) 
' } , under 


(3.2.37) 


(3.2.38) 


(3.2.39) 


We  find 
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Ep  [aisin(-nH(<t>i) ) sin ( <J> ')  ] 
n 

2 

= Ep  [sin  (<}>')I[0  £ <p  ' < ir]]+Ep  [-sin(<(> '-ir)sin(<}) ')I[tt  < $ < 2ir]] 
n n 1 

= 2 n f sin2(4> ' ) d ( 4> ')  - — 1 J sin ( 4>  '--it ) sin ( 4> ' )d ( 4> ') 

0 IT 

1 ? 2 

= i J sin^ (4, ')d (<(.') 

- h (3.2.40) 


and  using  3.2.20  again, 


Ep  [G(Zi)Z2v_1] 
n 


2v 


= 2uc  f G(z)z  ve 


2v  h 


dz 


2v  2v 

y z 

2irc  f (2irc  J ye  ^ dy)r2v  h 
0 0 


iz  e 


dz , 


1 1 


1 1 


Letting  y = xVsinV(u)  and  z = xVcosV(u),  we  find  that 


1-1  1 
dz  _ xv  cosv(u) 
dx  v 


1 1-1 

dz  _ -x  cos  (u)sin(u) 
du  v 


1-1  1 

dy  _ x sin  (u) 
dx  v 


1 1-1 

and  <*1  = xvsinv  (u)cos(u) 
du 


v 


v 


The  Jacobian  of  the  transformation  is 
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1 1 , 1 , 1 , 1 . 

v - +1  - -1  1 — +1 

- — (cosv  (u)sinv  (u)  + cosv  (u)sinv  (u)) 

v 


- -1  1 , 1 , 

v - -1  - -1 

- — ^-(cosv  (u)sinv  (u))(cos  (u)+sin2(u)) 


- -1  1 . 1 . 

v - -1  - -1 

- — (cosv  (u)sinv  (u)). 

v 


It  follows  that 


Ep  [G(Z1)Z^V_1] 
n 


4it  c 


2 2 » u/4  - -1  - -1 


1-1 


| J xv  cosV  (u)sinv  (u) (xsin(u) ) V(x2cos2(u) )e  hdudx 


v 0 0 


h 2 2 • - +1  - tt/4  1-1  1+1 

-( j xv  e dx)(  j sinv  (u)cosv  (u)du). 


v2  0 


The  first  integral  in  the  product  above  is 


- 1 +1  _ x_ 

J xv  e h dx 


* _3  _ I 

1 f 2v  h 
P J y e dy 


h2v 


+1 
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The  second  integral  is 


tt/4 


? 1 £ 1_ 

- -1  - + 1 1/2  , 1/2, v . 1 /2v  v +1 

/ sinv  (u)cosV  (u )du  = j — dw 

o o 2w  ^(i-w y,<L 


1/2  — -1 

= P { wv  ( 1-w)  v dw, 


using  the  transformation  w = sin  (u).  Hence, 


Ep  [G(Zi)Z^v_1] 
n 


2 2 — + 1 
iL_2_rl  +i'iv12v  wl  r 

2 ^2  r^2v  1 h 2 Bv  2 ’ 

v 


where 


(3.2.41 ) 


1/2  1-1  J- 


Bv(l)  = J wV  (1-w)2v  dw, 


It  follows,  using  3.2.39  - 3.2.41,  that 


Covp  [Sif  Li] 
n 


3 2 1/21 

llf  r(^4  ) (r(-) )n (1) 

vK  2v  v v 2 


/2 


( r ( — ) ) 

V 


5/2 


and 


(3.2.42) 


712  = CoVp  [1/6  C^/6  Sif  L±] 
n 

= l^/T>  CoVp  [S  , L/. 
n 
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We  have  found  by  Le  Cam's  lemma  that  Sn  = ^ ( 6n_1  )1 /2CMB+2,2( 6n_1 ) 1 /2SMg 

is  asymptotically  normal  with  mean  Y,  and  variance  J if+i,2  under  0 

12  12  n 

Since  Cov  [C  , L ] =0,  it  follows  that  ( 6n_1  ) 1 /2 (C , SUD)  is 
r Yi  1 1 mb  mb 

under  asymptotically  bivariate  normal  with  mean  vector 

(0,  /6  CoVp  [Sif  L ])  and  covariance  matrix  I . We  conclude  that 
n 2 

under  Qn 


MB  - (6n"')2(C^B»S^B) 


is  asymptotically  Chi-square  with  two  degrees  of  freedom  and 
noncentrality  parameter 


(v) 


(/6  CoVp  [§,  L^)2 

n 

vk/3  r(^  *,)(r(f))2  6v(l)'2 


5/2 


(r(-)) 

v 


□ 


(3.2.43) 


3.3  Asymptotic  Distribution  of  Blumen's  Sign  Test 
Blumen  (1958)  found  the  asymptotic  null  distribution  of  his 
statistic,  B,  to  be  Chi-square  with  two  degrees  of  freedom.  We  will 
derive  this  result  in  a different  fashion  so  that  we  can  again  apply 
Le  Cam's  third  lemma  to  find  the  asymptotic  distribution  of  B under 
the  contiguous  alternatives  considered  in  section  3.2.  Using  the 
notation  in  the  previous  section,  recall  that 


B = 2n  ( Cg+ Sg ) , 


where 


(3.3.D 
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n 


n 

l 

i=1 


a,cos(irH  (*.)) 
i n 1 


and 


SR  = ~ E a sin(irH  ($.)). 
ts  n l n i 

Let  us  define 


C 


B 


n 


and 


l 

n 


n 

I 

i = 1 


where  (3-3.2) 

CBi  = a^cos  (irH(  ) ) 
and 

SBi  - a^sin(irH(  (Ji^ ) ) . 

In  the  following  lemma  we  show  that  /n(C  , S ) and  /n(CD,  S_) 

B B B B 

have  the  same  asymptotic  null  distribution. 

Lemma  3.3.3*  Suppose  Y^,...,Y  are  independently  and  identically 

distributed  from  a continuous  bivariate  distribution  with  median 

0 = (0  , 0 )'  satisfying  1.1.4.  Under  H : 0=0, 

_ i d o 
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P 

^(VV  SB-SB)  - (0,  0), 

as  n + ». 

Proof.  As  in  the  proof  of  Lemma  3.1.15,  write 

eh  [<^<cb-cb»2] 

0 


= nE  [{1  z a [cos(irH  U ))-cos(irH(<t>.  ))]}2] 

H I n 1 1 


(3.3.-4) 


Recalling  that  a1t...,an  are  i.i.d.  under  HQ  with  [a.]  = 0,  and 

o 

that  a-, an  are  independent  of  4^,...,$  it  follows  that  the 

second  term  in  the  sum  above  is  equal  to  zero.  Hence, 


eh  [{/n(CB-CB)}2]  = Eh  [ {cos(trHn(4)i))-cos(TrH(4>  ))  }2].  (3-3.5) 

o o 

H 

o 

By  Glivenko-Cantelli,  as  in  3-1-20,  H (<fc  )-H(<|>  ) = o (1),  hence 

n 1 i p 


H 

o 

cos(irH  ( <f>  ) )-cos(irH(  A ) ) = o (1). 
n 1 ip 


Using  Lebesgue  Dominated  Convergence  Theorem  it  follows  that  the 
expression  in  3-3-5  converges  to  zero,  as  n -*•  ».  This  shows  that 


69 


H 

/"<V5b> 

and,  in  an  analogous  fashion,  it  is  seen  that 


P 


Lemma  3.3-3  readily  shows  that  the  asymptotic  distribution  of  B 
is  Chi-square  with  two  degrees  of  freedom.  It  is  easily  verified  that 
under  HQ  Cov[CB.,SB.]  = 0,  E[CB±]  - E[SBi ] = 0,  and  Var[CBi]  = 

Var C SBi ] = 2 • Hence  by  the  Multivariate  Central  Limit  Theorem  the 
asymptotic  null  distribution  of  /2n(CB>SB)  is  bivariate  normal  with 
mean  0 and  covariance  matrix  I2.  The  asymptotic  null  distribution  of 
B follows. 

In  the  following  theorem  we  use  Lemma  3.3.3  to  find  the 
asymptotic  distribution  of  B under  the  contiguous  alternatives  of 
section  3.2. 

Theorem  3.3.6.  Let  Y^,...,Y  be  independently  and  identically 
distributed  from  a density  function  given  by  3.2.1,  with  E = Var[Y.] 

= I2-  If  v > 1/2,  then  under  Hn:  0 = = (o,  Kn~1/2)',  B is 

asymptotically  Chi-square  with  two  degrees  of  freedom  and  noncentrally 
parameter 


where 


62(v) 


vK  r(-5~  + 1 ) ( r ( — ) ) 
2v v 

( r (— ) ) 3/2 

V 


(3.3.7) 


f(w) 


w-1  -u 
u e 


du , 


w > 0. 
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Proof.  We  will  follow  the  procedure  and  use  the  same  notation 
given  in  the  proof  of  Theorem  3.2.29.  Applying  Le  Cam's  third  lemma, 
we  need  to  show  that 


U /2nC  +£  /2nS_,  logL  ) 
i ts  d.  d n 

is  asymptotically  bivariate  normal  under  P^.  Using  Lemmas  3.2.8  and 
3-3.3,  it  suffices  to  show  that 

u/S  cB.y^  sB,  l*> 

£ 

is  asymptotically  bivariate  normal  under  P , where  L is  given  by 

3.2.24.  This  is  immediate  from  the  Multivariate  Central  Limit 

Theorem.  Having  already  verified,  in  Lemma  3.2.23,  that  the 

parameters  of  the  asymptotic  P distribution  of  logL  satisfy  the 

relationship  required  in  Le  Cam’s  lemma,  we  conclude  that 

i /2n  C + 8,  /2n  S and  therefore,  8./2n  C +8„/2n  S„,  are  each 
• ts  <L  ti  l B 2 B 

asymptotically  normal  under  Q with  mean 

n 

Ep  H^/2  CBi+«-2/2  SB.]  + CoVp  [(8-^2  C^  + l^/2  SgJ.L.] 
n n 

and  variance  (3.3.8) 

VaPp  Si*'/5  V’ 

n 

where  Li  is  given  by  3.2.24.  It  is  easily  verified  that 
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EP  [ V*  W5  V ' 0 


and 


(3.3.9) 


VarP  ll^  W5  V ‘ 

n 


We  find 


Covp  [5Bi-  Ll] 
n 


EP  CCBiLi] 


2vKh  1 Ep  [a^osUHt  4^ ) )sin(<f> ')  ]Ep  [Z2v  1 
n n 


- 0, 


(3.3.10) 


using  3.2.35  and  3.2.37.  In  addition, 


Covp  [SB1,  L;]  - Ep  [S  L ] 
n n 

= 2vKh_1Ep  [aisin(irH(4.i))sin(4»')]Ep  [Z2v_1] 
n n 

» 2vKh_1(1/2)(™  T(^  +1)  h^  1}),  (3.3.11) 


using  3.2.21  and  3.2.40.  Thus,  using  3-3.8  - 3.3.10,  l /2n  C+ZV 2n  SD 

B 2 B 

is  asymptotically  normal  under  Q with  mean  iV2Cov„  [S  L 1 and 

n 2 P Bi,  i 

2 2 n 

variance  Z +Z  . Therefore,  under  Q , /2n(CD,  S)  is  asymptotically 

n B B 
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bivariate  normal  with  mean  (0,  /2  Cov[SBi,  L.])  and  covariance  matrix, 

I . The  asymptotic  distribution  of  B = 2n(C2+S2)  under  Q is  then 
^ B B n 

noncentral  Chi-square  with  two  degrees  of  freedom  and  noncentrally 
parameter , 


6b(v)  = 2(Cov[Sb.,  L^])2 


— +1 

2vKh_1(l)(—  r(-l +1)  h2v  ) 
2 v 2v 


vK  r(^  +i)  (r(-))1/2  '2 

V 

( r (— ) ) 3/2 

V 


(3.3.12) 


CHAPTER  4 

COMPARISONS  OF  ONE-SAMPLE  TEST  STATISTICS 

1 Asymptotic  Relative  Efficiencies 

In  this  section  we  use  Pitman  asymptotic  relative  efficiencies  to 

2 

make  comparisons  among  Hotelling's  T , the  sign  statistics  of  Blumen 

and  Hodges,  and  the  proposed  modification  of  Blumen' s statistic.  (See 

Chapter  1 for  definitions  of  Hotelling's  T2  and  Hodges'  statistic, 

H.)  Because  these  statistics  are  all  affine-invariant  we  may,  without 

loss  of  generality,  make  the  simplifying  assumption  that  the 

population  covariance  matrix  is  the  identity.  Further,  alternatives 

of  the  form  0 = (K^  , K^) ' can  be  transformed  to  equivalent 

alternatives  of  the  form  0 = (0,  (K2+K2)172)  by  a careful  choice  of  an 

orthogonal  transformation  (rotation)  of  the  observations.  Hence  we 

will,  without  loss  of  generality,  consider  alternatives  of  the  form 

0 = (0,  A) ' and  we  will  apply  the  asymptotic  results  under  the 

contiguous  alternatives  discussed  in  Chapter  3. 

The  Pitman  approach  to  asymptotic  relative  efficiency  compares 

two  test  sequences  and  T^#  as  the  alternative  hypothesis, 

n n 

Hn:  0 = 0n»  aPPr°aches  the  null,  which  we  are  taking  to  be 

£ 

H ; 0=0.  The  subscripts  K and  K are  the  sample  sizes  for  tests 

o n n 

S and  T , respectively.  Let  Power  S (0  ) and  Power  T*(0  ) denote 
^ K _n  K*  _n 

n n n n 

the  probabilities  that  tests  S and  T , respectively,  reject 

K K* 

n n 
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H when  0=0.  Assume  that  K and  K are  such  that  the  limiting 
o ~ _n  n n 6 

significance  levels  are  the  same  for  each  test  sequence,  as  are  the 
limiting  powers: 

lim  Power  S (0  ) = lim  Power  T„„(0  ) = B, 
n+°°  n n+°°  n 

where  0 < a < B < 1,  and  a denotes  the  limit  of  the  significance 

levels  of  and  T , as  n + <*>.  The  Pitman  asymptotic  relative 
n n 

efficiency  of  S relative  to  T „ is 
K K* 

n n 

* 

K 

ARE(S,  T)  = lim  ^ , (4.1.1) 

n-*’®  n 


provided  the  limit  exists  and  is  the  same  for  all  {0  },  {K  } and  {K*} 

_n  n n 

satisfying  the  conditions  above.  (See,  e.g.,  Randles  and  Wolfe,  1979, 
p.  144.) 

Hannan  (1956)  shows  that  if,  under  Hn,  the  test  sequences  and 

n 

tk*  are  asymptotically  noncentral  Chi-square  with  v degrees  of  freedom 

" 2 2 
and  noncentrality  parameters,  6g  and  6 , respectively,  then 


ARE( S,  T)  = 


(4.1.2) 


It  is  well-known  that  the  asymptotic  distribution  of  Hotelling's  T2, 

-1  /2 

under  H^:  0 = 0^  = (0,  Kn  ) , is  noncentral  Chi-square  with  two 

degrees  of  freedom  and  noncentrality  parameter, 


6 


2 

T 


(4.1.3) 


(See,  e.g.,  Puri  and  Sen , 1 971  , p . 1 73- ) We  can  use,  then, 
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Theorems  3.2.29  and  3.3.6  to  find  the  asymptotic  relative  efficiencies 

p 

of  B and  MB  relative  to  T . These  results  are  stated  in  the  following 
theorem. 

Theorem  ^.1.^.  Suppose  Y^,...,Y  are  independently  and 
identically  distributed  from  a density  given  by  3.2.1,  with  I = I^. 

For  testing 

H : 0 = 0 against  H : 0 = 0, 

o a _ _n 

if  v > 1/2,  then 


ARE( B,  T ) = 


62bM 


.i)(r(f))1/2 

(r(-))3/2 

V 


and 


ARE(MB,  T2)  = 


(v) 

2 

T 


^r(5a+i)(r(i))1'2#v(i)| 

(r(l))5/2  • / 


2 


where 


T(w) 


r w-1  -u. 

= J u e du, 

0 


w > 0, 


6v(2} 


i 1 

f uV 


( 1-u)2vdu. 


and 
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Table  4.1  provides  the  Pitman  asymptotic  relative  efficiencies  of 
B and  MB  relative  to  T2,  for  selected  values  of  v. 


Table  4.1 

Pitman  Asymptotic  Relative  Efficiencies 


V 

ARE(B,  T2) 

ARE (MB , T2) 

10.0 

.50530 

1 .33231 

5.0 

.51872 

1.22101 

1.0 

.78540 

.98464 

.8 

.91839 

.99574 

.7 

1.02988 

1 .01488 

• 5 + 

1.50000 

1.12500 

We  see  that  although  Hotelling’s  T2  performs  best  for  normal 

alternatives  (v  = 1.0),  MB  appears  to  be  quite  competitive.  Except 

for  distributions  very  close  to  normal,  MB  is  asymptotically  more 

powerful  than  Hotelling's  T2.  While  Blumen's  statistic,  B,  performs 

poorly  for  lighter-tailed  distributions,  it  performs  better  than 
2 

Hotelling’s  T and  MB  for  heavy-tailed  distributions.  The  Pitman 

asymptotic  relative  efficiency  has  not  been  found  for  Hodges' 

statistic.  However,  Bahadur  (I960)  notes  that  under  certain 

conditions  the  limiting  Bahadur  efficiency  coincides  with  that  of 

Pitman.  Killeen  and  Hettmansperger  (1972)  found  the  limiting  Bahadur 

efficiency  of  Hodges'  statistic,  H,  relative  to  Hotelling's  T2  to  be 
2 

- under  normal  alternatives.  It  appears  then  that  under  normal 
alternatives  B,  MB,  and  t2  are  each  superior  to  H asymptotically. 


77 


We  have  also  displayed  in  Table  4.2  conjectured  Pitman  asymptotic 

relative  efficiencies  when  0 < v £ using  the  form  of  the 

efficiencies  given  in  Theorem  4.1.4.  For  heavy-tailed  distributions, 

further  evidence  of  the  superior  performance  of  B and  the  inferior 
2 

performance  of  T is  indicated. 


Table  4.2 

Conjectured  Pitman  Asymptotic 
Relative  Efficiencies 


V 

ARE( B,  T2) 

ARE (MB,  T2) 

.4 

2.09838 

1.27541 

.3 

3.69648 

1.64385 

.2 

11.59689 

3.00306 

.1 

366.57838 

25.03483 

4.2  Monte  Carlo  Results 

Additional  comparisons  of  several  one-sample  test  procedures  will 
be  made  in  this  section  using  Monte  Carlo  methods.1  The  affine- 
invariant  statistics  of  the  last  section,  namely  B,  H,  T2,  and  MB  are 
examined  along  with  two  statistics  which  are  not  affine-invariant. 
Bennett's  sign  test  statistic  B1  (see  1.1.5)  is  included  along  with 
his  signed-rank  test,  which  was  also  introduced  in  Chapter  1.  The 

1 All  calculations  and  random  variable  generations  were  performed 

with  an  IBM-4341  running  on  a vm/cms  operating  system.  Fortran  77 
and  several  IMSL  library  subroutines  were  used. 
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precise  form  of  the  signed-rank  test  is 

B2  = [(^(U*.  V+)]'(  nS)"1[(JT)(U+,  V+)  ] , 

, + + 

where  U and  V are  the  signed  rank  statistics  for  the  marginal 

observations  as  described  in  1.1.6,  and  S is  a consistent  estimator  of 

the  asymptotic  covariance  matrix  of  (u  + , V+).  The  estimator, 

/n (n+1 ) 

S,  we  use  has  diagonal  elements  equal  to  g2^  ^ , and  off  diagonal 
elements  equal  to 

1 n 

x I Rank(|z  |)Rank(|z  | )I[Z  > 0]I[Z_.  > 0]. 

n (n+1 ; i=1  11  2l  11  2i 

(See,  e.g.,  Hettmansperger , 1984,  p.  284.)  These  six  test  statistics 
were  compared  for  samples  taken  from  several  distributions,  most  of 
which  were  of  the  form  given  in  3.2.1.  Additional  comparisons  were 
made  using  samples  from  Pearson  Type  II  and  Type  VII  (see  Chapter  6 of 
Johnson,  1987),  and  mixtures  of  bivariate  normal  distributions.  (A 
bivariate  normal  mixture  is  obtained  by  selecting  at  random  one  of  two 
bivariate  normal  distributions.  Each  of  the  sample  elements  is 
selected  with  probability  p from  the  first  distribution  and  with 
probability  1-p  from  the  second  distribution.  See  Johnson,  1987, 
Chapter  4.)  In  each  Monte  Carlo  simulation,  the  proportion  of  times 
each  test  statistic  exceeded  the  upper  a-percentile  of  its  asymptotic 
null  distribution  is  reported.  The  rejection  proportions  for  T^  using 
critical  values  from  its  distribution  under  normality  (namely,  a 
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multiple  of  n_2)  are  given  in  parentheses.  The  alternatives 
considered  were  of  the  form  0 =(0,0)',  which  results  in  no  loss 

-3.  3 3 

in  generality  for  the  affine-invariant  statistics.  For  most  of  the 

simulations,  values  for  0 were  selected  so  that  0'Z_10  equaled 

a ~a  _a 

2 2 2 

•2  , .4  , or  .6  , where  I = Var[Y]  was  taken  to  be  the  identity. 

Heavy-tailed  distributions  were  treated  somewhat  differently.  For 

settings  in  which  Var[Y]  does  not  exist,  the  random  variables  were 

scaled  to  produce  readily  comparable  results.  For  the  heaviest-tailed 

distribution  examined  from  the  class  given  in  3.2.1  (v  = .10),  the 

amount  of  translation  was  altered.  Significance  levels  were  chosen  to 

agree  with  available  asymptotic  percentiles  for  Hodges’  null 

distribution.  (See  Joffe  and  Klotz,  1962.) 

In  Table  4.3,  the  results  from  the  first  Monte  Carlo  study  are 

presented.  As  indicated,  1000  samples  of  size  25  were  obtained  from 

five  members  of  the  class  of  distributions  given  by  3.2.1,  with 

Z = I2.  Heavy-tailed  distributions,  v = .10  and  .50,  the  bivariate 

normal  distribution,  v = 1.00,  and  light-tailed  distributions, 

v = 50.00  and  120.00,  were  included.  The  bivariate  medians  of  the 

underlying  distributions  and  the  nominal  significance  level,  a,  are 

indicated.  Although  it  is  difficult  to  compare  empirical  powers  for 

statistics  with  varying  rejection  proportions  under  the  null 

hypothesis,  it  appears  that  the  Monte  Carlo  results  tend  to  agree  with 

the  asymptotic  results  of  the  previous  section.  Hotelling’s  T2  has 

the  greatest  power  for  normal  alternatives.  For  distributions  close 
2 

to  normal,  T and  MB  perform  similarly,  both  performing  better  than 
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Table  4.3 

One-Sample  Monte  Carlo  Results  with  n = 25  and  a = .05318 


Statistics 

Bivariate 

medians 

B 

H 

T2  MB 

B1 

B2 

( .000, .000) 

.052 

.045 

.047 

v = .10 
( .01 8) 

.060 

.047 

.051 

( .01 4, .01 4) 

.144 

.136 

.061 

( .026) 

.087 

.128 

.086 

(.028, .028) 

.304 

.270 

.074 

(.041) 

.138 

.278 

.148 

( .042, .042) 

.478 

.437 

.096 

( .064) 

.21 1 

.419 

.246 

(.000, .000) 

.052 

.045 

.080 

v = .50 
(.049) 

.060 

.047 

.048 

( .141  , .141  ) 

.154 

.140 

.192 

( .142) 

.151 

.118 

.160 

(.283, .283) 

.485 

.449 

.515 

(.459) 

.443 

.398 

.453 

( .424, .424) 

.809 

.760 

.806 

(.745) 

.726 

.726 

.776 

(.000, .000) 

.052 

.045 

.084 

v = 1.00 
(.054) 

.058 

.047 

.048 

( .141  , .141  ) 

.107 

.105 

.181 

(.134) 

.135 

.091 

.119 

(.283,-283) 

.313 

.289 

.468 

(.381) 

.378 

.268 

.361 

( .424, .424) 

.61  2 

.574 

.798 

( .724) 

.708 

.528 

.704 

(.000, .000) 

.052 

.045 

.080 

v = 50.00 
( .062) 

.060 

.047 

.051 

( .141  , .141  ) 

.082 

.076 

.171 

( .127) 

.145 

.069 

.108 

(.283, .283) 

.202 

.198 

.448 

(.362) 

.441 

.157 

.306 

( .424, .424) 

.442 

.401 

.794 

(.721) 

.796 

.366 

.639 

(.000, .000) 

.052 

.045 

.088 

v = 120.00 
(.054) 

.061 

.047 

.048 

(.141,.  141) 

.076 

.072 

.170 

(.122) 

.166 

.057 

.104 

(.283, .283) 

.140 

.132 

.431 

(.351) 

.514 

.110 

.302 

( .424, .424) 

.298 

.265 

.785 

(.703) 

.864 

.234 

.624 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are  the 
rejection  proportions  for  Hotelling's  T2  using  the  F2  n_2 
distribution.  ’ 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 

Distributions:  Class  given  in  3.2.1,  £ = I . 
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Blumen's  statistic.  For  the  heaviest-tailed  distribution  examined, 

Blumen's  test  appears  to  perform  best  and  MB  performs  considerably 
2 

better  than  T . The  superiority  of  MB  for  very  light-tailed 
distributions  is  shown.  Hodges'  test  statistic,  H,  follows  the  same 
general  pattern  as  Blumen's  statistic,  B,  performing  well  for  heavy- 
tailed distributions,  but  not  well  for  light-tailed  distributions. 
Bennett's  sign  test  and  signed-rank  test  perform  well  except  for 
light-  and  heavy-tailed  distributions,  respectively.  However,  an 
examination  of  Table  4.4  indicates  that  the  performances  of  Bennett's 
nonaffine- invariant  tests  depend  highly  on  the  covariance  structure  of 
the  distribution.  As  the  correlation  decreases  (from  0.00  to  -.99), 
the  empirical  powers  tend  to  decrease,  becoming  rather  low  in  the  most 
extreme  case  in  which  the  correlation  is  -.99.  (The  direction 
that  0 takes  from  the  origin  also  plays  a role  in  the  performances  of 

the  statistics  that  are  not  affine- invar i ant . ) Hence,  we  focus  on  B, 

2 

H,  T , and  MB.  In  Tables  4.5  and  4.6  with  significance  levels  of 
.08865  and  .00881,  respectively,  the  same  general  results  are 
obtained . 

Likewise,  similar  results  are  seen  using  samples  of  sizes  20,  15, 
12,  10,  and  8.  (See  Tables  4.7  through  4.11.)  Since  the  asymptotic 
null  distributions  for  Hodges'  statistic  are  not  suited  for  moderately 
sized  samples,  we  have  used  the  exact  null  distributions  here.  Hodges 
(1955)  gives  the  null  distribution  of  H for  small  samples.  Klotz 
(1964)  gives  percentiles  of  the  null  distribution  of  Blumen's  and 
Hodges'  statistics  for  samples  of  sizes  8 through  12.  We  have 
included  the  available  small  sample  critical  values  for  Blumen's 
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Table  4.4 

One-Sample  Monte  Carlo  Results  with 
n = 25,  a = .05318,  and  Varying  Correlation  Structures 


Statistics 

Bivariate 

medians 

B 

H 

T2  MB 

B1 

B2 

Bivariate  normal 

(unit 

variances 

, correlation  = 

0.0) 

( .000, .000) 

.052 

.045 

.084 

(.054) 

.058 

.047 

.048 

( .1 41 , .141 ) 

.107 

.105 

.181 

(.134) 

.135 

.091 

.119 

(.283, .283) 

.313 

.289 

.468 

(.381) 

.378 

.268 

.361 

( .424, .424) 

.61  2 

.574 

.798 

(.724) 

.708 

.528 

.704 

Bivariate  normal 

(unit 

variances 

, correlation  = 

-0.5) 

( .000, .000) 

.051 

.050 

.077 

(.052) 

.054 

.039 

.054 

(.100,. 100) 

.116 

.100 

.177 

(.136) 

.140 

.087 

.114 

( .200, .200) 

.277 

.261 

.452 

(.368) 

.364 

.184 

• 311 

(.300, .300) 

.604 

.551 

.796 

(.722) 

.716 

.386 

.659 

Bivariate  normal 

(unit 

variances 

, correlation  = 

-0.99) 

( .000, .000) 

.047 

.047 

.075 

(.048) 

.056 

.01  9 

.054 

(.014, .014) 

.107 

.094 

.167 

( .121 ) 

.128 

.022 

.088 

( .028, .028) 

.301 

.270 

.468 

(.390) 

.388 

.022 

.244 

( .042, .042) 

.61  5 

.571 

.774 

( .724) 

.707 

.028 

.470 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are  the 
rejection  proportions  for  Hotelling's  T2  using  the  F?  _2 
distribution.  ,n" 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 

Distributions:  Bivariate  normal. 
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Table  4.5 

Results  with  n = 25  and  a = .08865 

Statistics 

Bivariate 
medians  B 

H 

T2  MB  B1 

B2 

(.000, .000) 

.087 

.050 

.076 

v = .10 
(.041) 

.088 

.113 

.088 

( .01 4, .01 4) 

.203 

.132 

.085 

(.047) 

.137 

.21  9 

.134 

(.028,. 028) 

.368 

.248 

.097 

(.063) 

.190 

.359 

.209 

( .042, .042) 

.525 

.395 

.126 

( .085) 

.266 

.506 

.306 

( .000, .000) 

.087 

.050 

.135 

v = .50 
(.091) 

.114 

.113 

.090 

( .141  , .141  ) 

.224 

.151 

.242 

( .202) 

.224 

.234 

.21  9 

(.283, .283) 

.575 

.414 

.580 

(.524) 

.534 

.533 

.562 

( .424, .424) 

.844 

.747 

.840 

(.799) 

.790 

.802 

.834 

(.000,. 000) 

.087 

.050 

.120 

v = 1.00 
(.085) 

.095 

.113 

.086 

( .141  , .141  ) 

.170 

.106 

.234 

(.188) 

.204 

.173 

.179 

(.283, .283) 

.379 

.276 

.543 

( .478) 

.482 

.368 

.449 

( .424, .424) 

.688 

.556 

.838 

(.793) 

.784 

.614 

.763 

( .000, .000) 

.087 

.050 

.119 

v = 50.00 
( .088) 

.093 

.113 

.088 

(.141, .141) 

.141 

.076 

.244 

(.197) 

.226 

.146 

.168 

(.283, .283) 

.269 

.186 

.510 

(.450) 

.539 

.253 

.405 

( .424, .424) 

.541 

.380 

.833 

(.789) 

.854 

.471 

.722 

(.000,. 000) 

.087 

.050 

.119 

v = 120.00 
(.086) 

.092 

.113 

.081 

(.141, .141) 

.122 

.067 

.227 

( .177) 

.255 

.140 

.171 

(.283, .283) 

.211 

.132 

.494 

(.423) 

.598 

.215 

.391 

( .424, .424) 

.387 

.263 

.812 

(.751) 

.902 

.338 

.679 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are  the 
rejection  proportions  for  Hotelling's  T2  using  the  F2  n_2 
distribution.  ' 

Nominal  significance  level:  .08865. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 

Distributions:  Class  given  in  3.2.1,  1 = I . 


Table  4.6 

One-Sample  Monte  Carlo  Results  with  n = 25  and  a = .00881 


Statistics 

Bivariate 

medians 

B 

H 

T2 

MB 

B1 

B2 

( .000, .000) 

.009 

.002 

.008 

v = .10 
( .004) 

.009 

.004 

.004 

( .01 4, .01 4) 

.030 

.01  4 

.01 1 

(.005) 

.01  4 

.025 

.01  5 

( .028, .028) 

.072 

.045 

.01  8 

(.007) 

.023 

.069 

.029 

( .042, .042) 

.172 

.092 

.020 

(.012) 

.039 

.141 

.063 

(.000, .000) 

.009 

.002 

.017 

v = .50 
(.005) 

.008 

.004 

.005 

(.141, .141) 

.041 

.023 

.071 

( .040) 

.038 

.030 

.041 

(.283, .283) 

.168 

.117 

.258 

(.164) 

.129 

.151 

.145 

( .424, .424) 

.440 

.298 

.586 

(.449) 

.338 

.365 

.428 

( .000, .000) 

.009 

.002 

.025 

v = 1.00 
( .006) 

.004 

.004 

.002 

( .141  , .141  ) 

.029 

.006 

.077 

( .042) 

.024 

.01  6 

.023 

(.283, .283) 

.094 

.043 

.247 

(.162) 

.131 

.072 

.113 

( .424, .424) 

.268 

.182 

.543 

(.401) 

.325 

.21  6 

.330 

(.000, .000) 

.009 

.002 

.026 

v = 50.00 
(.012) 

.010 

.004 

.005 

(.141  , .141) 

.01  8 

.006 

. 066 

( .040) 

.028 

.01  2 

.025 

(.283, .283) 

.051 

.029 

.224 

(.131) 

.129 

.042 

.078 

( .424, .424) 

.158 

.092 

.520 

(.379) 

.402 

.112 

.262 

(.000, .000) 

.009 

.002 

.030 

v = 120.00 
(.010) 

.007 

.004 

.003 

( .141  , .141 ) 

.01  6 

.004 

.062 

( .037) 

.032 

.011 

.01  9 

(.283, .283) 

.034 

.020 

.196 

(.112) 

.156 

.028 

.076 

( .424, .424) 

.088 

.046 

.462 

(.339) 

.501 

.070 

.207 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are  the 
rejection  proportions  for  Hotelling's  T2  using  the  F2  n_2 
distribution. 

Nominal  significance  level:  .00881. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 

Distributions:  Class  given  in  3.2.1,  I = I . 
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Table  4.7 

Results  with  n = 20  and  a = .03044 

Statisti cs 

Bivariate 
medians  B 

H 

T2  MB  B1 

B2 

( .000, .000) 

.023 

.032 

.028 

v = .10 
(.007) 

.030 

.032 

.025 

( .01 4, .01 4) 

.073 

.088 

.034 

(.009) 

.044 

.086 

.038 

( .028,  .028) 

.169 

.184 

.043 

(.018) 

.080 

.180 

.084 

( .042, .042) 

.282 

.287 

.072 

(.033) 

.118 

.285 

.152 

( .000, .000) 

.023 

.032 

.057 

v = .50 

(.023) 

.040 

.032 

.023 

( .141  , .141  ) 

.086 

.090 

.141 

( .082) 

.096 

.077 

.077 

(.283,-283) 

.281 

.265 

.359 

(.255) 

.235 

.247 

.259 

( .424, .424) 

.543 

.545 

.627 

( .524) 

.491 

.503 

.530 

( .000, .000) 

.023 

.032 

.066 

v = 1.00 
(.034) 

.030 

.032 

.023 

( .141  , .141  ) 

. 066 

.065 

.134 

(.068) 

.073 

.056 

.053 

(.283, .283) 

.174 

.174 

.332 

(.232) 

.207 

.162 

.188 

( .424, .424) 

.391 

.379 

.622 

(.508) 

.463 

.342 

.458 

(.000, .000) 

.023 

.032 

.060 

v = 50.00 
(.032) 

.026 

.032 

.023 

(.141, .141) 

.045 

.053 

.115 

(.069) 

.071 

.043 

.045 

(.283,-283) 

.114 

.118 

.311 

(.218) 

.245 

.110 

.159 

( .424, .424) 

.259 

.247 

.608 

( .474) 

.525 

.235 

.378 

(.000, .000) 

.023 

.032 

.066 

v - 120.00 
(.034) 

.033 

.032 

.023 

(.141, .141) 

.043 

.050 

.127 

(.077) 

.105 

.045 

.055 

(.283, .283) 

.089 

.088 

.301 

(.191) 

.296 

.080 

.152 

( .424, .424) 

.175 

.167 

.569 

( .430) 

.597 

.162 

.340 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Hodges'  rejection  proportions 
were  obtained  using  the  exact  small  sample  null  distribution 
of  H.  Values  in  parentheses  are  the  rejection  proportions 
for  Hotelling's  T2  using  the  F2  n_2  distribution. 

Nominal  significance  level:  .03044.  ’ 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  20. 

Distributions:  Class  given  in  3.2.1,  l = I^. 
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Table  4.8 

Results  with  n = 15  and  a = .07050 

Statistics 

Bivariate 
medians  B 

H 

T2  MB  B1 

B2 

(.000, .000) 

.065 

.064 

.065 

v = .10 
(.023) 

.071 

.073 

.046 

( .01  4,  .01  4) 

.122 

.113 

.070 

( .028) 

.093 

. 1 40 

.071 

( .028, .028) 

.21 1 

.207 

.097 

(.040) 

.127 

.237 

.135 

( .042, .042) 

.302 

.299 

.133 

(.063) 

.173 

.345 

.191 

(.000, .000) 

.065 

.064 

.117 

v = .50 
(.058) 

.085 

.073 

.059 

( .141  , .141  ) 

.129 

.130 

.197 

(.116) 

.132 

.149 

.112 

(.283, .283) 

.324 

.320 

.416 

(.294) 

.295 

.305 

.292 

( .424, .424) 

.592 

.548 

.664 

(.540) 

.528 

.562 

.556 

( .000, .000) 

.065 

.064 

.107 

v = 1.00 
(.057) 

.069 

.073 

.057 

( .141  , .141  ) 

.090 

.093 

.164 

(.102) 

.106 

.119 

.094 

(.283, .283) 

.205 

.213 

.363 

(.255) 

.257 

.230 

.233 

( .424, .424) 

.391 

.388 

.628 

(.503) 

.493 

.402 

.460 

(.000, .000) 

.065 

.064 

.116 

v = 50.00 
(.063) 

.086 

.073 

.057 

( .141  , .141  ) 

.075 

.080 

.174 

( .110) 

.134 

.105 

.092 

(.283,-283) 

.164 

.155 

.355 

(.249) 

.290 

.163 

.201 

( .424, .424) 

.303 

.293 

.607 

(.482) 

.557 

.292 

.421 

(.000, .000) 

.065 

.064 

.120 

v = 120.00 
(.065) 

.095 

.073 

.053 

(.141, .141) 

.078 

.073 

.176 

( .106) 

.153 

.094 

.092 

(.283,-283) 

.124 

.129 

.343 

(.243) 

.336 

.137 

.206 

( .424, .424) 

.222 

.222 

.588 

(.446) 

.599 

.21  9 

.384 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Hodges'  rejection  proportions 
were  obtained  using  the  exact  small  sample  null  distribution 
of  H.  Values  in  parentheses  are  the  rejection  proportions 
for  Hotelling’s  T2  using  the  ?2  n_2  distribution. 

Nominal  significance  level:  .07050.  ’ 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  15. 

Distributions:  Class  given  in  3.2.1,  I = I . 
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Table  4.9 

One-Sample  Monte  Carlo  Results  with  n = 12  and  a = .05859 


Statistics 

Bivariate 

medians 

B 

H 

T2 

MB 

B1 

B2 

( .000, .000) 

.050 

(.052) 

.058 

v = 
.062 

.10 

(.020) 

.059 

.049 

.045 

( .01 4, .01 4) 

.086 

( .105) 

.092 

.073 

( .020) 

.078 

.084 

.054 

(.028, .028) 

.1  32 

(.145) 

.147 

.095 

(.040) 

.116 

.149 

.091 

( .042, .042) 

.21  4 

(.230) 

.21  4 

.127 

( .058) 

.144 

.21  7 

.134 

( .000, .000) 

.050 

(.052) 

.058 

v *= 
.119 

.50 

(.055) 

.069 

.049 

.045 

(.141, .141) 

.090 

(.103) 

.104 

.187 

(.102) 

.115 

.079 

.091 

(.283, .283) 

.232 

(.248) 

.249 

.366 

(.219) 

.248 

.202 

.207 

( .424, .424) 

.423 

( .447) 

.440 

.590 

(.431) 

.402 

.388 

.406 

(.000, .000) 

.050 

(.057) 

.058 

v = 
.123 

1.00 

(.059) 

.065 

.049 

.050 

( .141  , .141  ) 

.078 

(.086) 

.084 

.179 

( .088) 

.097 

.081 

.074 

(.283,-283) 

.152 

(.161) 

.156 

.321 

(.195) 

.206 

.144 

.160 

( .424, .424) 

.295 

(.314) 

.31  9 

.542 

(.369) 

.359 

.231 

.311 

(.000, .000) 

.050 

(.056) 

.058 

v = 
.133 

50.00 

(.069) 

.073 

.049 

.054 

( .141  , .141  ) 

.062 

( .072) 

.075 

.182 

(.092) 

.105 

.063 

.064 

(.283, .283) 

.120 

(.126) 

.127 

.318 

(.186) 

.210 

.094 

.140 

( .424, .424) 

.225 

( .241 ) 

.242 

.530 

(.361) 

.383 

.167 

.270 

(.000,. 000) 

.050 

(.058) 

.058 

v = 
.144 

120.00 

(.068) 

. 066 

.049 

.046 

(.141, .141) 

.055 

(.059) 

.065 

.157 

( .092) 

.100 

.059 

.072 

(.283, .283) 

.090 

(.102) 

.099 

.295 

(.164) 

.223 

.079 

.119 

( .424, .424) 

.147 

( .159) 

.167 

.498 

(.311) 

.425 

.126 

.237 

Entries:  Proportion 

of  times 

each 

test  statistic 

exceeded 

its 

asymptotic  critical  value.  Hodges'  rejection  proportions 
were  obtained  using  the  exact  small  sample  null  distribution 
of  H.  Values  in  parentheses  are  the  rejection  proportions 
for  Hotelling's  T and  Blumen's  test  using  the  F2  n_2 
distribution  and  the  exact  small  sample  null  distribution 
for  B,  respectively.  Blumen's  test  was  randomized  when 
using  the  small  sample  distribution. 

Nominal  significance  level:  .05859. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  12. 

Distributions:  Class  given  in  3.2.1,  E = I . 
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Table  H . 1 0 

One-Sample  Monte  Carlo  Results  with  n = 10  and  a = .03906 


Statistics 

Bivariate 

medians 

B 

H 

T2 

MB 

B1 

B2 

(.000, .000) 

.039 

(.039) 

.045 

v = 
.045 

.10 

(.007) 

.039 

.051 

.013 

( .01 4, .01 4) 

.067 

(.067) 

.073 

.058 

(.011 ) 

.061 

.076 

.024 

( .028, .028) 

.106 

(.106) 

.104 

.075 

(.019) 

.072 

.1  21 

.043 

( .042, .042) 

.164 

(.164) 

.151 

.110 

(.033) 

.109 

.163 

.072 

(.000, .000) 

.039 

(.039) 

.041 

v = 
.111 

.50 

(.027) 

.050 

.051 

.017 

(.141, .141) 

.068 

(.068) 

.068 

.158 

(.064) 

.071 

.083 

.043 

(.283, .283) 

.155 

(.155) 

.139 

.301 

(.153) 

.150 

.1  61 

.108 

( .424, .424) 

.294 

(.294) 

.257 

.500 

( .274) 

.241 

.283 

.229 

(.000, .000) 

.039 

(.039) 

.046 

v = 
.118 

1 .00 
(.048) 

.052 

.051 

.028 

( .141  , .141 ) 

.064 

( .064) 

.053 

.149 

( .059) 

.071 

.061 

.040 

(.283, .283) 

.112 

(.112) 

.108 

.272 

(.121) 

.115 

.110 

.077 

( .424, .424) 

.197 

(.197) 

.162 

.455 

( .242) 

.197 

.187 

.149 

( .000, .000) 

.039 

(.039) 

.042 

v = 
.129 

50.00 

(.052) 

.050 

.051 

.026 

( .141  , .141  ) 

.052 

( .052) 

.055 

.147 

(.067) 

.062 

.058 

.038 

(.283, .283) 

.087 

(.087) 

.079 

.246 

(.111) 

.106 

.086 

.063 

( .424, .424) 

.133 

(.133) 

.121 

.410 

( .202) 

.192 

.142 

.118 

( .000, .000) 

.039 

(.039) 

.047 

v = 

.126 

120.00 

(.050) 

.053 

.051 

.025 

( .141  , .141 ) 

.049 

(.049) 

.053 

.150 

(.066) 

.079 

.053 

.038 

(.283, .283) 

.072 

(.072) 

.059 

.231 

(.117) 

.132 

.071 

.060 

( .424, .424) 

.104 

(.104) 

.090 

.385 

(.188) 

.232 

.108 

.109 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Hodges’  results  based  on  a 
randomized  test  using  the  exact  small  sample  null 
distribution.  Values  in  parentheses  are  the  rejection 
proportions  for  Hotelling’s  T2  and  Blumen's  test  using  the 
F2  n_2  distribution  and  the  exact  small  sample  null 
distribution  for  B,  respectively. 

Nominal  significance  level:  .03906. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  10. 

Distributions:  Class  given  in  3.2.1,  l = I . 
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Table  *1.11 

One-Sample  Monte  Carlo  Results  with  n = 8 and  a = .06250 


Statistics 

Bivariate 

medians 

B 

H 

T2 

MB 

B1 

B2 

(.000, .000) 

.068 

(.068) 

.068 

v = 
.108 

.10 

(.024) 

.090 

.017 

.032 

( .01  4, .01 4) 

.086 

( .086) 

.086 

.118 

( .025) 

.098 

.028 

.045 

( .028, .028) 

.133 

(.133) 

.133 

.138 

(.034) 

.136 

.053 

.076 

( .042, .042) 

.167 

(.167) 

.167 

.160 

( .057) 

.152 

.085 

.114 

(.000, .000) 

.068 

(.068) 

.068 

v = 
.167 

.50 

(.051) 

.089 

.017 

.042 

(.141, .141) 

.095 

(.095) 

.095 

.226 

(.066) 

.114 

.032 

.064 

(.283,-283) 

.181 

(.181) 

.1  81 

.366 

(.153) 

.179 

.068 

.135 

( .424, .424) 

.306 

( .306) 

.306 

.525 

(.294) 

.262 

.156 

.259 

( .000, .000) 

.068 

(.068) 

.068 

v = 
.180 

1 .00 
(.075) 

.088 

.017 

.050 

( .141  , .141  ) 

.084 

(.084) 

.084 

.21  8 

( .081 ) 

.104 

.021 

.058 

(.283, .283) 

.137 

(.137) 

.137 

.317 

(.136) 

.171 

.044 

.095 

( .424, .424) 

.231 

(.231) 

.231 

.482 

( .243) 

.251 

.089 

.191 

(.000, .000) 

.068 

(.068) 

.068 

v = 
.183 

50.00 

(.079) 

.092 

.017 

.045 

( .141  , .141 ) 

.079 

(.079) 

.079 

.226 

(.091) 

.110 

.01  9 

.060 

(.283,-283) 

.115 

(.115) 

.115 

.313 

(.138) 

.158 

.036 

.095 

( .424, .424) 

.187 

(.187) 

.187 

.440 

(.219) 

.235 

.080 

.152 

(.000, .000) 

.068 

(.068) 

.068 

v = 
.178 

120.00 

(.076) 

.092 

.017 

.055 

( .141  , .141  ) 

.072 

(.072) 

.072 

.21  8 

( .094) 

.114 

.01  8 

.064 

(.283,-283) 

.093 

(.093) 

.093 

.313 

(.134) 

.166 

.024 

.082 

( .424, .424) 

.142 

(.142) 

.142 

.407 

(.214) 

.263 

.046 

.128 

Entries:  Proportion 

of  times 

each 

test  statistic 

exceeded 

its 

asymptotic  critical  value.  Hodges'  rejection  proportions 
were  obtained  using  the  exact  small  sample  null  distribution 
of  H.  Values  in  parentheses  are  the  rejection  proportions 
for  Hotelling's  T and  Blumen's  test  using  the  ?2  n_2 
distribution  and  the  exact  small  sample  null  distribution 
for  B,  respectively. 

Nominal  significance  level:  .06250. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  8. 

Distributions:  Class  given  in  3.2.1,  E = I . 
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statistic  in  our  Monte  Carlo  study  for  samples  of  size  8,  10,  and 
12.  The  significance  levels  selected  were  chosen  from  the  natural 
levels  of  Blumen's  and/or  Hodges'  null  distribution.  In  some  cases, 
either  Blumen's  test  or  Hodges'  test  was  randomized  so  that  the 
significance  levels  would  agree.  (There  are  relatively  few  natural 
significance  levels  for  Hodges’  statistic  for  small  samples.)  An 
examination  of  the  empirical  significance  levels  obtained  suggests 
that  samples  of  size  10  or  larger  are  needed  to  effectively  use  the 
asymptotic  null  distribution  of  MB. 

In  Tables  4.12  we  present  the  empirical  powers  for  samples  from 

light-tailed  Pearson  Type  II  distributions  and  heavy-tailed  Pearson 

Type  VII  distributions.  Again,  the  same  general  trends  are 

indicated.  Finally  in  Table  4.13  we  examine  samples  from  bivariate 

normal  mixtures  violating  the  assumption  of  elliptical  symmetry 

required  for  the  asymptotic  results  of  Chapter  3.  Nevertheless,  MB 

appears  to  be  rather  robust.  The  contours  of  the  first  two  mixtures 

are  illustrated  on  page  56  of  Johnson  (1987).  (The  parameters  of  the 

distributions  and  the  mixing  probabilities  are  indicated  in  Table  4.13 

for  each  mixture.  The  means  have  been  translated  from  that  shown  by 

Johnson  so  that  the  bivariate  null  mean  of  the  mixture  is 

0.)  For  these  asymmetric  distributions,  MB  performs  similarly  to 
2 

Hotelling's  T and  better  than  the  other  competitors.  The  bimodal 
distribution  used  in  Mixture  3 (Johnson,  1987,  p.  61  ) results  in  a 
relatively  large  significance  level  for  MB,  but  significantly  greater 
power  compared  to  the  power  of  its  competitors.  The  last  two  mixtures 
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Table  4.12 

One-Sample  Monte  Carlo  Results  for  Pearson  Type  II 
and  Type  VII  Distributions  with  n = 25  and  a = .05318 


Statistics 

Bivariate 

medians 

B 

H 

T2 

MB 

B1 

B2 

(.000,. 000) 

.056 

.043 

Bivariate  Normal 
.088  (.062)  .067 

.056 

.051 

( .141  , .141  ) 

.119 

.108 

.187 

(.128) 

.135 

.103 

.119 

(.283, .283) 

.326 

.295 

.478 

(.395) 

.394 

.248 

.363 

( .424, .424) 

.600 

.559 

.794 

(.727) 

.702 

.508 

.688 

( .000, .000) 

.056 

.043 

Pearson 

.096 

Type  II 
( .065) 

(m  = 1) 
.068 

.056 

.049 

( .1  41 , .1 41) 

.092 

.080 

.1  81 

(.131) 

.153 

.084 

.109 

(.283, .283) 

.221 

.205 

.463 

(.387) 

.445 

.186 

• 327 

(.424, .424) 

.441 

.405 

.778 

(.706) 

.787 

.344 

.632 

(.000, .000) 

.056 

.043 

Pearson 

.099 

Type  II 
(.063) 

(m  = .5) 
.068 

.056 

.051 

( .141  , .141  ) 

.085 

.075 

.173 

(.128) 

.197 

.073 

.112 

(.283, .283) 

.195 

.176 

.458 

(.377) 

.545 

.163 

.325 

( .424, .424) 

.379 

• 353 

.780 

(.710) 

.840 

.294 

.641 

( .000, .000) 

.056 

.043 

Pearson 

.068 

Type  VII 
( .047) 

(m  = 2) 
.068 

.056 

.055 

( .1  41 , .1  41) 

.167 

.154 

.160 

(.112) 

.155 

.138 

.146 

(.283, .283) 

.476 

.444 

.388 

(.325) 

.411 

.401 

.439 

(.424, .424) 

.791 

.750 

.624 

(.561) 

.666 

.713 

.741 

(.000,. 000) 

.056 

Pearson  Type 
.043  .046 

■ VII  (m 
(.025) 

= .5,  Cauchy) 
.071  .059 

.056 

( .141  , .141  ) 

.294 

.274 

.113 

(.069) 

.212 

.184 

.238 

(.283, .283) 

.712 

.669 

.273 

(.218) 

.491 

.463 

.650 

( .424, .424) 

• 933 

.910 

.443 

(.382) 

.708 

.671 

.899 

Entries:  Proportion  of  times  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are  the 
rejection  proportions  for  Hotelling’s  T2  using  the  ?2  n_2 
distribution. 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 
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Table  4.12 — continued. 


Distributions : 

Bivariate  Pearson  Type  II 


T (m+2) 
r (m+1  )tt 


e|  1 /2{1-(y-0)  'Z  Vy-e)}"1  if  (y-0) 'Z_1 (y-0)  £ 1,  m > -1. 


Bivariate  Pearson  Type  VII 


T (m) 

r(m-1  )tt 


|E|"1/2{1+(y-0) 'Z-1 (y-0) }"m, 


m > 1 , y e R 
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Table  4.13 

One-Sample  Monte  Carlo  Results  for  Bivariate  Normal  Mixtures 
with  n = 25  and  a = .05318 


Statistics 

Bivariate 

medians  B H T2  MB  B1  B2 


Mixture 

1 (P  = .5, 

y11  = 

“12  ' ■’•0l  °11 

11 

OJ 

O 

II 

• 0, 

P1  “ 

o.o,  p21  = 

U22  = 

1.°.  o2,  - »22 

C\l 

Ql 

O 

II 

= 0.0) 

( .000, .000) 

.054 

.051 

.082 

( .060) 

.054 

.050 

.051 

( .1  41 , .1 41) 

.060 

.061 

.115 

(.078) 

.088 

.057 

.071 

( .283, .283) 

.099 

.090 

.21  2 

(.152) 

.165 

.093 

.133 

(.424, .424) 

.170 

.149 

.356 

(.290) 

.293 

.158 

.256 

Mixture  2 (p  = 

.5,  yn 

y12 

-.5,  on 

= °12  = 1 

.0, 

p1  = o.o,  u21 

y22  = 

.5,  o21 

°22  = 1 

•0,  p2  = 

.9) 

( .000, .000) 

.054  .051 

.073 

(.043) 

.059 

.050 

.071 

( .1  41 , .1  41) 

.096  .085 

.140 

(.089) 

.099 

.076 

.083 

( .283, .283) 

.181  .172 

.278 

(.219) 

.231 

.163 

.207 

( .424, .424) 

.345  .311 

.507 

(.414) 

.409 

.300 

.392 

Mixture 

3 (p  = 

.5,  u,, 

y12  " 

-1.0,  on 

°1  2 = 

1.0, 

P1  " " 

.9,  u21 

P 22  = 

1.0,  o21 

°22  = 1 

CM 

CL 

o 

= 0.0) 

( .000, .000) 

.050 

.045 

.089 

( .062) 

.086 

.050 

.092 

( .1  41 , .1  41) 

.058 

.068 

.106 

(.070) 

.110 

.051 

.060 

(.283, .283) 

.086 

.106 

.199 

( .157) 

.235 

.115 

.135 

(.424, .424) 

.136 

.168 

.390 

(.298) 

.424 

.191 

.268 

Mixture  4 (p  = 

.9, 

P22  " 

0.0,  on 

= °12  = 1 

• 0, 

p1  >=  o.o,  u21 

‘ V22  " 

0.0,  o21 

°22 

20.0,  p2 

= 0.0) 

( .000, .000) 

.037  .039 

.022 

(.012) 

.044 

.035 

.033 

(.141, .141) 

.092  .085 

.050 

(.024) 

.102 

.067 

.083 

( .283, .283) 

.261  .232 

.120 

( .088) 

.260 

.198 

.250 

(.424, .424) 

.500  .466 

.231 

(.193) 

.494 

.402 

.469 
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Table  4.13 — continued. 


Bivariate 

medians 

Statistics 

B 

H 

T2 

MB 

B1 

B2 

Mixture  5 (p  = 

.9,  - »,2  - 0. 

0, 

II 

OJ 

o 

II 

.0, 

P1  = 

o.o,  y21 

" U22  ' °-°'  °21  ' 

ro 

o 

o 

o 

^r 

ii 

00 

00 

o 

0, 

P2  = 

0.0) 

( .000, .000) 

.037 

.039 

.018  (.008) 

.051 

.035 

.032 

( .1  41 , .1  41) 

.087 

.088 

.040  (.021) 

.097 

.067 

.081 

(.283, .283) 

.251 

.230 

.107  (.079) 

.268 

.197 

.248 

(.424, .424) 

.502 

.467 

.191  (.157) 

.485 

.401 

.468 

Entries: 


Proportion  of  times  each  test  statistic  exceeded  its 
asymptotic  critical  value.  Values  in  parentheses  are  the 


using  the  F 


2,n-2 


rejection  proportions  for  Hotelling's  T 
distribution. 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  size:  25. 

Distributions:  Bivariate  normal  mixture,  choosing  N1  with  probability 

p and  N2  with  probability  1-p,  where  Nj»  i - 1 , 2,  is 
bivariate  normal  with  null  mean  (u 


matrix 


i 1 ’ 


ui2)  and  covariance 


i 1 


pi°i1°i2 


P i°i 1 °i 2 


12 
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are  heavy-tailed  distributions.  In  each  case  a bivariate  normal 

random  variable  with  independent  standard  normal  components  is 

contaminated  with  a bivariate  normal  variable  with  independent 

components,  each  having  a mean  equal  to  zero  but  a relatively  large 

variance.  Again  the  superiority  of  MB  and  the  other  competitors  over 
2 

Hotelling’s  T is  demonstrated  for  heavy-tailed  distributions. 

Lastly,  we  note  that  we  attempted  to  approximate  the  null 
distribution  of  MB  for  small  samples  using  the  null  distribution  of 

* -1  *2  *2 

m ' 6"  (CHB+SMB>’ 

where 


* 1 n uRankU  ) # 

CMB  = - I aicos( - — — )(Rank(Di)), 

i=1 

* n nRank(ij) . ) # 

SMB  = n 1 aisin( n — — ) (Rank(Di ) ) , 


» * 

and  is  as  in  Chapter  3.  That  is,  we  replaced  by  D in  the 

definition  of  the  test  statistic  (see  2.1.5  and  3.1.2).  The 

* 

independence  of  4^  and  under  the  null  hypothesis  allows  us  to 

* 

obtain  the  null  distribution  of  MB  by  generating  all  possible 

* # 

assignments  of  ranks  to  4^  , . . . ,4^  and  D^...^  . For  larger  samples, 
this  is  quite  computer- intensive , so  we  instead  randomly  selected 

100,000  rank  assignments  from  all  possible  assignments  of  the  ranks 

* * 

to  4>,,...,4>  and  D,  ,...,D  . Monte  Carlo  simulations  were  then 

1 n 1 n 

performed  using  the  approximate  percentiles  obtained.  The  results 
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were  inconclusive.  In  many  cases  these  small  sample  percentiles 
appeared  closer  to  the  true  percentiles  than  that  obtained  from  the 
asymptotic  distribution;  however,  this  was  not  always  true.  In 
general,  though,  for  moderately  large  samples,  namely  n £ 10,  the 
large  sample  percentiles  appear  (in  Tables  4.3  - 4.13)  to  be 
reasonably  close  to  the  true  percentiles  of  the  null  distribution. 

In  summary,  of  the  affine-invariant  statistics,  Blumen’s  test  has 
the  greatest  power  for  heavy-tailed  distributions  while  MB  has  the 
greatest  power  for  light-tailed  distributions.  Unlike  Blumen's 
statistic  and  Hodges’  statistic,  MB  performs  very  well  relative  to 
Hotelling's  T2  under  bivariate  normality.  The  only  setting  for  which 
Blumen's  statistic  and  Hodges'  statistic  perform  better  than 
Hotelling's  T2  is  when  the  underlying  distribution  is  heavy-tailed. 

On  the  other  hand,  MB  performs  better  than  Hotelling's  T2  for  both 
light-  and  heavy-tailed  distributions.  Further,  MB  appears  robust 
with  respect  to  moderate  departures  in  the  assumption  of  elliptical 
symmetry.  The  wide  range  of  distributions  for  which  MB  has  relatively 
high  power  suggests  that  MB  is  a promising  procedure  for  the  one- 
sample  bivariate  location  problem. 


CHAPTER  5 

A BIVARIATE  TWO-SAMPLE  TEST  FOR  LOCATION 


5. 1 Definition  of  the  Test  Statistic 
A natural  extension  of  the  one-sample  statistic,  MB,  introduced 
in  Chapter  2,  will  be  described  in  this  section.  Assume  that 
Xl’’'',Xm  3005  *1’*’’,Yn  arS  indePendent  random  samples  from  bivariate 
populations  with  continuous  distribution  functions  ( u ^ , u^)  and 

^2^U1  ’ U2^  = U2~®2^  ’ resPecti vely • We  are  interested  in 

testing 


H : 0 = (0  , 0 ) = 0 against  H : 0 = (0,,  01'  * 0.  (5.1.1) 

o i d a „ 1 2 _ 

We  mentioned  in  the  introduction  that  Mardia’s  two-sample 
statistic  extends  Blumen's  work  to  two  samples.  Mardia  uses  the 
center  of  gravity  of  the  standardized  versions  of  the  centered  vectors 
from  the  first  sample.  Specifically,  define 


Z 


m+n 


m 

I 

i = 1 


n 


X_i 


+ I 
j = 1 


(5.1.2) 


the  mean  vector  for  the  combined  sample.  Let  4>  . (<(>  ) denote  the 

X 9 1 I f J 

angle  between  X^Z  (Yj-Z)  and  the  positive  half  of  the  horizontal  axis, 
where  0 £ ^ ^ j)  < 2ir,  i = 1,...,m,  j = 1,...,n.  The  center  of 

gravity  (C.,  Su)  of  Mardia's  statistic  is  defined  by 

MM 
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98 


1 


m 


CM  ' S Z 003<' 
1=1 


2TTRank(*  ) 

^ f 1 


•) 


and 


(5.1.3) 


. m 2irRank(({)  .) 

SM  = ~ E sin( X>1 

Mm  N 


), 


where  N = m+n,  and  Rank(<|>  ) denotes  the  rank  of  4>  among 

^ 9 X j i 

♦x,i'”*'*x  m’^Y  1 n’  The  sign  statistic  of  Mardia  is  then 


M = 2 S (n_1){c2+s2} 

n MM 


(5.1 .H) 


The  signed-rank  statistic  we  propose  uses  the  distances  of  the 
centered  observations  from  the  origin,  in  addition  to  the  directions 
of  the  vectors  from  the  origin.  Define 

- 1 m n 

Z = -{  I (X  -Z)(X  -Z)'  + I ( Y -Z) ( Y . -Z ) ' } , 
i-1  ~ ~ -1  ~ j-1  ~ ~J  ~ 


,2-1 


DX,i  = (Xi~Z)'E  (xi-z).  i = 1 m, 


(5.1.5) 


Dy>J  = ( Yj-Z)  'E_1  (Yj-Z) , j = 1 n, 


and  let  Rank(D  ) (Rank(D  .))  denote  the  rank  of  Dv  . (D  .)  among 

1 »J  X , l Y , j 


°X,r  '••'DX,m,DYf1 DY,n*  Let 


m 


2irRank(<t>  ) Rank(Dv  .) 


'MB2 


- E cos(- 
m i-1 


X,i 


N 


-)(- 


XjlL 


N 


) 


n 


..  2irRank(<{)  ) Rank(D  .) 

E cos( — )(  ,J  % 

n j=1  N N ; 
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and 


(5.1.6) 


1 


m 


2irRank(<j>  ) Rank(Dv  .) 


SMDO  = — I sin(- 
MB2  m 


rX,i 


N 


-)(- 


X,i 


N 


•) 


1 


n 


- — I sin( 

J-1 


2irRank(<|>  .)  Rank(D  .) 

1^-)( _I^_) 


N 


N 


The  proposed  two-sample  statistic  is 


MB2  . 


(5.1.7) 


Notice  that,  unlike  the  components  of  Mardia's  statistic  (in  5.1.3), 

CMB2  30(1  SMB2  h3Ve  6aCh  been  constructed  t0  have  an  asymptotic  mean  of 
zero  under  HQ.  This  centering  is  useful  when  developing  asymptotic 
properties  of  MB2.  Likewise,  it  will  be  convenient  to  examine  a 
"centered"  version  of  Mardia's  statistic.  Let 


1 


m 2irRank(  <f>  ) n 2irRank(<t>v  .) 

CM2  • s °°3< N ' ^ - 5 «•< tr^ 


), 


m 


- 2wRank(*  ) n 2irRank(<fc  ,) 

^ ■ s sl"< ir^)  - n jf,  31"< ir^). 


and 


(5.1.8) 


1 1-12  2 
M2  = 2(—  + 1)  '{cl  + s"  }. 
m n M2  M2 


This  centered  modification,  M2,  of  Mardia's  two-sample  statistic  will 
be  investigated  and  compared  to  MB 2 in  the  following  chapters. 
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5.2  Affine-Invariance  of  MB2 
Like  MB,  the  statistic  MB2  is  invariant  with  respect  to 
nonsingular  linear  transformations  of  the  observations.  In  verifying 
this,  we  use  the  fact  that  for  all  i and  j, 


A j A 

(LX  -LZ)'(LEL'f  \lx  -LZ)  = (X  -Z)  'E_1  (X. -Z)  = Dy 


and 


(5.2.1) 


A j A 

(LY  -LZ)'(LEL')~  (LY  -LZ)  = (Y  -Z) 'E_1 ( Y -Z)  = Dv  . , 
~ ~ —J  ~ ~J  — ± »J 


for  any  2x2  nonsingular  matrix,  L.  That  is,  the  standardized 

distances  D ,...,D  and  D ,...,D  are  themselves  invariant 
X,1  X,m  Y, 1 Y,n 

with  respect  to  nonsingular  linear  transformations  of  the 
observations.  Again,  we  prove  the  invariance  of  MB 2 through  a 
sequence  of  three  lemmas. 

Lemma  5.2.2.  The  statistic,  MB2,  is  invariant  with  respect  to  a 
rotation  of  all  of  the  observations. 

Proof.  Let  X X_  , Y, Y_  , 


_m  _1 


’ _n ’ TX,1 


"Vm*  ♦y.I *Y,n’ 

3nd  DX,1 Dx  m’  DY  1 °Y  n 1)6  as  already  defined.  Now  rotate 

the  observations  clockwise,  yielding  the  transformed  vectors, 

# # * * 

= PX1 ^ = PXm,  and  Y1  = PY1 Yn  = PYn,  where  P is  a 2x2 

orthogonal  matrix  of  the  type  discussed  in  the  proof  of  Theorem 


-# 

Z 


{ 

i=1 


* 

Z Y }, 
J-1 


2. 2. 7.  If  we  let 
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£ £ 

then  Z = PZ.  Hence  if  we  find  the  angles  <j>  * and 

~ X , 1 X ,m 

* * , » * * * 

4>y  i » - • • » 4>  Y n associated  with  X , . . . ,X  and  Y ,...,Y  to  construct 
1 > 1 if*1  ~ i _m  _ 1 _n 

# * 

MB2,  then  <j>  ( ) represents  the  angle  between  P(X  -Z)  (P(Y.-Z)) 

A » 1 1 > J 

and  the  positive  half  of  the  horizontal  axis.  The  vectors 

P(X  -Z) , . . . , P(x  -Z)  and  P(Y  -Z) , . . . ,P(Y  -Z)  are  obtained  from 
- 1 ~ - ~ i ~ _m  _ 

Xi"Z Xm-Z  and  Y^-Z Y ' - Z,  respectively,  by  rotating  each  vector 

clockwise  via  P.  Assume  that  only  one  centered  vector  is  rotated 

across  the  horizontal  axis.  Without  loss  in  generality,  assume  that 

this  vector  is  X, -Z.  Then 
_k  _ 


Rank(<f>x  ^ = 


|Rank(<f>x  ^-1, 


if  i * k, 


if  i = k , 


and 


Rank((j)y  ) 


Rank(4>Y  ) - 1,  for  all  j . 


We  now  construct  MB2  from  the  rotated  observations,  yielding 


m 


MB2 


- I cos(- 
m i-1 


2irRank(<|i  ) Rank(Dv  . ) 


N 


JL± 


-)(- 


X , i 


N 


■) 


1 


n 


I cos( 

n j=1 


2irRank(<f>  ) Rank(D  .) 
i .J  ) ( 


N 


N 


), 
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m 


2irRank(<}>  ) Rank(Dv  ) 


MB2 


- Z sin(- 
m i-1 


Xi 


-)(- 


X,i 


N 


•) 


1 


n 


2irRank(<)>  •)  Rank(Dv  .) 


Z sin(- 

" J-1 


Xi 


N 


-)(- 


Id. 


N 


■), 


and 


1 1 -1  * * 

MB2  = 6(-  + -)  {C  + S 

m n 1 MB2  MB2  ’ 


using  5.1.6,  5.1.7,  and  5.2.1.  Noting  that 


2TiRankU  ) 2ir(Rank(*  J-1) 

COS  ( ;; ’■ ) = COS(— jj—  ) = cos(0)  = COS  ( - — 1 ), 


N 


and 


sin(- 


2irRank((f>x  J 
N 


) = sin(- 


2ir(Rank(<f>x  J-1) 
N 


We  can  replace  Rank(*  ) and  RankU  ) by  Rank(4>v  . )-1  and 

I f J X , 1 

Rank(4>x  J-1,  respectively,  for  all  i and  j in  the  expressions  for 


CMB2  antl  W yleldi"8 


m 


'MB2 


1 r.  , 

- Z cos  (- 
m i=1 


2ir(Rank(i(i  )-1)  Rank(D„  J 


X , i 


N 


-)(- 


X,i 


n 


- Z cos ( 

J-1 


2ir  ( R an  k ( <f>  J-1)  Rank(D„  J 

^ )( 


N 


N 


and 


(5.2.3) 
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SM2 


1_ 

m 


m 

Z 

i=1 


2ir(Rank(<{>  )-1 ) 
sin( Jbl )( 


Rank (D  ) 

X > 1 N 

N ’ 


n 

Z 

j = 1 


2ir(Rank(4>  )■ 
sin( — 


1) 


N 


-)( 


Rank(D  .) 

hi) 

N ; ‘ 


Using  the  trigonometric  identities  given  in  2.2.4,  we  can  write 


"MB  2 


m 2irRank(4>  ) 2irRank(<|>v  ) _ Rank(D  ) 

S^'003' N — ‘ ^)cos(-J)*sin( 


n 


2irRank(ij)  .) 


1 ” ‘i-v  4/  ?_  2nRank(<fi  ) Rank(Dv  .) 

1 z { cos  ( n — jJ-)cos(^)+sin( --Y-t.J..)sin(^) } 


n 


j = 1 


N 


and 


m 


MB2 


i1_ysin(- 


2irRank(<j>  ) 


2irRank(<j)  . ) 


N 


y ■] ' i-nnuumy  . / _ Rank(D  . ) 

)cos(^)-cos( rAi_)sin(2l)  } 


n 


2iTRank(<t>  ) 


- - Z {sin( 

j = 1 


2irRank(<|>  ) 


1 ' )cos  (^)-cos  ( 


N 


N 


)sin(— ) }- 


Rank(D  ) 

— hll 

N 


That  is, 


4, 


cos (— — )C 

IT  MB  2 


+ sin(^)S,m„ 
N MB2 


and 


4a 


003(^)SMB2 


sin(— )C 


MB2 ' 


Hence, 
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CMB2+SMB2 


(0os2(2l).sln2(2i)){c^2.s^2) 


2 2 
= C +S 
MB2  MB2* 

As  in  the  one-sample  setting,  the  invariance  with  respect  to  counter- 
clockwise rotations  and  further  clockwise  rotations  follows. 

Lemma  5.2.4.  The  statistic,  MB2,  is  invariant  with  respect  to 

reflections  of  the  observations  about  the  horizontal  axis. 

* * * # 

Proof.  Let  X ,,..., X and  Y ,Y  be  the  reflected  vectors. 

-1  ~m  _1  _n 

That  is. 


and 


where 


Then 


Xi  RX^ , i 1 , . . . ,m , 


Yj  = RYj,  j = 1 n, 


R = 


0 -1 


-#  i ra  * n * 

z = =•  { Z X + Z Y } = RZ, 
i-1  ~ J-1 


# -* 


Xj-Z  = R(Xi-Z), 


and 


Y*-Z*  = R(Y  -Z). 
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So,  the  centered  vectors  are  also  reflected  about  the  horizontal 

# # — # # _ # 

axis.  Letting  <f>  (<j>  ) denote  the  angle  between  X -Z  (Y  -Z  ) and 

~i  ~ ~j  .. 

the  positive  half  of  the  horizontal  axis,  it  follows  that  the  ordering 

^ & $ $ 

of  the  angles  <p  <p  , <p  is  the  reverse  of  the 

a , i a ,m  i,i  i ,n 

ordering  of  the  angles  <J>  <p  <f>v  . That  is, 

A » 1 a i,i  i ,n 


RankU  ) = N+1 -RankU  ) , 

A t 1 X , 1 

and 


RankU  ) = N+1-RankU  ), 
1 >J  * *J 


for  all  i and  j.  As  in  5.2.3,  we  compute  C*B2  and  S*g2  from  the 
reflected  observations  and  find  that 


"'MB  2 


1 


m 


= — I cos 
mi  = 1 


(2irRank  ( 41  )URank(D  )l  n i 2 

k-f 


irRank(4)  .)  l,Rank(Dv.  .) 


+ > J 


N 


m 


1 i eojiimu  . 2l,lian'<(*X.l)lRimk(Dx.i 


N 


n 


f 


)|Rank<°x.i)j 


1 i 003|M#ii) _ ^k<*Y,.1)ira"k(Di.i) 


N 


N 


N 


J 


cos<^^)c_.sin(2^)s. 


N ' MB2 


N MB2  ’ 


hi 


N 


;J 


using  the  identities  in  2.2.4.  Also 
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MB2 


_ i-2it(  N+1 ) ,2ir(N+1)._ 

31n(— T~ : icmb2-°03(— >W 


So , 


*2  *2 
CMB2+SMB2 


c2  +S2 

MB 2 MB2' 


□ 


Lemma  5.2.5.  The  statistic,  MB2,  is  invariant  with  respect  to 
changes  in  scale. 

Proof.  Apply  the  proof  of  Lemma  2.2.6  to  the  vectors 

^ 1 * • • • » X^-Z , Y^-Z, . . . , Y^-Z . 

Theorem  5.2.6.  The  statistic,  MB2,  is  invariant  with  respect  to 
nonsingular  linear  transformations  of  the  observations. 

Proof.  Using  Lemmas  5.2.2,  5.2.4,  and  5.2.5,  the  proof  follows 
the  pattern  of  Theorem  2.2.7.  □ 

It  is  also  important  for  the  two-sample  problem  to  show  that  MB2 
is  invariant  with  respect  to  translations  of  the  observations.  That 
is, 


MB2(X1 


,Yn)  = MB2(X*, 


), 


where  X*  = X.+d,  Y*  = Yj+d  for  all  i and  j,  and  d = (d, , d ) 
constant  vector.  Note  that  MB2  is  a function  only  of  X^-Z,. 


is  a 

.X-Z, 
_m  _ 


-* 

Z 


{ 

i=1 


n 

+ I 

J-1 


} 


= Z+d, 
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x*-z*  = (Xi-d)  - (z+d)  = x -z, 


and 


* _# 

Tfz- 

for  all  i and  j.  Hence  the  translation  invariance  follows. 

Lemma  5.2.7.  The  statistic,  MB2,  is  invariant  with  respect  to 
translations  of  the  observations. 

Slight  modifications  in  the  proofs  of  this  section  will 
demonstrate  that  the  two  versions  of  Mardia's  statistic  (given  by 
5.1.^  and  5.1.8)  are  also  affine-invariant. 


CHAPTER  6 

ASYMPTOTICS  FOR  THE  TWO-SAMPLE  PROBLEM 


6.1  Asymptotic  Null  Distribution  of  MB2 


In  this  section,  we  find  the  asymptotic  null  distribution  of  MB2 
under  the  assumption  that  the  underlying  populations  are  elliptically 
symmetric  (see  3.1. 11)).  As  in  the  one-sample  case,  we  introduce  a 
statistic  which  has  the  same  asymptotic  null  distribution  as  MB2.  Let 

^1'***’jn  *ncJePendent  random  samples  from  continuous 

elliptically  symmetric  bivariate  distribution  functions,  G1 , u ) 
and  G2(u1  , u2)  = G^u^O^,  U2~Q2^’  resPecti velY • Note,  under 


d 


H0:  ® _ (®i>  Xi  = Yj  » for  all  i and  j.  Let 


m, 


and 


(6.1.1) 


where 


and 


(6.1.2) 


o 
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Let 


F2(t)  PH  (DX,i  * t}  " PH  (DY,j  “ t)f 


(6. 1 .3) 


and  for  i = 1 , . . . ,m  and  j = 1 , . . . ,n  define  a (<j>  ) to  be  the  angle 

x 1 1 Y , j 

between  X^y  (Y  -y)  and  the  positive  half  of  the  horizontal  axis, 


0 S *x  i <+y  j)  < 2ir.  Let 


H2(t)  ' PH  <*X,1  S ° ' PH  (»Y,j  S » 


(6.1.4) 


Finally,  define  two  empirical  distribution  functions 


1 m n 

F (t)  = rr  { l I H D < t]  + I I[Dy  < t]}, 
1=1  A>1  J-1 


and 


(6.1.5) 


H (t)  = ^ { z IC 4>x  i s t]  ♦ E I[*  < t]}, 

i=1  X>1  j = l Y’J 

where  D^,  Dy>J,  ^ and  <f>Y  j are  as  given  in  Chapter  5. 
Using  5.1.6,  we  can  write 


'MB2 


1 m , n 

5 jf,  °°3<2,'"n(*X,1),FN(DX.1)  ~ K jf,  00S(2’'HN<*Y,j)FN(DY,j) 


and 


(6.1.6) 


m 


^MB2  m 2 3 i n ( 2ttHm ( <f>v  ,))FH(DV  4)  - 


i=1 


N ^X , i FN  °X , i 5 sin(2*V  *Y,  j > } FN(DY , j > ' 
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Define 


"MB2 


1 m * * 

5 1f,OOS(2'H2<Vi,)F2<DX,i) 


i jEtaoS(2,H2(*;>J))F2(D;>J) 


and 


(6.1.7) 


§MB2  = m if13in(27rH2(<,>X,i))F2(DX>i) 


n j=Vin(2-H2(<J))F2(D;fJ). 


1 , K-1/2, 


We  will  show  that  (-  + -)  ''^(C  , S ) and  (—  + q 1 

m n'  MB2  MB2;  Vm  n;  'SlB2*  plB2 ' 

have  the  same  asymptotic  null  distribution.  To  do  this  we  will  need 
to  show  that 


lHN(<J>X,i) 

- V*X,1> 

I - opO), 

(6.1.8) 

IVVj) 

H 

I opO), 

(6.1.9) 

lFN(DX,i) 

- F2(DX.1>I 

H 

■°  V"- 

(6.1.10) 

and 


iF»(DY,j>  - F2(DY,J>I  ‘ %<’>• 

d 


(6.1.11) 


Since  under  Hq,  X.  = for  all  i and  j,  6.1.8  is  equivalent  to  6.1.9 
and  6.1.10  is  equivalent  to  6.1.11.  Hence  we  will  verify  explicitly 
6.1.8  and  6.1.10.  Using  Lemma  3*1.^,  it  will  suffice  to  establish 


Ill 


N 


-1/2 


(Z-y)  = 0(1), 


(6.1.12) 


(6.1.13) 


H 

o 


weg 


2 


(6.1.14) 


and 


H 


o 


(6.1 .15) 


* 

where  <px  ^ ( A)  is  the  angle  between  Xj-A  and  the  positive  half  of  the 
horizontal  axis,  is  a bounded  sphere  in  Ft2,  and  b2  is  a class  of 


2x2  symmetric  matrices  with  uniformly  bounded  elements.  (Conditions 

6.1.12  and  6.1.15  establish  conditions  a and  b,  respectively,  of  Lemma 
3.1.4,  showing  that  6.1.8  holds.  Similarly  conditions  6.1.12  and 

6.1.13  jointly  verify  condition  (a)  of  Lemma  3.1.4  and  6.1.14  verifies 
condition  (b)  of  the  lemma,  establishing  6.1.10.)  Conditions  6.1.12- 
6.1.15  are  established  in  the  following  lemma. 

Lemma  6.1.16.  Let  X^,...,)^,  Y^,...,Y  be  i.i.d.  as  Z,  where 
Z = ^Z1’  Z2^  iS  from  a bivariate  continuous  distribution  function 
with  finite  mean  vector  y and  positive  definite  covariance  matrix 
I^.  Let  Z,  Z,  and  <f>^  ^(A)  be  as  defined  previously.  Assuming  the 
covariance  matrix  of  (Z^  Z2>  Z2,  Z^,  Z^)  is  positive  definite, 
conditions  6.1.12  through  6.1.15  hold. 
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Proof . 

Condition 

6. 1 . 

theorem.  We 

use  Lemma 

3.1. 

/ 

s, . 

s 

1 1 

12 

E = 

S12 

S22 

Z2  “ 


12 


12  ' 


22 


(6.1 


? ’ (S11  ’ 3 22'  312)''  and  ? ’ (4,,'  V V'- 


1 /2 

We  first  verify  that  N (S-A)  is  asymptotically  tri variate  normal 
Let 


V - (Vr  v2,  V v4,  V )' 


_■>  m n m n m n 

N (II  .n  n . I I , z x,ix2/  1 Y11Y21' 

1-1  J-1  J 1-1  j-1  J 1-1  1 41  j-1  U 2J 


where 


and 


m 2 " 2 m 2 n 2 

EX1i+ZY1i’  ZX2i+EY2ir’ 
i-1  j-1  i-1  21  j-1  2j 


Xi  = (X1i’  X2ir  1 = 


Ij  = <V  V' 


j 1 n , 


m+n, 


.17) 


By  the  Multivariate  Central  Limit  Theorem  (recalling  X , 

_1  _m 

Y1’*’*’^n  are  here)  ^(  V-E[V])  is  asymptotically 

multivariate  normal  with  mean  0 and  covariance  matrix  = Var[V], 
Define 


and  let 


g(V)  = (gl(V),  g2(V),  g (V))' 


2 2 

= (V, -V  . V -V  . V -V  V ) ' 
4 r 5 2*  3 1 2‘ 


(S11  ’ 9 22  ’ 3^2y’ 


^y  = y P^y)  ~ E[  V]  , 


then  the  partial  derivatives  of  g.^  evaluated  at  are 


9g1 

avjv-uy  = _2Vi  I v=yv  = ~2uir 

^ 3g1  3g1 

3v^v-pv  = 9v^U=uv  = 9v^ly-uv  = °’ 

and 


Since  these  partials  are  continuous  and  not  all  zero  at  y^,  g has 
nonzero  differential  at  y . (See  remark  3.1.11.)  Similarly 


1 1 it 


and 


!?£,  _ % ff2. 

3V1  _=^V  8V3T^V  ~ ^it  ^v  ’ 

9g2  a«2 

9V^U=yv  = “2p2V’  T5vHv-p  = 1 ; 


9g3 

1 

Sg3, 

!=^v = 

9V,. 1 V> 
o _ 

= o, 
"Hv 

3g3, 

3g3, 

3V1  1 

Y-Uv  = 

~W2V’ 

av2  !=1!v  “ 

!!3, 

3V3 1 

1 • 

So,  g2  and  g3  each  have  nonzero  differentials  at  y Hence,  by  Lemma 
1 /2 

3.1.10,  N (S-A)  is  asymptotically  trivariate  normal.  The  proof  is 
now  easily  completed  by  referring  to  the  proof  of  Lemma  3.1.12,  in 

which  we  showed  that  /^(E_1-E~1 ) is  0 (1)  if  the  elements 

1/2“  ^ 
of  N (E-Eq)  are  jointly  asymptotically  multivariate  normal.  (Given 

the  asymptotic  multivariate  normality  of  the  elements  of  n^2(E-E  )f 

the  exact  structures  of  I and  I are  unimportant,  provided  z is 

o o 

positive  definite.) 

In  verifying  condition  6.1.14,  let 


M = (u1 , v2)  ',  l"1 
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Choose  M such  that  |r.  | < M,  i = 1 , 2 and  |w  | < M,  1 < i < j < 2,  for 
r e and  w e B2*  Then 

sup  | { (X .-(y+N  1/2r))'(E"1+N_1/2w) (X  -(u+N_1/2r))  - (X.-y) 'l"1 (X. -y) } I 
reB1  ~ ~ ~ -l  ~ - ~i  - 2 -1  - 

weB2 

= sup  |{(X.-y)'E21(X.-y)-(X1-y)'Z21(N'1/2r)  - (N_ 1 /2r ' Jz"1 (X^ (y+N_1 /2r ) ) 
re  B ^ 

weB2 


+ (X±- (p+N  1/2r))'N  1/2w(X.-(y+N  1/2r))  - (X . -y ) 'E^  ^ -y ) } | 

= sup  |{-(X  -y)'z"1(N_1/2r)-(N-1/2r')z"1(X.-(y+N"1/2r)) 
reB1  * ~ “ * -1  ~ 

weB2 

+ (X1-  ( y+N~ 1 /2r  ) ) V 1 /2w  (X^  ( y+N_  1 /2r ) ) } | 


S sup  | (X  -y)'z“1(N_1/2r)|  + sup  | (N_ 1 /2r ' )l' (X  -(y+N"1 /2r ) ) I 
reB1  -1  ~ ^ - re6i  ~ 2 ~i  ~ 

we&2  weB2 


+ sup  | (X  -(y+N  1/2r)) 'N~1/2w(X.-(y+N~1/2r)) | . 
reB1  ~ ~ -1  ~ 

weB2 


Considering  the  first  term,  write 
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sup  | (X  -y)'z"1(N_1/2r)  | 
re61  ~ ~ 

wcB- 

- 1 /2  # # # ^ 

■ 3UP  H K«1-u1)l11r1*(X2-u2)4l2r,.(x)-u1)412V(X2-u2M22r2}| 

weB2 

“ 1 /2  # # * « 
s N MI|(Xru1)411|*|(X2-u2)412|.|(X,-ut)412|.|(X2-u2)422|} 

P 

+ 0,  as  N ♦ », 


Similar  arguments  show  that  the  other  terms  also  converge  in 
probability  to  zero,  as  N -*•  «,  establishing  6.1.14. 

As  for  condition  6.1.15,  define 

* -i  /? 

g(h(r))  = <j>x  1 (y+N  r), 

taking 

h(r)  = (X1.-u1-N"1/2r1,  X21-y2-N"1/2r2) 

and  g (h (r ) ) equal  to  the  angle  between  h(r)  and  the  positive  half  of 
the  horizontal  axis.  Now 


sup 

reg 


1 


|h(r)-h(0) | 


i , »_1/2  „-1/2  , | 

= sup  | (-N  r , -N  r )| 
re61 


-1/2,  2 2,1/2 

sup  N (r  +r  ) 

re61 


-1/2  2 1/2 
£ N '^(2M  ) d 


0,  as  N -*■  ®. 
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That  is, 


a.c . 

sup  |h (r  )-h(  0)  | -*■  0. 

re61 

Since  g is  a continuous  function  of  h a.c.,  g o h is  uniformly 
continuous  a.c.  for  re^.  We  have,  then,  that 

a.c. 

sup  | g (h (r ) )-g(h( 0) ) | + 0, 

reB 

as  N “.  This  verifies  6.1.15  and  Lemma  6.1.16  is  established. 

We  now  have  the  preliminary  results  needed  to  establish  the 

asymptotic  approximation. 

Lemma  6.1.18.  Let  X , . . . ,X  and  Y , . . . ,Y  be  independent  random 

~ 1 ~m  _ 1 _n 

samples  from  continuous  distribution  functions  G^u^  u ) and 

°2^U1  ' U2^  = G1^U1~01’  U2~Q2)  ’ resPectively-  If»  under  Hq:  0 = 0 the 

conditions  of  Lemma  6.1.16  hold  and  lira  ^ - X,  0 < X < 1,  where 
N = m+n,  then 


(“  + “)' 
m n 


1 /2 


(CMB2~CMB2’  SMB2~SMB2} 


(0,0)  , 


as  N -*■  ®. 


Proof.  We  first  consider 
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^ ^(m  + n)  (CMB2~EMB2^  ■* 
o 


1 1.-1. 


m 


-1;*;1  - _iicos(2,hnuYiJ)>fn<dyj): 

-[^Ii00s(2,H2(^il))F2(Dx,>1)-iji003(2,H2(»Y*ij)F2(D;j)]}2] 
l ll  1^1 

’ (5*  K>‘  EHoCI51£/CM(2l'H»<*X,l)>FN(DX,i)-°03<2’H2(*X,i))F2(DX,i)] 


* ,,,2. 


-jIi[cos(2,Hn(*IiJ)Fn(DYiJ)-oos(2,H2(*YiJ))F2(DIjJ):)^ 


111  1 m 1 11 
« (1  + [{1  Z U - - E VJ2], 

m n mi=i  i nj  = 1 j 


(6.1.19) 


where 


U1  ‘ °°S(2,'HN(*X,l))FN(I’x,l)-003(2'lH2(+X,l))F2<DX,i)' 


and 


vj  ” cos( 2irHN( 4>y , j ) ) fn (dy, j )~c°s ( 2^H2( 4>Y ^ j ) ) f2(e>y ^ j ) , 

for  i = 1 , . . . ,ra  and  j = 1 n.  Expanding  the  expression  in  6.1.19, 

we  can  write 


EH  ^(m  + n)  (CMB2-EMB2^  ^ 
o 


1 . 1.-1 


m 


■ 


m n 


)}] 


(»  * X U2*  z u u „)*n'2(  z v2*  z v.vj-  — z r u.v. 

m " % i-i  1 i.k  1 k j.i  J j.k  J k "Vi  >i  1 J 


m n 


}] 


1 1.-1, _ ,12  m-1 


1 2 n-1, 


(-  + ~)  {E:  [-U.+  U.UJ  + E C-V  + V.V  ] -E  [2U  V ]},  i * 

m n 1Q  m l m l k Hq  n j n i k Hq  i j 
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Since  under  Hq  U1 U^,  are  exchangeable,  we  can  simplify 

the  above  expression  to  find  that 


^ [{m  + n)  (CMB2~Sb2^  ^ 
o 

■ <;*  b'\  1 * k 

0 

Er  [U^J-Er  CUiUk],  i * k. 
o o 

Since  |U  | £ 2 a.c.,  for  all  i,  we  can  apply  the  Lebesgue  Dominated 

Convergence  Theorem  to  show  that  E [{1  + -)_1/^2(C  -C  )}2] 

H m n MB2  MB2 

o 

converges  to  zero  as  N ■»  °°.  To  do  so,  we  need  to  show  that 
P P 

U.  -*•  0.  Referring  to  6.1.19,  it  is  clear  that  -*•  0 if 


H 


H 


H 


IV*X,i)'H2(*X,i>l  ' V': 


and 


FN(DX,i)"F2(DX,i)l  ' °p(1)> 


As  we  mentioned  these  conditions  (see  6.1.8  and  6.1.10)  are 

established  via  6.1.12-6.1.15  which  we  verified  in  Lemma  6.1.16.  By 

the  Markov  inequality,  it  follows  that 

H 

o 

= Op(1). 


(m  + n} ' 
m n 


1 /2 


(CMB2_CMB2) 


Similarly , 


(m  + n)  (SMB2~^MB2)  * °p(1)* 


Lemma  6.1.18  readily  yields  the  asymptotic  null  distribution  of 
MB2  as  given  in  the  following  theorem. 
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Theorem  6.1.20.  If  the  conditions  of  Lemma  6.1.18  hold  and  the 
distribution  functions  and  are  elliptically  symmetric  (see 
3.1.14),  then  MB2  is  asymptotically  Chi-square  with  two  degrees  of 
freedom  under  HQ. 

Proof.  Using  Lemma  6.1.18,  we  need  only  find  the  asymptotic  null 
distribution  of  the  statistic  6(-  + -)-1(C2  + S2  ).  Let 

m ki  MP  O MD  O ' 


T 


N 


M6(l  + ^)~V/2C  U (6(1  + 

• m n MB2  2 m 


n ' bMB2’ 


be  an  arbitrary  linear  combination  of  (6(-  + -)-1)1/2C  and 

m n MB2 

1 1-11 /2- 

(6(~  + n^  ^ SMB2*  We  wil1  show  that  is  asymptotically  normal 
under  HQ . Write 


'MB2 


m 


m 

Z 5 

i=1 


X , i 


1 n - 
" j = 1 CYJ 


and 


MB2 


1 m - 

“ I Sy  , 
m ^ X , i 


n 


- Z S 
nJ-1  Y'J 


where,  for  all  i and  j, 


(6.1.21) 


X.i 


cos(2,H2(^>1))F2(Dx*>i), 


cos(2tiH2(4,*jJ))F2(D*  ), 


3X,i  ~ sin<  2^h2<  4>x  , i ) > F’2(Dx , i 5 » 


and 
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3Y,j  = sin(2»H2(*y^))F2(Dy^), 


Then  we  can  write 


TN  ” (6(m  + n)  ) *m ^ (*1EX,i+i2EX,i)-  nJf1 


Under  H , ( 4>  , D ) = (<f>  , D ) for  all  i and  j,  and  <j>  . is 

O A,1A,1  I > J t J ''■>1 

* 

independent  of  since  and  are  elliptically  symmetric.  Also 

* * 

H2^X  i^  and  F2^DX  i^  are  eaCh  unif’ormlY  distribution  on  (0,  1),  under 
Hq . Therefore, 


Eu  [Cv  J 


X,i]  ‘ ^ I-CY,j]  ^ [cos(2lTH2(<,>X,i))F2(DXfi)] 


Eh  [cos(2ttH2(<()x  ))]Eh  [F2(Dx  j)] 

o o 

1 1 

= J cos(2iru)du  j udu 
0 0 


= 0 


and 


i 1 

EH  ^EX  i-'  = ^ = J sin(2iru)du  j udu 


(6.1.22) 


Also , 


-2 


VarH  [CX,i-1  = VarH  ^CY,J]  ^ [CX,i-1 


,2  . * 


Eh  [cos  (2,H2(^x  ))]Eh  CF2(Dx  )] 

o o 

r 2 } 2 

I cos  (2iru)du  J u du 

0 0 


1 /6, 
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Var H ['SX,i]  ' VarH„['SY,j’  ' ^ ["SX.il 


1 


J sin2(2iru)du  j u2du 


0 

1 /6, 


and 


C°VHo[CX,i*  SX,i]  = C°VHoCCY,j’  SY,j]  = ^ [CX,iSX,i] 


(6.1.23) 


Eh  Ccos(2TrH2(<J,x>i))sin(2TTH2(^X(.))]EH  [F2(DX  i>] 


= j cos(2im)sin(2iru)du  j u2du 
0 0 


Using  two-sample  U-statistic  theory  (see,  e.g.,  Randles  and 

Wolfe,  1979,  pp.  92-93),  under  Hn  T is  asymptotically  normal  with 

2 2 

mean  0 and  variance  i^  + 8,  . It  follows  that,  under  HQ, 

^6^m  + n^  1^/2^MB2’  ^MB2^  is  cally  bivariate  normal  with 

mean  vector  0 and  covariance  matrix  I2*  Therefore,  by  Lemma  6.1.18, 

the  asymptotic  null  distribution  of  (6(-  + -)”V/2(C  . S ) is 

m n'  ' v MB2  MB2; 

also  bivariate  normal  with  mean  vector  0 and  covariance  matrix  I . 

2 

Hence  the  asymptotic  null  distribution  of 


is  Chi-square  with  two  degrees  of  freedom. 


6. 2 Asymptotic  Distribution  of  MB 2 Under  Contiguous  Alternatives 
The  asymptotic  distribution  of  MB2  under  certain  contiguous 
alternatives  is  easily  obtained  by  extending  the  results  for  the 
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one-sample  statistic,  MB.  We  assume  we  are  sampling  from  the  class  of 
densities  given  by  3.2.1,  namely, 


-[(yen'1  (y-0)]v 


f0(y ) = c exp {- 


} , for  y e R , 


(6.2.1) 


where  c = c(v),  and  h = h(v),  are  given  by  3.2.2,  and  v > 0.  Again, 
since  MB 2 is  invariant  with  respect  to  nonsingular  linear  transforma- 
tions we  may  take  = I2*  Since  MB2  is  in  addition  invariant  with 
respect  to  translations  of  the  observations,  we  may,  without  loss  in 
generality,  take  the  mean  of  the  first  population  to  be  (0,  0)  and  the 
mean  of  the  second  population  to  be  of  the  form  (0,  A).  (See  Section 
4.1.)  Le  Cam's  third  lemma  (Lemma  3.2.4)  will  be  used  to  establish 
the  following  theorem. 

Theorem  6.2.2.  Suppose  X , . . . ,X  and  Y,,...,Y  are  independent 

~ l ~m  _ 1 _n 

random  samples  from  density  functions  f and  f , respectively,  where 

: ?N 

these  densities  are  of  the  form  given  by  6.2.1  with  l2  = I2>  N = m+n, 
and  0^  represents  a sequence  of  alternatives: 

Hn:  0n  = (0,  KN~1/2)\  (6.2.3) 

If  — ■*  A as  N -+  °°  with  0 < A < 1,  then  under  this  sequence  of 
alternatives  the  asymptotic  distribution  of  MB2  is  noncentral  Chi- 
square  with  two  degrees  of  freedom  and  noncentrality  parameter 


,1  1 .-12  , . 

A + 1 -A  SMB(v)‘ 


1 21) 


2 

Here , 6 (v),  given  by  3.2.30,  is  the  noncentrality  parameter  of  the 

asymptotic  distribution  of  the  one-sample  statistic  MB  under  similar 
contiguous  alternatives. 

Proof.  Let  p and  q denote  the  joint  density  functions  of 
N N 

’ ' " ’ jn,Yr  " ”Yn  Under  Ho:  0 = 0 and  hn:  0 = 0N  = (o,  KN~1/2)  , 

respectively.  That  is, 


m 


PM  (X.  , . . • , X , Y , 
N _i  _m 


..*  ) = * W 11  fn(Y-) 
-n  , . 0 _i  0 _j 


i-1 


and 


(6.2.1)) 


m 


3N  ~ qN(^^’‘•,^m,  II !n)  " f0(^i)  ^ *0^^ 


Denote  the  corresponding  probability  functions  by  P and  Q . Let 

N N 

L = ^ 

N PN  ’ 

and 


logL^  = log 


" 9n  (V 
,1-1  -J 


n f (tj. 


Lj-i 


-1  /2 


Letting  A = KN  and  transforming  to  polar  coordinates  via 
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and 


(6.2.5) 


Yy  ' zjSIn<Vj> 


for  all  1 £ j £ n,  where  Z > 0 and  0 £ 4>  < 2ir,  we  can  write 

J * t J 


logL„,  = 


'N  = i MZj2v-(zJ-2Zjsin(^>j)A+A2)V} 


Again  note  that  Z and  41  are  independent  under  P with  respective 
J * >J  N 

densities 


gi(V 


2ircZje 


h 


2v 


0 < Zj  < ®, 


and 


(6.2.6) 


82(*y,J)  " 2it  ' 


0 £ *yiJ  < 2,. 


where  c,  h and  v are  as  previously  defined.  The  Taylor  expansion  used 
in  3.2.11  shows  that 


= — l 2vsin(*  JZ2V  - Z R 
^ /Nh  J-1  Y*j  J Nh  j.i  NJ 


(6.2.7) 


where 


* 5.*  *2.v-2  »2 


Nj 


= 2v(v-1)(Zj-2ZJsin(*yj)Aj  + Aj  ) (Aj  -Z.sin(4>YJ)) 


+ v(zJ2-2zJ.i„(#;iJ,A;+A;2,-ii 
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* 

for  some  A between  0 and  A,  J - 1 n.  In  the  proof  of  Lemma  3.2.8 

we  show  that 


nh 


n 

Z 

i = 1 


,n  k2w 


m 


where  W and  R are  given  in  the  statement  of  that  lemma.  Noting  that 


£ 

(<p't  Z{)  of  Lemma  3.2.8  is  equal  in  distribution  to  ( ^ , z ) under 
P^,  for  all  i and  j,  here  it  can  be  seen  that 


id  J r / k!w 

nh  .L.  n j h ’ 
j = 1 


and  hence 


K_ 

Nh 


n 

Z 

J-1 


Nj 


n Kf  " 
N nh  ^ 


Nj 


( 1 -x) 


k2w 


Therefore 


logL  - -r-  Z 2vsin($*  1)Z2v"1- 

N /Nh  j-1  Y'J  J h 

+ o (1).  (6.2.8) 

N 

Using  the  results  of  Lemma  3.2.23,  it  follows  that,  under  P , 

N 

is  asymptotically  normal  with  mean  (1-A)p  and  variance 
2 - 1 

-2(1-A)p,  where  p = -K  Wh  . So  Q is  contiguous  to  P„. 

N N 
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If  we  verify  that 


tn-  1i(6(S^)"1>,/2WV6(1*s)‘1),/2s! 


MB2 


(6.2.9) 


and  logL  are  jointly  asymptotically  normal  under  p , for 
” N 

(^1>  ^2^  * ^<3,  0) * then  it  will  follow  by  Le  Cam's  third  lemma 

that  Tn  is  asymptotically  normal  under  and  hence 
1 1-11/2 

+ n^  ^ ^CMB2’  SMB2^  wil1  be  ^y^Ptotically  bivariate  normal 
under  Q . We,  therefore,  show  that 


SN  ■ ((6(S  * b~')U2cm2-  <<  * S>',)1/2sHB2-  1o«ln>  (6-2-'0) 


is  asymptotically  multivariate  normal  under  P„.  Note  that 

N 


,N  N.  ,1  1 

+ n>  * (T  + W}’  as  N 


and 


(6.2.11) 


l°gLv 


— I 2vKh“1sin(<t>*  )Z‘7v_1-(  1-A)K2Wh"1+oD  (1) 
/N  j = 1 I,J  3 PN 


(J-  + /N  n J-1 


-(1-A)K2Wh_1+o  (1) 

PN 


/N  r/1-A.1/2  1 n -1  . * 2v-1 

1 T72  t(— : ) - l 2vKh  sin( ())v  '] 


* » j j 


■( 1-A)K2Wh_1+o  (1) 

PN 
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Using  6.2.11  and  Lemma  6.1.18,  the  asymptotic  distribution  under  P 

N 

* 

of  S..  is  the  same  as  that  of 
N 


S„  = 


N 


1/2 


N (1  + _L)1/2 

VA  1 -A ; 


^ ^MB2  ’ ^ ^MB2  ’ ^ 


- (0,  0,  (1-A)K2Wh_1), 


(6.2.12) 


where 


1 n 

— E L , 
nj  = 1 J 


-1  * 2v-1 

Lj  = 2vKh  sin(*y  )Zj  , 


j 1 , • • • ,n y 


and  CMB2  311(1  SMB2  are  given  in  6-1-7.  Consider  then  the  linear 
combination 


.1/2 


N (1  + J_)1/2aiv^  ^MB2+*'2V^  ^MB2+£3  ^ ^ 


V..  = 


'A  1 -A ' 


N1/2  1 m 

7 ; — -rio  1 U./6  CY  .+1./6  Sv  .) 

1 + 1 ) 1 /2  m 1 X,i  2 X,i 


A 1-A 


(6.2.13) 


where  ^ , sx  i » Cy  j » and  j are  glven  in  6.1.21.  Let 
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°C  = VarP  ^ C i ^ = Varp  C J, 

L PN  X>i  PN  Y,J 

°s  - VarPN[/5  §X,i3  • Va,'pN[/5 
°CS  = C0VPN[CX,i’  SX,i]  = COVPNCCY,j’  SY,j]’ 


°CL  = C°VP  CCY,j’  Lj]’  °SL  “ C°VP  [SY , j ’ Lj]’ 


'?  ' VarPNU1A  EX,l*l2A  ^X.i3 

o = Var  [£  /6  C +lJZ  Sy  -1  L.] 

^ 1 ^ ^>J  3 /T  J 


o.  = Var  [L  ]. 
L PN  J 


It  follows,  using  6.1.23,  that 


2 „ 2 2 
°1  ai+a2 


and 


and 

Also,  by  6.1.22  and  3.2.26, 


✓w 

/I  ^ * 3 A ”SL' 


EPnU1/?  h,i*l2ri  \,i3  - 0 


and 


W'l*  J*.3*‘i/55r.j-,3^Lj]  -°- 

Hence  the  asymptotic  distribution  under  P of 

N 


(6.2.14) 


(6.2.15) 


(6.2.16) 
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V.,  = 


_N 

,,,  , , 2 2,1/2 

( ( 1-A ) o1 +Ao2) 


(6.2.17) 


is  standard  normal,  as  seen  by  applying  two-sample  U-statistic  theory 
to  the  representation  of  given  in  6.2.13.  (See,  e.g.,  Randles  and 
Wolfe,  1979,  pp.  92-93.)  Let 

G(t ) = P [Z  < t] 

N J 


and  note  (from  6.1.1)  that  under  P%t 

N 

dy,j  = “ V J = 1 n 

and 

F (t2)  = P (D*  S t2)  = Pp  ( Z £ t)  = G(t ) , 
N ,J  PN  J 


(6.2.18) 


for  t > 0. 


Using  the  independence  of 


* 

4> v . and  Z.  under  P..,  it  follows 

* 9 J J N 


that 

°CL  " C0VPN[5Y,j*  Lj] 

- 2vKh"1Ep^[cos(21rH2(4>*>j)F2(D*>j)sin(4>*  j)Z2v"1] 

= 2vKh_1E  [cos(2ttH  (/  ,))sin(<j>*  ,)]E_  [G(Z  . ) Z?v_1  ] . 
N c J J 

Since,  using  6.2.6, 

H2(t)  - pPN(*;tJ  * t)  - ^ , 


(6.2.19) 


(6.2.20) 


we  find 
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Ep  [cos ( 2irH?( <p  ))sin(<(>  )] 

r^J  C 1 9 J 1 t J 


Ep  [cos(<J>Y  j )sin(<t»Y  j)] 


r 


j cos  ( cj>)sin(  cf>)d4> 


0, 


(6.2.21) 


implying  o = 0.  So 

L/Li 


(l-X)cyAOp 


(l-\)U^lp^\(Z^!Lp^(l-X)Z^-2Z2l3/6/\(T-X)aSL 


,2  .2  .2 


l1+l2+A3(1-X)oL-2H28,3/6/X(1-X)  osl> 


and 


(6.2.22) 


N 


N 


1 /2 


2 , 2.1/2  J 1 ,1/2 


^MB2+l2^  ®MB2+13  yj  ^ 


((1-X)o![  + Xo2)  (j  + YTjj-) ((Jl^  + J!,2+il3(1“X)oL~2*'2X'3’/^Y/X(1“A)aSL) 


4 


is  asymptotically  standard  normal  under  P„.  Hence,  under  P,T 

N N 


N 


I /2 


(i+  _L)1/2 

v A 1-X 


^ SiB2 ’ ^ ^MB2  * yj  *") 


is  asymptotically  multivariate  normal  with  mean  vector  0 and 
covariance  matrix 


1 


0 


0 


-/6  /A ( 1 -A ) o 


SL 


-/F  /A ( 1 -A ) o 


SL 


( 1 — A ) of 


The  asymptotic  distribution  under  P of 

N 
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s»  ■ "<♦  < * y>U\B2.  loeLN> 


is  therefore  also  multivariate  normal,  with  mean  vector 
2 -1 

(0,  0,  — ( 1 -A ) K Wh  ) and  covariance  matrix  r^.  Thus, 


is  under  asymptotically  bivariate  normal  with  mean  vector 
2 -1 

(0,  — ( 1 — A ) K Wh  ) and  covariance  matrix 


2 2 
Z +Z 
1 2 


-%2JZ  ✓( 1 -A) A ogL 


-i2sg  /( i-m  oSL 


( 1 -X) of 


Via  Le  Cam's  third  lemma,  we  conclude  that,  under  Q , the  asymptotic 

N 

distribution  of 


o f ^ K-K1/2n  „ ,,,1  1 .-1.1/2- 

*-1t6(iT  * n>  > > SMB2 


^ ^ o o 

is  normal  with  mean^l^/E  /( 1 — A ) X and  variance  Z^+Z^.  Hence,  under 


,,, 1 1 . “I . 1 /2._ 

6 m + n ) (CMB2’  SMB2) 


is  asymptotically  bivariate  normal  with  mean  vector  (0,  -/6  /( 1 -X) A oc  ) 

oLi 

and  covariance  matrix  I„.  Finally,  it  follows  that  under  Q.t 

£ N 


«B2  - 6(1. 


is  asymptotically  Chi-square  with  two  degrees  of  freedom  and 
noncentrality  parameter, 
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6MB2^  = 6(1-A)Xa|L-  (6.2.23) 

Referring  to  6.2.18,  we  find 

°SL  " C0VPN[§Y,j*  Lj] 

= EPN[SY,j’  Lj3 

= 2vKh  V [sin(2irH  ($*  ,))sin(cf>*  ,)]E_  [G(  Z . ) Z2v" 1 ] . 

N 1 »J  kn  J J 

Since,  using  6.2.6  and  6.2.20, 

E [sin(2irH  (♦*  ) )sin( 4,*  )]  = Ep  [sin2(/  )] 

1 2?  2 * * 

“ pZ  J sin  (41  )d4 

0 

- 1/2,  (6.2.211) 


it  follows  via  3.2.39  and  3.2.40  that  0 = CovD  [S. , L.],  as  given  in 

SL  PN  i 1 

3.2.42.  Hence 


6MB2^v3  = A ( 1 -A ) 6mb ( v ) , 


where  6Mg(v)  is  given  in  3.2.43. 


□ 


6. 3 Asymptotic  Distributions  of  Mardia's  Statistic 
and  a Modification  ~ 

Mardia  (1967)  shows  heuristically  that  his  statistic,  as  defined 
in  5.1.3  and  5.1.4,  is  under  the  null  hypothesis  asymptotically  Chi- 
square  with  two  degrees  of  freedom.  A heuristic  derivation  of  the 
asymptotic  distribution  under  alternatives  is  also  given  by  Mardia. 

In  this  section,  we  will  find  the  asymptotic  distributions,  under  the 


null  hypothesis  and  under  contiguous  alternatives,  of  the  modification 
of  Mardia's  statistic  defined  in  5.1.8.  With  HN  as  defined  in  6.1, 
write 


1 m n 

CM2  = m cos(2irHNUx>i))  " “ ^ cos ( 2uHn( ♦y , j ) ) 


(6.3.1) 


m 


M2 


m ^ sin^2TrHN^Xfi))  ~ n ^ sin(2TrHN(<J>Y^)), 


and  recall  that  the  modification  of  Mardia's  statistic  that  we  will 
use  is 


M2  = 


2(~  + V 1 

m n 


Using  the  definitions  of  <p  , $ , and  H„  as  given  in  section  1 of 

^ 9 * I 9 J <- 

this  chapter,  we  introduce  a statistic  asymptotically  equivalent  to  M2 
under  the  null  hypothesis.  Let 


M2 


1 m # 

in  if1  cos(2TrH2(^x>i)) 


n ^ cos(2wH2(^iJ)) 


and 


(6.3.2) 


m 


§M2  ' I ^ s‘n(2"H2<*x*,1>) 


1 11  * 

- I sin(2iTH  ( 4>  )) 

n j _ i *-  1 « J 


The  following  lemma  shows  that  (-  + -)""^2(C  . S ) and 

m n'  v M2’  M2; 

(^  + ^)  have  the  same  s^yroptotic  null  distribution. 
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Lemma  6 . 3 . 3 - Let  X,  ,...,X  and  Y ,Y  be  independent  random 

-I  _m  _1  _n 

samples  from  bivariate  populations  with  continuous  distribution 

functions  G^ (u^  , u^)  and  G^u^  , u^)  = G^Cu^-G^,  U2_02^’  resPect* velY* 

If  under  Hq:  0 = (0^  ©2) ' = 0,  ECX^  = E[Y  ] = y is  finite, 

Var[X.]  = Var[Y  ] * Z is  positive  definite,  and  lim  ^ = X, 

_ l 12  N 

N+“ 

0 < X < 1 , then 


(m  + n}  (CM2”Si2’  SM2~^M2) 


P 

- (0,  0), 
H 

o 


as  N -*•  ®, 

Proof.  Referring  to  the  proof  of  Lemma  6.1.18,  we  see  that  we 
can  write 


EH  [{(S  4 Hr1/2(CH2-5M2))2]  • EH  [W’  1 • k 

O 0 0 


where 


(6.3.10 


Ui  *=  cos  (2irH^(  (Ji^  ^ )-cos  (2ttH2(<|)x  i)). 


for  i = 1,...,m.  An  examination  of  the  proof  of  Lemma  6.1.16  shows 
that  the  conditions  stated  here  are  sufficient  for  6.1.12  and 
6.1.15.  These  two  conditions  were  shown  to  imply,  via  Lemma  3 • 1 . ^ , 


that 


iVViWVi5 


Op(1) 


(6.3.5) 
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7T 

Hence,  IL  = o (1)  and » by  Lebesgue  Dominated  Convergence  Theorem, 

EH  [lW  -«<’>• 

o o 

Therefore , 


(m  + n)"1/2(CM?"^M?)  = ° (1) 

m n M2  M2  p 


and,  similarly, 


i i -i /?  _ u 

(£  + b (sM,-sM,)  - od) 

m n M2  M2  p 


□ 


Using  Lemma  6.3.3  we  establish  the  following  theorem. 

Theorem  6.3.6.  Under  the  conditions  of  Lemma  6.3.3,  the 
asymptotic  null  distribution  of  M2  is  Chi-square  with  two  degrees  of 
freedom . 

Proof.  Let 


tn- <2<5*i)'1),/2<t,EM2*12V 


1 . 1X-1. 1/2,1 


m 


(2(m  n)  ) {m  J 1 (£1  CMX  , i+S,2SMX  ,i) 


" n Jf1U1CMY,j+£2^MY,j)}’ 


where,  for  1 ;£  i £ m and  1 < j < n, 


CMX,i  = cos(21;H2(4>x>i)), 


(6.3.7) 
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CMY,j  = cos  ( 2ttH2(  4>y  )), 


SMX, 1 ‘ 3ln(2»H2(*xl)), 


and 


SMY,j  ain(2irH2(^y^)). 

* * 

Since  and  h2^y  each  uniformly  distributed  on  (0,1) 

under  the  null  hypothesis  it  follows  that,  for  all  i and  j, 


EH  [SlX,i]  “ ^ [CMY,j]  " ^ [EMX,i]  " ^ [EMY,j]  " °’ 


VaPH  [^MX,i]  = VarH [SlY,j]  VarH  [^MX,i]  _ VarH  [EMY,j]  = 2 


and 


(6.3.8) 


C°1 * * * VH  [EMX,i’  EMX , i ^ * C0VH  ^MY.J*  EMY,j]  " °* 


Hence,  by  U-statistic  theory,  is  under  HQ  asymptotically  normal 

2 2 

with  mean  0 and  variance  £ +t  . The  asymptotic  null  distribution  of 

1 1-11/2- 

(2(-  + — ) ) ^2^  is>  thus,  bivariate  normal  with  mean  vector 

0 and  covariance  matrix  1^.  Lemma  6.3.3  shows  that  this  is  also  the 

asymptotic  null  distribution  of  (2(-  + -)_1)1/2(C  S ).  The 

m n M2  M2 

asymptotic  null  distribution  of  M2  follows.  I I 

The  asymptotic  distributions  under  continguous  alternatives  of  M2 
and  B,  Blumen's  one-sample  statistic,  are  related  in  the  same  way  as 
that  of  the  two-sample  statistic,  MB2,  and  the  one-sample  statistic, 
MB.  We  next  obtain  the  asymptotic  distribution  of  M2  under  contiguous 


alternatives . 
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Theorem  6.3.9.  Under  the  conditions  of  Theorem  6.2.2,  the 

asymptotic  distribution  under  H„:  0 = 0.  = (0,  KN~1/2)  of  M2  is 

N _ _N 

noncentral  Chi-square  with  two  degrees  of  freedom  and  noncentrality 
parameter 

r2  . .1  1-1  2,  . 

^M2  V = A + 1 - A 6B 

2 

where  5g(v),  given  by  3.3.7,  is  the  noncentrality  parameter  of 
Blumen's  one-sample  statistic  under  Pitman  alternatives. 

Proof.  We  proceed  exactly  as  in  the  proof  of  Theorem  6.2.2. 
Replace  in  6.2.13  by 


V.,  = 


N1  /2  1 m 

1 , 1 jl/2  {m  ^l^  Sx,i  + l2,/2  ^MX,i) 


‘A  1-A 


n ^MY,j  + *'2’/2  ^MY , j 


1 - A 


LJ} 


A J 


Let  Lj  and  L be  as  defined  in  6.2.12  and  define 


°CM,L  = C°VP^[CMY, j ’ LJ]» 


(6.3.10) 


°SM,L  ‘ C°VPN[SMY,j’  Lj]’ 

.2 


and 


°1  = VarPN[il1,/2  ^MX,i  + l2,/2 


°2  “ VarPNl-ill’/2  ^MY,j  + *'2l/2  ®MY,J-13  ^ Lj]' 
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By  6.3.8  a'2  * 8,^  + 8, ^ and  write 


2 2 2 2 -1  2 
0'  = 

-2^1  f.^/2  — r-  °CMiL 
‘“zV5  ^ °SM,L* 


Appealing  to  6.2.19-6.2.21,  it  is  seen  that  o„w  , = 0.  Following  the 

CM,  L 

approach  given  in  the  proof  of  Theorem  6.2.2  (see  6.2.16  and  6.2.17) 
we  find  that 


,1/2 


l,/2  C^t/2  SH2.l3 


/I -A 

/I 


«i-i>„,'2.ja'2>  (i  * ^P'2  (i2n2n2(,-no^2t2J3/2/rrTyT0sH_L) 


1 /2 


is  asymptotically  standard  normal  under  P . This  implies  that 


,1  /2 


1 1 ,1/2  Sl2’  ^ ^M2’ 


/I -A  T- 


(A  + 1-A} 


/A 


L) 


is,  under  P , asymptotically  trivariate  normal  with  mean  vector  0, 
and  covariance  matrix 


M2 


,1  -/2  /(I -A) A 0 


SM,  L 


-/2  /( 1 -A ) A 0 


SM,  L 


( 1 - A ) of 


/ 


The  asymptotic  distribution  under  P of 

N 


1 40 


(<2<m  * n>~,>1/2cM2’  10«Ln> 


is  therefore,  using  Lemma  6.3.3  and  the  proof  of  Theorem  6.2.2,  also 

2 - 1 

trivariate  normal  with  mean  vector  (0,  0,  - ( 1 — A ) K Wh  ),  and 

covariance  matrix  r . Under  P , 

2 N 


is  asymptotically  bivariate  normal  with  mean  vector, 
2 -1 

(0,  - ( 1 — X ) K Wh  ),  and  covariance  matrix, 


2 2 
9 Si 

1 2 


°SM,L 


-£2/2  /TPTTT  »SHiL 


(1  -»)o 


Le  Cam's  third  lemma  implies  that,  under  the  sequence  of  alternatives 
Qn,  the  asymptotic  distribution  of 


V2(m  + n^^WM^m  + ^"V^S 

i m n M2  2 m n M2 


^ O O 

is  normal  with  mean  -2.-/2  /( 1 -A)  A a . and  variance  27  + 2-.  Therefore 

2 oM  | Li  12 


,-,1  1.-1. 1/2.„  „ . 

(2(m  + ^ ) (C«-  W 

m n M2  M2 


is  asymptotically  bivariate  normal  under  (0  with  mean  vector, 

N 


ni 


(0,  -/2  /( 1 — X ) A a ),  and  covariance  matrix,  I„.  Under  Q . the 
bM, L 2 N 

asymptotic  distribution  of  M2  is  thus  noncentral  Chi-square  with  two 
degrees  of  freedom  and  noncentrality  parameter 


(v) 


2( 


1 — X ) X a 


2 

SM,  L* 


Furthermore 

°SM, L = CovPNC^MY,j’  Lj] 

= 2vKh_1E  [sin(2TrH  ($*  . ))sin($*  .)Z,V_1] 

^ 1 1 »J  J 

■ CovPNc§Bi-  Li3 

6„<v)  1/2 

- (-V">  ■ 

using  3.3.11,  3.3.12,  3. 2. MO,  and  6.2.2M. 

It  can  be  checked  that  the  heuristically  developed  expression  for 
the  noncentrality  parameter  of  Mardia's  statistic,  given  by  12.19  of 

p 

his  1967  paper,  agrees  with  6^0(v)  for  the  contiguous  alternatives 


under  consideration. 


CHAPTER  7 

COMPARISONS  OF  TWO-SAMPLE  TEST  STATISTICS 


7 . 1 Asymptotic  Relative  Efficiencies 

In  this  section  we  develop  asymptotic  comparisons  of  the  proposed 

statistic,  MB2 , with  the  modification  of  Mardia's  statistic,  M2,  and 

2 

with  the  two-sample  Hotelling's  T 2 (see  1.1.11).  Again  due  to  the 
affine-invariance  of  these  statistics,  we  may,  without  loss  of 
generality,  assume  that  the  common  population  covariance  matrix  is  the 
identity  and  we  need  only  consider  alternatives  of  the  form 

V ! = ®n  = (0*  AK- 

Under  , with  A = KN_1/2,  if  jj  -*•  \ , 0 < X < 1 , as  N -*■  « , 

2 

then  the  asymptotic  distribution  of  T^  is  Chi-square  with  two  degrees 
of  freedom  and  noncentrality  parameter 


6 


(7.1.1) 


(See,  e.g.,  Puri  and  Sen,  1971,  p.  212.)  This  result  in  conjunction 
with  Hannan's  work  (see  4.1.2)  and  Theorems  6.2.2  and  6.3.9  verifies 
the  following  result. 

Theorem  7.1.2.  Suppose  X , . . . ,X  and  Y,,...,Y  are  independent 

~ 1 _m  _ 1 _n 

random  samples  from  density  functions  f and  f , respectively,  where 

° ?N 

these  densities  are  of  the  form  given  by  6.2.1  with  » I^  and 
N = m+n.  For  testing 


1 42 


H : 0=0  against  H : 0 = 0 . 

o a _ _N 

if  v > 1/2  and  lim  ^ = X,  0 < X < 1,  then 
N+»  N 


ARE(M2,Tp 


4m 


and 


ARE(MB2 , Tp 


6MB2(v) 


6MB(v) 


That  is,  the  Pitman  asymptotic  relative  efficiencies  of  MB2  and 
2 

M2  relative  to  are  precisely  those  of  the  one-sample  versions,  MB 
and  B,  respectively,  relative  to  the  one-sample  Hotelling’s  T2.  Table 
4.1  gave  one-sample  Pitman  asymptotic  relative  efficiencies  for 
selected  values  of  v.  For  convenience,  we  display  these  efficiencies 
again  in  Table  7. 1 . 


Table  7.1 

Pitman  Asyptotic  Relative  Efficiencies 


V 

ARE(M2,  T§) 

ARE(MB2,  T|) 

10.0 

.50530 

1 .33231 

5.0 

.51872 

1 .221  01 

1.0 

.78540 

.98464 

.8 

.91  839 

.99574 

.7 

1.02988 

1.01488 

. 5+ 

1.50000 

1 .12500 

Corresponding  to  the  results  in  the  one-sample  setting,  MB 2 is  quite 

2 

competitive  with  for  all  the  distributions  considered.  It 

2 

dominates  M2  and  T2  for  light-tailed  distributions.  M2,  on  the  other 
hand,  performs  best  asymptotically  for  heavy-tailed  distributions. 

The  conjectured  efficiencies  from  Table  4.2  for  0 < v < 1/2  are  given 
in  Table  7.2. 


Table  7.2 

Conjectured  Pitman  Asymptotic 
Relative  Efficiencies 


V 

ARE(M2,  T|) 

ARE(MB2,  T|) 

.4 

2.09838 

1.27541 

.3 

3.69648 

1.64385 

.2 

11.59689 

3.00306 

.1 

366.57838 

25.03483 

7.2  Monte  Carlo  Results 

In  this  section  we  compare  several  two-sample  procedures  using 
Monte  Carlo  simulations.1  The  affine  invariant  statistics  M,  T^,  and 
MB2  are  examined  along  with  the  bivariate  Wilcoxon  Rank  Sum  statistic 
proposed  by  Chatter jee  and  Sen  (1964).  As  discussed  in  Chapter  1,  the 
statistic  of  Chatterjee  and  Sen  is  not  affine  invariant.  It  is  of  the 
form: 


1 


All  calculations  and  random  variable  generations  were  performed 
with  an  IBM-4341  running  on  a vm/cms  operating  system.  Fortran  77 
and  several  IMSL  library  subroutines  were  used. 


1 45 


[FT(sr  s2>]'(nv)',[n7T<!V  V]' 

where  and  are  the  coordinate-wise  Wilcoxon  Rank  Sum  statistics, 

as  given  in  1.1.12,  and  V is  a consistent  estimator  of  the  asymptotic 

covariance  matrix  of  (S  , S ).  Using  Hettmansperger  (1984, 

/N(N+1)  1 2 

pp.  289-290),  we  chose  ■)  g( nI( N+1 ) f>or  the  dia8°nal  elements  of  V and 


mn 

N2(N-1) 


m 
{ I 
i = 1 


Rank(X1.)Rank(X2i ) 
( N+1  ) 2 


Rank(Y,  . )Rank(Y„.  ) 

1_i 2i 

(N+1)2 


for  the  off-diagonal  elements. 

2 

The  four  statistics  M,  T , MB2,  and  CS  were  compared  using 
samples  from  the  class  given  in  6.2.1,  Pearson  Type  II  (light-tailed) 
and  Type  VII  (heavy-tailed)  distributions,  and  bivariate  normal 
mixtures.  (These  same  distributions  were  examined  in  comparing  the 
corresponding  one-sample  statistics.  See  Section  4.2  of  Chapter  4 for 
a more  detailed  discussion  of  these  distributions  and  the  generation 
procedure . ) 

2 

The  performances  of  M,  T , and  MB2  parallel  the  performances  of 

2 

their  one-sample  counterparts,  namely,  B,  T , and  MB,  respectively. 

Looking  at  the  empirical  powers  in  Table  7.3  (m  = n = 15, 
significance  level  = .05813  as  in  the  one-sample  case),  we  see  that 
Mardia's  statistic,  M,  has  the  greatest  power  for  the  heaviest-tailed 
distributions  (v  = .10  and  v = .50)  and  Hotelling's  T2  performs  quite 
poorly  in  the  heaviest-tailed  situation.  For  the  bivariate  normal 


Table  7.3 

Two-Sample  Monte  Carlo  Results  with 
m = 15,  n = 15,  and  a = .0531  8 


Bivariate 

Statistic 

median 

difference 

M 

T2 

2 

MB  2 

CS 

(.000,. 000) 

.057 

.040 

v = .10 
(.025) 

.048 

.047 

( .01  4,  .01 11) 

.068 

.044 

(.025) 

.055 

.052 

(.028, .028) 

.100 

.045 

(.026) 

.066 

.075 

( .042, .042) 

.146 

.051 

(.029) 

.086 

.108 

(.000,. 000) 

.052 

.075 

v = .50 
(.049) 

.065 

.051 

( .141  , .141  ) 

.070 

.102 

( .070) 

.083 

.079 

(.283,-283) 

.162 

.191 

(.149) 

.178 

.1  61 

( .424, .424) 

.328 

.355 

(.279) 

.315 

.308 

( .000, .000) 

.059 

.086 

v = 1.00 
(.060) 

.067 

.063 

(.141  , .141) 

.068 

.104 

( .082) 

.087 

. 066 

(.283, .283) 

.118 

.183 

(.143) 

.164 

.133 

( .424, .424) 

.204 

.334 

(.276) 

.293 

.247 

( .000, .000) 

.049 

.078 

v = 50.00 
(.063) 

.073 

.046 

(.141, .141) 

.057 

.105 

( .073) 

.090 

.064 

(.283, .283) 

.090 

.188 

(.146) 

.175 

.129 

( .424, .424) 

.142 

.31  4 

( .250) 

.305 

.21  3 

( .000, .000) 

.050 

.075 

v = 120.00 
(.051) 

.059 

.042 

( .141  , .141  ) 

.058 

.095 

(.074) 

.083 

.071 

(.283, .283) 

.082 

.172 

(.134) 

.175 

.117 

( .424, .424) 

.110 

.289 

( .238) 

.337 

.21  6 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  using  the  F?  N_p 
distribution.  ’ 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 15. 

Distributions:  Class  given  in  6.2.1,  I = I . 

2 2 


distribution  (v  = 1.0)  and  distributions  close  to  the  bivariate 
2 

normal,  MB2  and  perform  similarly.  For  very  light-tailed 
distributions  (v  = 120),  MB2  has  greater  power  than  each  of  the  three 
competitors  included  in  the  study.  The  statistic  of  Chatterjee  and 
Sen  performs  well,  although  again  problems  due  to  its  nonaf fine- 
invariant  nature  can  be  seen  in  Table  7.M.  A dramatic  loss  in  power 
occurs  as  the  correlation  decreases  to  -.99. 

Changes  in  significance  levels  and/or  sample  sizes  (Tables  7.5  - 
7.10)  produce  similar  results.  For  some  of  the  sample  sizes  examined, 
the  exact  null  distribution  of  Mardia's  statistic  was  available.  (See 
Mardia,  1967,  and  Mardia,  1969.)  In  these  cases,  critical  values  were 
included  from  the  exact  and  the  asymptotic  null  distributions. 
Significance  levels  were  selected  from  the  naturally  occurring 
percentiles  and  nonrandomized  tests  were  performed  using  both  critical 
values.  An  inspection  of  the  simulated  powers  under  the  null 
hypothesis  suggests  that  samples  of  size  1 2 or  so  may  be  needed  to  use 
the  asymptotic  null  distribution  of  MB2. 

Simulations  using  the  light-tailed  Pearson  Type  II  and  heavy- 
tailed Pearson  Type  VII  distributions  (Table  7.11)  yield  the  same 
general  pattern  of  results.  The  bivariate  normal  mixtures  studied 
(Table  7.12)  suggest  that  MB2,  like  MB,  is  robust  under  departures 
from  elliptical  symmetry.  (Recall  Mixtures  1,  2,  and  3 are  not 

elliptically  symmetric,  with  Mixture  3 being  bimodal.)  Lastly,  the 

2 

poor  performance  of  Hotelling's  T^  relative  to  MB2  is  demonstrated  for 
heavy-tailed  mixtures  (Mixtures  ^ and  5). 
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Table  7.4 

Two-Sample  Monte  Carlo  Results  with  m = 15,  n = 15, 
a = .0531  8,  and  Varying  Correlation  Structures 


Bivariate 

median 

difference 

Statistic 

M 

T 

•2 

2 

MB  2 

CS 

Bivariate 

! Normal  (unit  variances 

, correlation  = 

0.0) 

(.000, .000) 

.041 

.072 

(.048) 

.049 

.051 

( .141  , .141  ) 

.055 

.105 

(.065) 

.075 

.081 

(.284, .284) 

.117 

.1  83 

(.145) 

.153 

.172 

( .424, .424) 

.21 1 

.325 

( .256) 

.280 

.290 

Bivariate 

Normal  (unit 

variances , 

correlation  = 

-.50) 

( .000, .000) 

.048 

.063 

(.048) 

.057 

.048 

(.10, .10) 

.061 

.094 

(.056) 

.082 

.057 

( .20, .20) 

.094 

.167 

(.128) 

.141 

.111 

(.30, .30) 

.1  94 

• 315 

(.254) 

.265 

.21  3 

Bivariate 

Normal  (unit 

variances , 

correlation  = 

-.99) 

( .000, .000) 

.055 

.082 

(.054) 

.056 

.049 

( .01 4, .01 4) 

.064 

.103 

( .077) 

.080 

.059 

(.028, .028) 

.120 

.184 

(.145) 

.157 

.084 

( .042, .042) 

.203 

.334 

(.261) 

.273 

.132 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  using  the  F2 
distribution. 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 15. 

Distributions:  Bivariate  normal. 


Table  7.5 

Two-Sample  Monte  Carlo  Results  with  m = 15,  n = 15,  and  a = .10 


Bivariate 

Statistic 

median 

difference 

M 

f2 

l2 

MB  2 

CS 

( .000, .000) 

.100 

.099 

v = .10 
(.064) 

.117 

.092 

( .01  if,  .01  if) 

.125 

.101 

( .068) 

.122 

.108 

( .028, .028) 

.172 

.104 

(.071) 

.147 

.144 

( . 042 , . 042) 

.241 

.114 

( .075) 

.180 

.202 

( .000, .000) 

.100 

.133 

v =*  . 50 
( .098) 

.118 

.102 

( .141  , .141  ) 

.152 

.168 

( .138) 

.167 

.147 

(.283, .283) 

.279 

.317 

(.260) 

.294 

.277 

( .424, .424) 

.474 

.496 

( .437) 

.461 

.467 

(.000, .000) 

.100 

.121 

v = 1.00 
(.098) 

.118 

.101 

(.141  , .141) 

.126 

.173 

( .137) 

.154 

.133 

(.283,-283) 

.21 1 

.277 

(.234) 

.242 

.234 

( .424, .424) 

.334 

.451 

( .400) 

.410 

.374 

( .000, .000) 

.103 

.135 

v = 50.00 
(.104) 

.127 

.103 

(.141  , .141) 

.124 

.166 

( .139) 

.161 

.133 

(.283,-283) 

.182 

.280 

(.237) 

.277 

.220 

( .424, .424) 

.256 

.438 

(.389) 

.443 

.353 

(.000, .000) 

.107 

.128 

v = 120.00 
(.095) 

.112 

.099 

( .141  , .141 ) 

.121 

.173 

( .135) 

.168 

.136 

(.283,-283) 

.150 

.270 

(.225) 

.306 

.229 

( .424, .424) 

.201 

.425 

(.367) 

.500 

.352 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  using  the  F2  N_o 
distribution.  ’ 

Nominal  significance  level:  .10. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 15. 

Distributions:  Class  given  in  6.2.1,  l = I . 
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Table  7.6 


Two-Sample  Monte 

Carlo 

Results  with  m = 15,  n = 

1 5,  and  a = 

.01 

Bivariate 

Statistic 

median 

difference 

M 

t2 

2 

MB2 

CS 

( .000, .000) 

.009 

.008 

v = .10 
(.002) 

.009 

.008 

( .01  4, .01 4) 

.010 

.009 

( .003) 

.01 1 

.013 

( .028, .028) 

.01  9 

.009 

(.005) 

.017 

.015 

( .042, .042) 

.034 

.011 

(.007) 

.023 

.027 

( .000, .000) 

.010 

.020 

v = .50 
(.010) 

.004 

.007 

(.141  , .141) 

.01  4 

.029 

( .011  ) 

.01  7 

.015 

(.283, .283) 

.048 

.077 

(.046) 

.045 

.043 

( .424, .424) 

.119 

.172 

(.104) 

.100 

.095 

( .000, .000) 

.01  2 

.025 

v = 1.00 
(.010) 

.01  4 

.010 

( .141  , .141  ) 

.01  6 

.039 

( .022) 

.024 

.01  6 

(.283, .283) 

.033 

.075 

(.043) 

.046 

.033 

( .424, .424) 

.065 

.171 

(.099) 

.094 

.071 

(.000, .000) 

.01  3 

.033 

v = 50.00 
(.015) 

.01  4 

.012 

(.141  , .141  ) 

.01 1 

.048 

( .020) 

.020 

.014 

(.283,-283) 

.025 

.082 

(.047) 

.048 

.032 

( .424, .424) 

.040 

.153 

( .097) 

.103 

.069 

( .000, .000) 

.01  2 

.031 

v = 120.00 
(.012) 

.01  4 

.010 

( .141  , .141 ) 

.01 1 

.040 

( .025) 

.01  8 

.01  7 

(.283,-283) 

.01  4 

.070 

(.037) 

.057 

.028 

( .424, .424) 

.028 

.148 

( .094) 

.135 

.069 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling’s  T?  using  the  F2  N_o 
distribution.  * 

Nominal  significance  level:  .01. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 15. 

Distributions:  Class  given  in  6.2.1,  £ = I . 

2 2 
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Table  7.7 

Two-Sample  Monte  Carlo  Results  with 
m = 15,  n = 10,  and  a = .0531  8 


Bivariate 

Statistic 

median 

difference 

M 

t2 

2 

MB  2 

CS 

( .000, .000) 

.044 

.032 

v = .10 
(.015) 

.058 

.035 

( .01 4, .01 4) 

.059 

.031 

(.019) 

.071 

.040 

(.028, .028) 

.092 

.040 

(.020) 

.086 

.078 

( .042, .042) 

.132 

.053 

( .022) 

.116 

.112 

(-  000,  .000) 

.038 

.061 

v = . 50 
(.036) 

.059 

.033 

(.141, .141) 

.076 

.096 

(.063) 

.083 

.061 

(.283,-283) 

.152 

.190 

(.127) 

.154 

.129 

( .424, .424) 

.289 

.335 

( .254) 

.285 

.271 

(.000, .000) 

.039 

.067 

v = 1.00 
(.039) 

.046 

.029 

( .141  , .141  ) 

.061 

.097 

( .061 ) 

.077 

.063 

(.283, .283) 

.092 

.169 

(.125) 

.1  41 

.108 

( .424, .424) 

.175 

.310 

( .225) 

.235 

.199 

(.000, .000) 

.034 

.078 

<L 

• II 
O 

-tr  t_n 

o 

w • 

o 

o 

.064 

.037 

( .141  , .141  ) 

.052 

.091 

( .067) 

.083 

.052 

(.283,-283) 

.078 

.163 

(.113) 

.136 

.098 

( .424, .424) 

.135 

.285 

(.212) 

.252 

.178 

(.000, .000) 

.037 

.071 

v = 120.00 
(.044) 

.055 

.046 

(.141  , .141) 

.050 

.093 

( .065) 

.080 

.056 

(.283,-283) 

.067 

.170 

(.111) 

.167 

.106 

( .424, .424) 

.101 

.301 

(.214) 

.308 

.192 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  T^  using  the  F2 
distribution. 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 10. 

Distributions:  Class  given  in  6.2.1,  E = I . 


Table  7.8 

Two-Sample  Monte  Carlo  Results  with 
m = 1 2,  n = 1 2,  and  a = .0531 8 
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Bivariate 

Statistic 

median 
dif f erence 

M 

t2 

12 

MB  2 

CS 

( .000, .000) 

.040 

.042 

v = .10 
(.019) 

.051 

.046 

( .01  4, .01 4) 

.055 

.048 

(.021) 

.062 

.053 

( .028, .028) 

.088 

.053 

( .028) 

.074 

.071 

( .042, .042) 

.120 

.063 

(.033) 

.089 

.095 

(.000,. 000) 

.049 

.072 

v = . 50 
(.039) 

.067 

.043 

( .141  , .141  ) 

.081 

.106 

( .067) 

.082 

.067 

(.283, .283) 

.145 

.1  92 

(.145) 

.168 

.147 

( .424, .424) 

.277 

.317 

( .242) 

.275 

.278 

(.000,. 000) 

.058 

.077 

v = 1.00 
(.052) 

.051 

.045 

( .141  , .141 ) 

.059 

.102 

( .065) 

.067 

.063 

(.283, .283) 

.099 

.181 

(.129) 

.150 

.116 

( .424, .424) 

.181 

• 313 

(.243) 

.250 

.227 

( .000, .000) 

.052 

.077 

v = 50.00 
(.051) 

.061 

.039 

( .141  , .141  ) 

.051 

.095 

(.061) 

.078 

.051 

(.283,-283) 

.078 

.169 

(.107) 

.144 

.087 

( .424, .424) 

.119 

.288 

( .212) 

.262 

.182 

(.000,. 000) 

.051 

.071 

v - 120.00 
(.048) 

.057 

.042 

(.141  , .141  ) 

.058 

.089 

(.059) 

.082 

.051 

(.283,-283) 

.067 

.157 

(.103) 

.156 

.096 

( .424, .424) 

.097 

.273 

(.198) 

.288 

.174 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  using  the  F2  N_o 
distribution. 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m - 12,  n = 12. 

Distributions:  Class  given  in  6.2.1,  = 


Table  7.9 

Two-Sample  Monte  Carlo  Results  with 
m = 10,  n = 10,  and  a = .04936 
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Bivariate 

Statistic 

median 
differ ence 

M 

T2 

MB2 

CS 

( .000, .000) 

.051 

(.052) 

.051 

v = .10 
(.023) 

.058 

.037 

( .01  4,  .01  4) 

.057 

( .062) 

.052 

( .023) 

.067 

.046 

( .028, .028) 

.069 

(.078) 

.056 

(.028) 

.075 

.061 

( .042, .042) 

.087 

(.097) 

.062 

( .032) 

.086 

.086 

( .000, .000) 

.053 

(.057) 

.078 

v = .50 
(.042) 

.058 

.042 

(.141  , .141) 

.055 

(.064) 

.117 

( .059) 

.077 

.052 

(.283, .283) 

.105 

(.117) 

.175 

(.117) 

.133 

.109 

( .424, .424) 

.208 

(.216) 

.283 

( .206) 

.241 

.202 

(.000, .000) 

.055 

( .060) 

.093 

v = 1.00 
(.047) 

.064 

.043 

( .141  , .141  ) 

.065 

( .072) 

.104 

(.063) 

.078 

.048 

(.283, .283) 

.086 

(.096) 

.165 

(.100) 

.113 

.084 

( .424, .424) 

.144 

( .156) 

.246 

(.176) 

.179 

.154 

( .000, .000) 

.054 

(.059) 

.108 

v = 50.00 
(.065) 

.075 

.049 

(.141, .141) 

.054 

( .060) 

.114 

(.077) 

.086 

.064 

(.283, .283) 

.075 

( .080) 

.155 

(.097) 

.122 

.087 

( .424, .424) 

.111 

( .119) 

.244 

( .162) 

.198 

.126 

(.000, .000) 

.057 

(.065) 

.106 

v = 120.00 
(.062) 

.067 

.058 

( .141  , .141 ) 

.058 

(.065) 

.116 

(.076) 

.085 

.064 

(.283, .283) 

.069 

(.077) 

.164 

(.105) 

.133 

.091 

( .424, .424) 

.100 

(.103) 

.234 

(.165) 

.21 1 

.149 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  and  Mardia's 
statistic  using  the  F2  distribution  and  the  exact  null 
distribution  of  M,  respectively. 

Nominal  significance  level:  .04936. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m - 10,  n = 10. 

Distributions:  Class  given  in  6.2.1,  E = I . 
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Table  7.10 

Two-Sample  Monte  Carlo  Results  with 
m = 8,  n = 8,  and  a = .04973 


Bivariate 

Statistic 

median 
differ en ce 

M 

t2 

2 

MB2 

CS 

( .000, .000) 

.051 

( .044) 

.052 

v = .10 
(.011) 

.065 

.035 

( .01 4, .01 4) 

.049 

( .040) 

.053 

(.009) 

.065 

.040 

(.028,  .028) 

.072 

(.064) 

.055 

(.013) 

.075 

.047 

( .042, .042) 

.100 

(.090) 

.059 

( .021 ) 

.092 

.067 

( .000, .000) 

.047 

(.043) 

.076 

v = .50 
( .042) 

.063 

.042 

( .141  , .141  ) 

.054 

(.048) 

.108 

(.049) 

.079 

.042 

(.283,-283) 

.092 

(.084) 

.158 

(.086) 

.109 

.084 

( .424, .424) 

.160 

( .147) 

.259 

(.155) 

.195 

.150 

( .000, .000) 

.056 

(.052) 

.109 

v = 1.00 
(.047) 

.075 

.040 

( .1 41 , .141  ) 

.063 

( .058) 

.111 

( .065) 

.083 

.053 

(.283, .283) 

.080 

(.072) 

.155 

(.094) 

.108 

.075 

( .424, .424) 

.123 

(.108) 

.234 

(.146) 

.166 

.125 

( .000, .000) 

.045 

(.041) 

.098 

v = 50.00 
(.057) 

.063 

.046 

(.141, .141) 

.047 

( .042) 

.114 

( .062) 

.069 

.055 

(.283, .283) 

.061 

(.056) 

.148 

(.082) 

.102 

.065 

( .424, .424) 

.095 

(.083) 

.227 

(.133) 

.163 

.100 

(.000, .000) 

.045 

(.039) 

.085 

v = 120.00 
(.052) 

.050 

.036 

( .1  41  , .141  ) 

.050 

( .042) 

.091 

( .059) 

.060 

.044 

(.283,-283) 

.059 

(.052) 

.133 

(.075) 

.087 

.057 

( .424, .424) 

.073 

(.068) 

.203 

( .119) 

.158 

.096 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling’s  t|  and  Mardia's 
statistic  using  the  F2  distribution  and  the  exact  null 
distribution  of  M,  respectively. 

Nominal  significance  level:  .04973. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 8,  n = 8. 

Distributions:  Class  given  in  6.2.1,  = 1^. 
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Table  7.11 

Two-Sample  Monte  Carlo  Results  for  Pearson  Type  II  and 
Type  VII  Distributions  with  m = 15,  n = 15,  and  a = .05318 


Bivariate 

median 

difference 

Statistic 

M 

t2 

MB2 

CS 

Bivariate  Normal 

(I  = I) 

(.000, .000) 

.058 

.088  (.063) 

.068 

.057 

( .01 4, .01 H) 

.072 

.113  (.083) 

.086 

.085 

(.028, .028) 

.122 

.1  90  (.1  46) 

.152 

.128 

( .042, .042) 

.224 

.350  (.277) 

.291 

.255 

Pearson  Type 

II 

(m  = 1 ) 

(.000, .000) 

.055 

.091  (.058) 

.065 

.058 

( .141  , .141  ) 

.062 

.111  (.086) 

.088 

.078 

(.283, .283) 

.096 

.180  (.1  43) 

.170 

.127 

( .424, .424) 

.151 

.340  (.269) 

.31  8 

.21  7 

Pearson  Type 

II 

(m  = .5) 

( .000, .000) 

.057 

.093  (.058) 

.069 

.080 

( .141  , .141  ) 

.062 

.109  (.087) 

.095 

.112 

(.283, .283) 

.094 

.182  (.1  42) 

.198 

.225 

( .424, .424) 

.136 

•334  (.266) 

.390 

.427 

Pearson  Type 

VII 

(m  = 2) 

( .000, .000) 

.061 

.066  (.044) 

.069 

.061 

(.141, .141) 

.082 

.090  (.055) 

.085 

.091 

(.283,-283) 

.172 

.148  (.116) 

.178 

.157 

( .424, .424) 

.335 

.252  (.198) 

.301 

.31  6 

Pearson  Type  VII 

(m  ■ 

= .5,  cauchy) 

( .000, .000) 

.051 

.029  (.01  9) 

.064 

.061 

( .1  41 , .1 41) 

.096 

.052  (.029) 

.093 

.103 

(.283, .283) 

.238 

.101  (.076) 

.185 

.237 

(.424, .424) 

.431 

.153  (.126) 

.332 

.432 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  Tp  using  the  F? 
distribution.  ’ ^ 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m = 15,  n = 15. 
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Table  7.11 — continued. 


Distributions: 

Bivariate  Pearson  Type  II 


|z|‘1/2{1-(y-0)'E_1(y-0)}rn  if  (y-0)  'I_1  (y-0)  < 1,  m > -1. 


Bivariate  Pearson  Type  VII 


Tim) 
T(m-1  )tt 


|r|_1/2{1 +(y-0)  (y-0) }“m,  m > 1,  y e R2. 


157 


Table  7.12 

Two-Sample  Monte  Carlo  Results  for  Bivariate  Normal  Mixtures 
with  m = 15,  n = 15,  and  a = .0531  8 


Bivariate  Statistic 

median 

difference  M t|  MB2  CS 


Mixture  1 (p  = 

.5,  »„ 

y12  ' 

-i.o, 

°11  °12  = 

1.0, 

P1  = 0.0,  y 21 

= u22  = 

t.o,  021 

°22 

= i.o,  p2  = 

0.0) 

( .000, .000) 

.053 

.076 

( .054) 

.059 

.065 

(.014, .014) 

.054 

.074 

(.057) 

.059 

.070 

( .028, .028) 

.066 

.105 

( .070) 

.080 

.089 

( .042, .042) 

.083 

.151 

( .107) 

.112 

.129 

Mixture  2 (p  = 

.5,  P„ 

= y1  2 = 

-.5,  o 

11  = °12  = 1 

.0, 

P1  = 0.0,  y21 

y22  = 

• 5,  o21 

°22 

= 1.0,  p2  = 

.9) 

( .000, .000) 

.045 

.068 

( .042) 

.057 

.048 

( .1  41 , .1  41) 

.055 

.069 

(.052) 

.062 

.053 

(.283, .283) 

.071 

.114 

( .085) 

.096 

.081 

(.424, .424) 

.126 

.1  94 

(.147) 

.159 

.140 

Mixture  3 (p  = 

.5,  P„ 

y12  = 

-1.0, 

°11  = °1 2 = 

1.0, 

P,  ■ -.9,  p21 

V22  = 

1.0,  02, 

°22 

= i.o,  p2  = 

0.0) 

( .000, .000) 

.051 

.071 

( .052) 

.057 

.049 

( .1  41 , .1 41) 

.058 

.078 

(.056) 

.065 

.049 

(.283, .283) 

.061 

.111 

(.077) 

.097 

.072 

(.424, .424) 

.086 

.166 

(.129) 

.156 

.114 

Mixture  4 (p  = 

.9, 

= y1  2 = 

0.0,  o 

11  ' °12  ' 1 

.0, 

P1  = 0.0,  y21  = 

y22  = 

0.0,  o21 

°22 

= 20.0,  p2  = 

0.0) 

( .000, .000) 

.051 

.021 

(.012) 

.060 

.042 

(.141, .141) 

.059 

.021 

(.011) 

.072 

.053 

( .283, .283) 

.091 

.035 

(.021) 

.108 

.103 

(.424, .424) 

.153 

.060 

(.033) 

.166 

.177 

Table  7.12 — continued. 


Bivariate 
median 
dif f erence 

Statistic 

M 

E- 

MB2 

CS 

Mixture  5 (p  = 

*9,  y11 

= v12  = 0.0, 

°11  = °12 

= 1.0, 

p1  = o.o,  u21 

y22  = 

0.0,  o = 20 

•o,  a22  - 

40.0, 

P2  = 0.0) 

( .000, .000) 

.053 

.021 

(.012) 

.060 

.041 

( .1  41 , .1  41) 

.063 

.021 

(.011) 

.071 

.054 

( .283,-283) 

.088 

.031 

( .020) 

.096 

.102 

(.424, .424) 

.130 

.049 

(.031) 

.147 

.179 

Entries:  Proportion  of  times  that  each  test  statistic  exceeded  its 

asymptotic  critical  value.  Values  in  parentheses  are 
rejection  proportions  for  Hotelling's  t|  using  the  F? 
distribution.  ’ 

Nominal  significance  level:  .05318. 

Number  of  samples:  1000  from  each  distribution. 

Sample  sizes:  m - 15,  n - 15. 

Distributions:  Bivariate  normal  mixture,  choosing  N1  with  probability 

p and  N2  with  probability  1-p,  where  N1 , i = 1,  2,  is 
bivariate  normal  with  null  mean  (y..,  y.5)  and  covariance 
matrix  ■ 11 


159 


In  conclusion,  the  two-sample  affine-invariant  procedures  perform 
similarly  to  the  corresponding  one-sample  procedures.  Mardia's 
statistic  performs  best  for  heavy-tailed  distributions  and  MB2 
performs  best  for  light-tailed  distributions.  Although  Mardia's 
statistic  does  not  perform  well  under  bivariate  normality,  MB2 

p 

performs  quite  well  relative  to  Hotelling's  T . Mardia's  statistic 

2 

performs  better  than  Hotelling's  7^  only  for  heavy-tailed 
distributions,  but  MB 2 performs  better  than  Hotelling's  both  for 
light-  and  heavy-tailed  distributions.  Also  MB2  continues  to  perform 
well  for  mixed  bivariate  normal  distributions  in  which  the  assumption 
of  elliptical  symmetry  is  violated.  As  is  true  for  the  one-sample 
statistic,  MB,  MB2  performs  well  for  a wide  variety  of 
distributions.  Hence,  the  use  of  MB2  in  the  bivariate  two-sample 
location  problem  is  well  supported. 
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